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Introduction 



It has long been known that the dynamics of a p-dimensional gravitational 
theory is captured by quantities on {p — l)-dimensional hypersurfaces [1]. It 
was argued by Damour [2J, based on an analogy by Hartle and Hawking [3], 
that in the case of certain black hole solutions these surface quantities describe 
the flow of a viscous {p — 1) -dimensional fluid. Perhaps the most concrete 
realization of this idea is the one to one map between large wavelength features 
of asymptotically AdSp black brane solutions and (p— l)-dimensional conformal 
fluid flows presented recently in Refs. |1] and [5]. 

Indeed, within the framework of string theory's AdS/CFT correspondence, 
an initially surprising relationship between the vacuum equations of Einstein 
gravity in an asymptotically locally AdS^+i space and the equations of hy- 
drodynamics in d dimensions has been unearthed. In particular a one to one 
correspondence has been found between, on the one hand, a class of regu- 
lar, long wavelength locally asymptotically AdS^+i black brane solutions to 
the vacuum Einstein equations with a negative cosmological constant and, on 
the other, all long-wavelength solutions of the d dimensional hydrodynamical 
equations 

V^T'^^ = 

of conformal fluid flows. The AdS/hydrodynamics correspondence provides an 
explicit black brane solution for every history of a particular conformal fluid so 
long as the fluid variables are constant over distances which are large compared 
with the inverse temperature. 

The nonrelativistic scaling limit - long distances, long times, low speeds and 
low amplitudes - of the relativistic equations of hydrodynamics, connects the 
relativistic conservation equation to the incompressible non-relativistic Navier- 
Stokes equations 

+ V.\7y = _Vp + uV^V + f , 

at 

that have been investigated for nearly two centuries but whose extremely rich 
dynamics still remainsl to be completely clarified. A particularly interesting 
phenomenon is turbulence. In fact, most fluid flows become turbulent under 
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a wide range of conditions. Even though turbulent flows are complicated phe- 
nomena of a statistical nature, it has been proposed that they could be actually 
governed by a new and simple universal mathematical structure analogous to 
a fixed point of the renormalization group flow equations {e.g. [6]). 

Given the fluid-gravity correspondence, it is of great interest to investigate 
the gravity dual solutions of turbulent flows, and the conditions under which 
turbulence may be expected. These turbulent solutions in the gravity dual 
may cast light on the possible turbulent decay of gravitational solutions near 
spacelike singularities - where indeed chaotic evolution is expected |7j. 

It is well known, for example from the Richardson's cascade model [S], that 
to realize steady state turbulence, one has to inject energy into a system. 

A first possible way to do this is to apply an external perturbation - a forcing 
function - deforming the fluid. Another way is to consider boundary conditions. 
In the context of the fluid-gravity correspondence, this first approach has been 
investigated in Ref. [H] where it was argued that a laminar fluid flow and the 
dual gravity solution could decay to turbulent configurations. The problem 
with this method is that forcing functions are necessarily inhomogeneous and 
bring rather complicated solutions. 

The second approach is closer to the intuition deriving from classical hy- 
drodynamics. Indeed many solutions of the Navier-Stokes equations describing 
fluids subject to hard wall type boundary conditions are known, even though 
it is still not clear how to generate gravitational duals of these boundary con- 
ditions. 

In this Thesis we study boundary conditions in the AdS/hydrodynamics 
correspondence as a preliminary investigation to the implementation of this 
second approach. In particular we study the gravity dual of a boundary layer 
separating two nonrelativistic, incompressible, fluid solutions with different 
properties: a stationary fluid on the left side and a moving fluid on the right. 
This simple setting is sufficient to produce turbulent motion in the fluid. For 
this configuration we derive gravitational duals of the boundary, following the 
prescription of Ref. While the fluid-gravity correspondence maps fluids on 
the left and on the right of the layer into two vacuum gravity solutions, the 
non trivial point is how to glue these two regions together. 

We consider two distinct methods, that have different limits of validity, 
with respect to the "thickness of the layer". The first one consists in gluing 
the two gravitational solutions by applying the Israel matching conditions [TT] 
to determine the stress tensor on the surface layer that separates the two 
sides. One can understand this approach as letting the gravitational solution 
continuously interpolate between the two solutions over a finite distance d and 
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then taking the hmit as this distance tends to zero while keeping the extrinsic 
curvature of the interpolating layer finite. 

The second possibility is to consider the gravity dual of a velocity function 
interpolating between the two fluid solutions. Clearly this function will not 
be a a solution of the Navier-Stokes equation in this region and therefore the 
dual will not be a solution of the vacuum Einstein equations in this region. 
Nevertheless one can follow this idea to sec where it leads. In this case, one 
cannot take the interpolation distance d to zero, because the dual is not defined 
when derivatives are large with respect to the inverse of the temperature T. 

A first interesting result that we will describe in detail in the last chapter 
is that both of these methods yield stress tensors that do not depend on the 
"thickness of the layer" d. Moreover the two methods yield bulk stress tensors 
which differ by a finite amount. In particular, we will sec that the disagreement 
between the two calculations of the stress tensor arises entirely from the higher 
derivative terms. Of course the fluid map is not deflned at small d, as it yields 
a divergent series, and so no divergences appear within the range of validity 
of either approach. 

The solutions that we find have to be handled with care. Indeed, as we 
shall point out, the stress tensors derived using the two different methods do 
not satisfy null energy condition and thus, as we will discuss, it is not clear 
whether such a boundary layer may exist. 

This Thesis is organised as follows. In Chapter 1 we shall give a short review 
of "classic" work on black holes of the seventies which led to the membrane 
paradigm, that is a picture of black holes as analogues to dissipative branes 
endowed with finite electrical resistivity, and finite surface viscosity. We shall 
review the derivation of the classical STirfacc viscosity of black holes, that has 
been re-derived in the quantum context of the AdS/CFT duality. We shall 
also provide an introduction to black hole thermodynamics and a (sketchy) 
derivation of the phenomenon of Hawking radiation, that of its crucial impor- 
tance in fixing the coefficient between the area of the horizon and black hole 
"entropy" . 

Chapter 2 contains a presentation of properties of relativistic fiuids and 
a discussion of the special case of conformally invariant fiuids that will be 
relevant in the construction of gravity dual solutions. We shall review the ex- 
tremely useful Weyl invariant formalism that we will apply in the perturbative 
study of dissipative fiuid solutions up to the second order in the derivative 
expansion. Finally, we will introduce the non-relativistic scaling limit that re- 
duces relativistic conservation equations of a conformal fiuid to the celebrated 
Navier-Stokes formula and we will study the residual conformal symmetry of 
this equation. 
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In Chapter 3 we shall introduce the basic scheme for constructing gravita- 
tional solutions dual to fluid flows. We shall present - in various forms - the 
AdS/hydrodynamic map up to the second order in the derivative expansion, 
and we shall discuss its conformal invariance. Then we shall turn to some of the 
physical properties of these solutions. We shall conclude the chapter with the 
discussion of the non-relativistic limit of the AdS/hydrodynamic map under 
the non-relativistic scaling limit introduced in Chapter 2. 

In Chapter 4 we shall present some techniques and issues that are relevant 
in the study of black holes and singular hypersurfaces in gravity. In partic- 
ular we shall start by reviewing various possible energy conditions that are 
relevant for black holes physics. Then we discuss the issue of the violation of 
the loosest among these conditions, the null energy condition, and its physical 
consequences. A relevant part of the chapter will be devoted to the introduc- 
tion of some geometric notions useful for the description of hypersurfaces and 
in particular the extrinsic curvature. These concepts will be used in applying 
the seminal work of Israel on singular hypersurfaces in General Relativity to 
our solutions. 

Chapter 5 contains a derivation and discussion of the main results reported 
in Ref. ^lOj. In particular we consider boundaries between nonrelativistic flows, 
applying the usual boundary conditions for viscous fluids. We find that a naive 
application of the correspondence to these boundaries yields a surface layer in 
the gravity theory whose stress tensor is not equal to that given by the Israel 
matching conditions. In particular, while neither stress tensor satisfies the null 
energy condition and both have nonvanishing momentum, only Israel's tensor 
has stress. The disagreement arises entirely from corrections to the metric 
due to multiple derivatives of the flow velocity, which violate Israel's finiteness 
assumption in the thin wall limit. 

Finally we summarise the findings and we discuss some open issues. An 
Appendix is devoted to the notation and conventions that we have adopted. 



Chapter 1 

Black holes as dissipative branes 



Solutions of general relativity describing spherically symmetric objects were 
obtained soon after the formulation of Einstein's theory. The singular be- 
haviour of these solutions - named black holes - was interpreted, following 
Subrahmanyan Chandrasekhar's proposal, as the description of a potential well 
due to a gravitational collapse of very massive astrophysical objects of mass 
M when their radius shrinks below the limit 2GM/c^, known as Schwarzschild 
radius. Early works assumed black holes to be passive objects, i.e. as given 
geometrical backgrounds. This viewpoint changed in the 1970's, when black 
holes started being considered as dynamical objects, able to exchange mass, 
angular momentum and charge with the external world. The study of the 
global dynamics of black holes was pioneered by Penrose [12] , Christodoulou 
and Ruffini [Ti] . Hawking [TB], and Bardeen, Carter and Hawking jTB] . 
Later on this approach evolved with the study of the local dynamics of black 
hole's horizons in what is called the "membrane paradigm" [2.1]) thanks to the 
works of Hartle and Hawking [3] , Hanni and Ruffini [1^ , Damour [21 UHl |T9] , 
and Znajek ^Oj. According to this point of view the horizon of a black hole is 
interpreted as a brane with dissipative properties described for instance by an 
electrical resistivity and a surface viscosity. 



1.1 Black hole solutions 

A few months after the publishing of the General Theory of Relativity, Karl 
Schwarzschild proposed a solution of Einstein's equations in the case of a spher- 
ically symmetric object with mass M, that can be regarded as the general rel- 
ativistic analog of the gravitational field of a mass point. In 3 + 1 dimensions, 
the metric of the Schwarzschild solution can be written as 



ds^ = -A{ry&t^ + 5(r)dr2 + (d^^ + sm^Odip^) (1.1) 
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where t denotes the time coordinate measured by a stationary clock at infinity, 
and where the coefficients A (r) and B (r) have the form 

A{r) = 1 - ^ , (1-2) 
Bir) = ^ . (1.3) 

This solution has an apparently singular behaviour at the so called Schwarzschild 
radius 

rs = ■ 1-4 

For example, the sphere located at r = rs has the observable characteristic of 
being an infinite-redshift surface. Indeed, the redshift of a clock at rest in the 
metric, whose ticks are read from infinity, via electromagnetic signals, 

dsatr=oo ^/-9oo{r = oo) 1 

— = , — = — , (1.5) 

dSatr V-^Oo(r) ,/l_2GM 

goes to infinity for r — )■ r^. Similarly the force needed to keep a particle at 
rest at a radius r > rs goes to infinity as r ^ r5. 

The 2-dimensional surface at r = r^, can be regarded as a 3-dimensional 
hypersurface T-L in spacetime. It is possible to see that the hypersurface is a 
fully regular submanifold of a locally regular spacetime via a change of coordi- 
nates near r = rg. Introducing the ingoing Eddington-Finkelstein coordinates 
(v, r, 6*, ip), where f = t + r*, with r*, the so-called tortoise coordinate, defined 

by 

/■ dr f dr 2GM ^ f ch- \ 



the line element takes the form 



ds^ = - (l - ^) dv^ + 2dvdr + r\de^ + sin^ Od^''] 



(1.7) 



In the new coordinates, the metric presents a manifestly regular geometry at 
r = rs- An important property of T-L is that it is a null hypersurface, i.e. a 
co-dimension- 1 surface locally tangent to the light cone. 

The vector normal to the hypersurface, i^, such that i^dx^ = for all 
infinitesimal displacements dx^j, on the hypersurface, is a null vector i.e. = 
0. Since £^ is a vector tangent to the hypersurface, we can consider its integral 
fines £^ — dx^/dt, which he within the horizon. These integral curves are called 
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the generators of the horizon. They are null geodesies curves, lying entirely 
within the horizon. The physical consequence of this fact is that, classically, T-t 
is a boundary between two distinct regions of spacetime: an external one from 
which light can escape to infinity and an internal one where light is trapped and 
out of which cannot escape. The internal region of is an infinite spacetime 
volume which ends, in its future, in a spacelike singularity at r = 0, where the 
curvature blows up as r~^. 

This kind of matter configuration is called black hole - since from the clas- 
sical point of view it absorbs all the light that hits it, just like a perfect black 
body in thermodynamics -, while the hypersurface T-L is called the (future) 
horizon. Black holes are conjectured to be the final, stationary state reached 
by any type of matter configuration undergoing a gravitational collapse. 

The Schwarzschild solution was generalized in independent works by Reiss- 
ner, and by Nordstrom, considering electrically charged spherically symmetric 
objects. In this case, the coefficients A (r) and B (r) in (11. ip are given by 
(setting for simplicity G = c = 1 here and in what follows) 

A{r=l + \, 1.8 

B(r) = -^^ . (1.9) 

It is easy to realise that now there exist two different horizons: an outer and 
an inner one, defined by 

r± = M ± ^M^ - Q2 ^ (1.10) 

which are the two roots of A(r) = 0. 

A more general black hole solution, due to Kerr and Newman, is obtained 
considering in addition to gravity {gnu) also electromagnetic interactions me- 
diated by long range fields {A^). The final, stationary configuration of such a 
black hole is described by three parameters, its total mass M, its total angular 
momentum J, and its total electric charge Q, and is given by the following 
Kerr-Newman solution 

ds^ = -| ujI + I dr^ + Sd^^ + ^ , (1.11) 
Qr 

- F^^ dx^ A dx" = ^{r'^ - cos^ 6*) dr A + ^ ar cos Q sin 6* d6' A a;^^ , 



where 



^ , A = - 2Mr + + Q^ S = + cos^ Q , 
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ojf = dt — a sin^ 9 dip , 



(r^ + a^) dip — adt . 



(1.12) 



The Kerr-Newman solution describes a black hole with regular horizon if 
and only if the three parameters M, J, Q satisfy the inequality 



A special class of black holes, called extremal, are those for which the inequality 
is saturated. A Schwarzschild black hole {a = Q = 0) can never be extremal, 
while a Reissner-Nordstrom black hole (a = 0) is extremal when \Q\ = M, and 
a Kerr black hole {Q = 0) is extremal when J = M^. 

1.2 Global dynamics of black holes 

Until the end of the 60 's black holes were seen purely as geometric backgrounds. 
A famous gedanken experiment was proposed by Penrose in 1969 [12] showing 
that in principle it would be possible to extract energy from a black hole. The 
idea of Penrose was to consider a charged test particle - with energy Ei, an- 
gular momentum p^^ , and electric charge Ci - coming from infinity and falling 
into a Kerr black hole, moving on generic nonradial orbits. By Noether's theo- 
rem, the time-translation, axial and U{1) gauge symmetries of the background 
guarantee the conservation oi E,p^ and e during the "fall" of the test particle. 
Now it is possible to imagine a process in which the test particle splits near 
the horizon of the black hole into two particles characterised respectively by 
energy, angular momentum and charge E2, p^p^i ^2, and £"3, p^g, 63. 

The black hole is described as a kind of gravitational soliton, that is as a 
physical object, localized "within the horizon "H", possessing a total mass M, 
a total angular momentum J and a total electric charge Q. Upon dropping 
a massive, charged test particles one expects a change in the values of M, 
J and Q. Therefore in the Penrose experiment the black hole should evolve, 
due to absorption of particle 3, from an initial Kerr-Newman black hole state 
(M, J, Q) to a final one (M + 6M, J + 5J,Q + 5Q) where 



Penrose found that , under certain conditions, one finds that particle 3 can be 
absorbed by the BH, and that particle 2 may come out at infinity with more 
energy than the incoming particle 1, indeed the process can lead to a decrease 
of the total mass of the black hole: 6M < 0. 



(1.13) 



SM = E3 = El- E2 
5J = J3 = Ji - J2, 
SQ = es = ei - 62- 



(1.14) 
(1.15) 
(1.16) 
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Figure 1.1: Penrose process: a particle "1" splits in two particles "2" and "3" 
while particle 2 escapes at oo, particle 3 is absorbed by the black hole. 



Christodoulou and Ruffini [131 [E] proposed a detailed analysis of the Pen- 
rose work, casting light on the existence of a fundamental irreversibility in 
black holes dynamics. Let us consider for the sake of simplicity a Reissner- 
Nordstrom black hole, and a process of the kind described in fll.l4p . Consid- 
ering an on-shell particle of mass /x, the Hamilton- Jacobi equation reads 

g'^'' {p, - eA^) {p, - eA,) = -/i^, (1.17) 

where a ( — h ++) signature is adopted. The momentum p^ is equal to 

p^ = dS/dx'' , (1.18) 

where S is the action that in an axisymmetric^ and time-independent back- 
ground can be taken as a linear function of t and cp 

S = -Et + p^^ + S{r,d). (1.19) 

E = —pt = —po is the conserved energy, p^p is the conserved y^-component of 
angular momentum while the last term contains contributions depending on 
the radial distance r and angular coordinate 9. 



black hole is said to be axisymmetric if there exists a one parameter group of isometries 
which correspond to rotations at infinity. 
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The Hamilton- Jacobi equation using eq. (I1.19P and the inverse metric for 
a Reissner-Nordstrom black hole reads 

" ^ - ^^o)' + + \(pe + -At^pI) = -/^' ' (1-20) 

A{r) \ sm 9 ^ ) 

that can be recast as 

bo - eA^f = A{r)'pl + A{r) L' + ^) • (1-21) 



Solving for E = —po, and inserting the expression for the electric potential 

= -- , (1.22) 

r 

one finds that the energy (where the sign discriminates between particles and 
antiparticles) is 



E=^± jA(r)2p2 + A{r) ( + ^ ] . (1.23) 



This expression generalises the formula for flat spacetime E = iy^/x^ + to 
a black hole background. 

It is possible to pin down the conserved energy of particle 3 when it crosses 
the horizon r+ of the Reissner-Nordstrom black hole (assuming that the con- 
dition Q < M for a regular horizon is fulfilled), i.e. the limit r — )■ r_|_ of £^3. 
The expression for E^ is 

E,= '-^ + \f\, (1.24) 
r+ 

where p'^ = g'^^Pr = A{r)pr has a finite limit on the horizon and the absolute 
value of p^ comes from the limit of a positive square-root. Now remembering 
that SM = Es and SQ = 63, the expression above can be re-written as 

and from the positivity of \p^\ follows the inequality 

6M > — (1.26) 
-r+(M,Q) ^ ' 

Inequality f ll.26p spells out the irreversibility of black holes energetics. Indeed, 
there can exist two types of process: reversible ones in which a particle of 
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charge +e with \p^\ =0 is absorbed for which the inequahty (I1.26P is saturated; 
and irreversible ones for which it is a strict inequahty. Reversible processes 
are clearly quite peculiar and difficult to obtain, therefore one expects that 
irreversibility will occur in most black holes processes. The problems of re- 
versibility of process in black holes physics replicate in some sense the relation 
between reversible and irreversible processes in thermodynamics. 

The calculation performed for a Reissner- Nordstrom black hole can be repli- 
cated in the more complex case of a Kerr-Newman black hole, yielding 

aSJ + r+QSQ + cos' 9 
— — = — — - — p\ 1.27 

r| + r+ -I- 

where as before a = J/M, and the bound Q' + {J/M^ < M' must be fulfilled. 
Once again we can write the variation in mass as an inequality 



6M > -^Jy^QJQ . (1.28) 



Processes for which this bound is saturated are called "reversible" because, 
after having produced a change SM, 6 J and 6Q, it is possible to perform 
a new process such that 6' J = —6J, 6'Q = —6Q (and the corresponding, 
saturated 6'M = — 5M)return the system to its initial state. On the contrary, 
any elementary process for which equation ( ll.28p holds as a strict inequahty 
cannot be reversed. 

Christodoulou and Ruffini considered a sequence of infinitesimal reversible 
changes in the state of the black hole, obtained by dropping in the black hole 
particles for which — )■ 0, and studied states which are reversibly connected 
to some initial state with defined mass M, angular momentum J and charge 
Q. Equation fll.28p simplifies to the partial differential equation for SM, 



5M = -^J + ^^QJQ ^ (1.29) 



which is integrable and that has solution in the Christodoulou-Ruffini mass 
formula 

The mass formula is composed of two terms: the square of the sum of the 
irreducible mass and of the Coulomb energy (oc Q), and the the rotational 
energy (oc J). Irreducible mass, defined as 



1 



Mi„ = -\/rl + a^ (1.31) 
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enters as an integration constant. To understand its role, it is useful to differ- 
entiate the mass expression and insert it in fll.28p . One obtains 

5Mi„ > , (1.32) 

where the equality holds only for reversible transformations, while the relation 
holds as strict inequality for all irreversible processes. 

The irreversible increase of the irreducible mass has a striking similarity to 
the second law of thermodynamics. Therefore we can interpret the quantity 
M — Mi„ as the free energy of the black hole, i.e. the maximum extractable 
energy by depleting (in a reversible) way J and Q is M — Mirr. Indeed the free 
energy has both Coulomb and rotational contributions (it vanishes in the case 
of a Schwarzschild black hole (J = = Q)). As a consequence of this relation, 
black holes can no longer be considered to be passive object since they actually 
store energy - up to 29 % of their mass as rotational energy, and up to 50 % 
as Coulomb energy - that can be extracted. 



There exists a link between the irreducible mass and the area of the horizon 
of a Kerr-Newman black hole. Indeed, the metric of a Kerr- Newman black hole, 
when using r = r+ (A = and dr = 0) for the inner geometry of the horizon, 
is 



dcr^ = 7AB(x'^)dx^dx^ 



.2 /^C\ a^.Aa^B 

(ri + cos^ 



9 9 9 ^N ,^9 sm U[rX + a , , ^ , ^ 

r\ + cos' g) dg' + . „ ^ + (1-33) 



Therefore the area of a timeslice of the horizon is proportional to irreducible 
mass 

Akn = y y (4 + sin ed9dip = 4:n{rl + a') = WttM^^ . (1.34) 



More generally. Hawking proved [15] a theorem stating that the area A of 
successive time sections of the horizon of a black hole cannot decrease 

5A>0. (1.35) 

Moreover the the sum of the area of a system of separated black holes also 
cannot decreases 




(1.36) 
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These properties are consequences of Einstein's equations, when assuming the 
weak energy condition. 

Hawking's result suggested the possibihty of a formulation of "thermody- 
namic laws" for black holes, in particular it suggests the possibility of defining 
an entropy for black holes proportional to the area of the horizon. 

While the area law is reminiscent of the second law of the thermodynamics, 
the analog of the first law can be obtained by varying the mass formula fll.30p 

dM (Q, J, A) 

where 

V 
Q 

9 

V can be interpreted as the electric potential of the black hole, and Q as its 
angular velocity. The coefficient g, called surface gravity, in the Kerr-Newman 
case, can be rewritten as 

and is zero for extremal black holes. The surface gravity of a Schwarzschild 
black hole reduces to the usual formula for the surface gravity of a massive 
star, g = GM/rl = M/{2Mf = 1/(4M). It is also possible to prove that the 
horizon has constant surface gravity for a stationary black hole, this result is 
known as the "zeroeth law of black holes thermodynamics"^. 

The expression (11.371) . once again, suggests the possibility of interpreting 
the term as some kind of entropy. 



^In the "membrane" approach to black-hole physics - that will be discussed later on -, 
the uniformity of g for stationary black hole states can be viewed as a consequence of the 
Navier-Stokes equation since the "surface pressure" of a black hole happens to be equal to 
P = s/Stt. 



VdQ + ndJ + —dA 
Stt 



:i.37) 



a 

1 r+ — r_ 
2rl + a? 



:i.38) 
:i.39) 
:i.40) 
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1.3 Black hole thermodynamics 

In 1974, Bekenstein tried to push the thermodynamical analogy for black hole 
physics further by taking into account quantum effects. In particular he pro- 
posed Carnot-cycle-type arguments: it is possible in principle to extract work 
from a black hole by slowly lowering into it an infinitesimal a box of radiation, 
in this way all the energy of the box of radiation, mc^, would be converted 
into work. The efficiency of this Carnot cycle is defined as 

^ = 1-^, (1.42) 

J- hot 

where Tcom and Thot are respectively the two source temperatures between 
which the "engine" operates. Classically, Tcom is expected to be zero and 
therefore rj is expected to be 1. Bekenstein observed that quantum effects 
should arise modifying the classical result. Indeed the uncertainty principle 
poses restriction to the size of box of thermal radiation at temperature T: since 
the typical wavelength of radiation is A ~ the box will have a minimum 
finite size ~ A. From this limit on this size of the box, Bekenstein derived an 
upper bound a bound on the efficiency rj and therefore concluded that there 
will exist a nonzero temperature of the black hole Tbh 7^ 0. 

A second argument proposed by Bekenstein consider a reversible process 
in which a particle with p"^ = hits the horizon of a black hole at r = r_|_. 
The absorption of such a particle should not increase the surface area of the 
black hole. However, for this to be true both the (radial) position and the 
(radial) momentum of the particle must be exactly fixed: namely, r — r+ and 

= 0. This is in contradiction with limits posed by Heisenberg's uncertainty 
principle. To state this argument formally, we need the covariant component 
Pr of the radial momentum 

p'- = g'-'-p^ = A{r)pr 

= t^lAt^p^ ^ S^I±^Pr (1.43) 
r r+ 

- (^-^J ^rpr. (1.44) 

It is useful to notice that the surface gravity g is proportional to the partial 
derivative of A with respect to r. 
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Therefore we can re-write the expression for p*" as 

~ 2g5r5pr. (1.46) 

From the uncertainty relation 

SrSpr >lh, (1.47) 

we get 

p'>gh, (1.48) 

Substituting the bound on p'^ in the expression for 6M, we find the following 
inequality 

SM-^= l/l > gh, (1.49) 

r+ 

which using expression fll.34p for irreducible mass, can be written as 

5A > 8nh. (1.50) 

This result show that there is a strong limitation arising from the quantum 
level on the possibility of reversible processes. A particle falling into a black 
hole always produces an irreversible process that increases the area of the black 
hole by a quantity of order h. Such a process represents a loss of information - 
the information about the particle - for the world outside the horizon. For this 
reason, Bekenstein suggested the introduction of an entropy for black holes 
j22] of the form 

SB}i = a—A, (1.51) 

whose coefficient - which Bekenstein was not able to fixed in a unique, and 
convincing, manner - had to be dimensionless and of order one, d ~ (9 (1) 
(Bekenstein proposed a = ln2/87r). 

An important consequence of the definition of an entropy for black holes is 
the natural definition of a temperature 

Tbh = ^-g. (1.52) 

The existence of a temperature implies the possibility that black holes may 
emit radiation, contrasting the assumption about their "black" behaviour. In- 
deed, soon after Bekenstein's original suggestions, Stephen Hawking - who was 
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skeptical about the possibility that black holes could radiate - in 1974 discov- 
ered the universal phenomenon of quantum radiance of black holes [23j and 
was capable to fix in an unambiguous way the coefficient a to be 

a = \. (1.53) 
We shall discuss the Hawking computation in detail in the next section. 



1.4 Black hole quantum radiation 

The surprising phenomenon of Hawking radiation, that is the thermal radia- 
tion with a black body spectrum predicted to be emitted by black holes due 
to quantum effects, was first proposed by Hawking in 1974 [23]. In this sec- 
tion we will follow a simplified derivation due to Damour and Ruffini pl| [25]. 
considering for simplicity a 3 + 1-dimensional Schwarzschild black hole. 

Let us consider a massless scalar field v?(x) propagating in a Schwarzschild 
background. The Klein-Gordon equation coupled to gravity spells out 

= ^ d^{^gg^- if) = 0. (1.54) 

The solutions of this equation - given the symmetries of the background - 
can be decomposed into mode functions, given by the product of a Fourier 
decomposition into frequencies, spherical harmonics and a radial dependent 
factor, namely 

Vu^nmit, r, e, if) = -^=== "'^^"'^^^ Y(,^{e, if) . (1.55) 
V27r|a;| r 

Introducing a "tortoise radial" coordinate r^,, substituting the explicit form of 
the metric in the new coordinates and using the generic solution (11.551) . the 
Klein-Gordon equation boils down to a radial equation for Uuiimi^)'- 

— + (a;2-V,[r(r.)])« = 0, (1.56) 

where the effective radial potential has the form 
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Figure 1.2: Representation of the effective gravitational potential Ve{r) in the 
neighbourhood of a black hole. Spacetime is essentially flat both at infinity and 
near the horizon. In the central region the tidal-centrifugal barrier produce a 
grey body factor. 

The effective potential Vi{r) vanishes both at r — t- oo (which corresponds to 
r=K — )■ +oo) as the massless centrifugal potential l)/r'^, and at the horizon 
r = 2M (r* — > +oo). Therefore, the potential is negligible in these two regions 
and the solution of the wave equation is expected to behave essentially as in 
fiat space. The general solutions (p in the asymptotic region will be 



Conversely, the potential generated by the gravitational coupling is effective 
only in the intermediate region where it represents a positive potential bar- 
rier combining the effect of curvature and centrifugal effects. Modes generated 
near the horizon must penetrate this barrier to escape to infinity, hence a grey 
body factor which diminishes the amplitude of the quantum modes will appear 
in the solution of the Klein-Gordon equation. 

To quantize the scalar field if we will follow a rather standard procedure 
decomposing the field in the eigenf unctions of the Klein Gordon equation, 
with coefficients given by creation and annihilation operators. In the case of a 
black hole background the definition of positive and negative frequencies must 
be handled with care. 

Let us start reviewing the general formalism in the case of a quantum 
operator ip{x) describing real massless particles in a background spacetime 




(1.58) 
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which becomes stationary both in the infinite past, and in the infinite future. 
The operator (^(,x) can be decomposed both with respect to some "in" basis 
of modes, describing free, incoming particles, 

0ix) = «i>r(^) + , (1-59) 

i 

and with respect to an "out" basis of modes, describing outgoing particles, 

i 

where {d"\ (a™) + } and {a°"*, (a°"*)+} are the two sets of annihilation and 
creation operators, such that 

(«?)^] = ■ (1-61) 

The two sets of annihilation and creation operators correspond to a decom- 
position in modes which can physically be considered as incoming positive- 
frequency ones incoming negative-frequency ones which can 
be taken to be the complex conjugate oi pf^{x) in our case, outgoing positive- 
frequency ones and outgoing negative-frequency ones Mode 
functions can be normalized as 

{pTJD = , (pr, ri?) = , (rC, <) = -5,j , (1.62) 
where ( , ) denotes the standard (conserved) Klein-Gordon scalar product 

((/?!, (^2) ~« J da^ {(pld^ip2 - d^(p\(f2) . (1.63) 

Vacuum states |in), |out) are defined as respectively the states annihilated by 
a" and a°^*. The mean number of i-type "out" particles present in the vacuum 
I in) is given by 

(iV,) = (in| (ar)+ ar |in) = \T^3? (1-64) 

where = {pf^^in™) is the transition amplitude (also named Bogoliubov 
coefficients), between the incoming negative- frequency mode n™ and the out- 
going positive-frequency one 

In the case of a black hole background, the application of the general for- 
malism, detailed above, can be problematic since the background is not asymp- 
totically stationary in the infinite future - in the interior of the black hole the 
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Killing vector d/dt is spacelike -, while it can be considered asymptotically 
stationary in the infinite past only by ignoring the formation of the black hole. 
In his work, Hawking pointed out that it is possible overcome these problems 
by considering the high-frequency modes coming from the infinite past and 
the outgoing modes, viewed in the asymptotically flat region and in the far 
future. Outgoing modes can be unambiguously decomposed into modes with 
positive and negative frequency, because, as explained above, their asymptotic 
behaviour is given by a sum of essentially flat-spacetime modes as seen in 
(USE]). 

We want to calculate the average number of outgoing particles seen in the 
I in) vacuum, that is the transition amplitude J2j l^iiP; defined in (11.641) . from 
an initial negative frequency mode n*"" into a final outgoing positive frequency 
one observed at infinity. To compute the sum |Tjjp, we have to find 
out what is an initial negative frequency mode nj^. As pointed out before, 
there is a physically infinite redshift between the surface of the horizon and 
asymptotically fiat space at infinity. Therefore, particle with finite frequency 
at infinity will correspond to a wave packets with very high frequency near 
the horizon, hence very localized. A second important observation is that the 
near horizon geometry is regular and with a finite radius of curvature, hence 
it is sensible to approximate it locally by a fiat spacetime. In conclusion, the 
technical issue that we need to discuss is how to describe a negative frequency 
mode nj^ in a small neighbourhood of the horizon, that we can assume to be 
a Minkowski fiat vacuum. 

Let us introduce the technical criterion, that we will use later on, for char- 
acterizing positive and negative frequency modes in a locally fiat space-time. 
Consider, in Minkowski space, a wave packet 

^^{x) = (fk^{k)e'^^'"^ (1.65) 

made only of negative frequencies, i.e. such that the 4-momenta k^ in the 
wave packet are all contained in the past light cone C~ oi k^. A convenient 
technical criterion for characterizing, in x-space, such a negative-frequency 
wave packet is the well-known condition that ^-{x) be analytically continu- 
able to complexified spacetime points + iy^ with lying in the future light 
cone, G . If one performs a complex shift of the spacetime coordinate, 
— )■ + iy^, where y^ is timelike-or-null and future-directed [y^ G C"*"), 
then, the e*'^''^'' term will be suppressed by a e'^'^^'" term, where the scalar 
product k^y^^ is positive because it involves two timelike vectors that point in 
opposite directions (we use the "mostly plus" signature). This ensures that a 
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negative-frequency wave-packet can indeed be analytically continued to com- 
plex spacetime points of the form + iy^, with G C"*". 

Now we are ready to attack the calculation of the transition amplitude in 
the Schwarzschild background. The first step is to eliminate the "singularity" 
at the horizon r = 2M by adopting the Eddington-Finkelstein coordinates 
introduced above fll.7l) . Near the horizon, an outgoing mode from 7i, which 
will have a positive frequency at infinity, is given by expression f ll.SSp with a 
minus sign in front of r*, and with a; > 0. Introducing Eddington-Finkelstein 
coordinates, this mode function can be re-written as 

[^T{v, r)Urn oc 6+^^ ( ) , (1.66) 



where we used 



t-r^ = t + n-2r^ = v-2r^ = v-2r-AM\ni j . (1.67) 

The outgoing mode fll.66p appears to be singular on the horizon where oscilla- 
tions have shorter and shorter wavelengths as r — )■ 2M, and not to be defined 
inside the horizon. 

The local, x-space criterion for characterizing positive and negative fre- 
quency modes, discussed above, can be applied - in a local frame near the hori- 
zon - to the wave packet (11.661) . Since the infinitesimal displacement r — t- r — e, 
f — )■ f is seen to be future directed and null, this criterion finally tells us that 
the following new, extended wave packet, defined by applying the analytic 
continuation r ^ r — ie to (11.661) . 

fr -2M - jrV^^^"^ 
n^{v, r) ^ (^', r-ie)^ e'^^^' e^-'^ ( — j (1.68) 

is, when Fourier analyzed in the vicinity of "H, a negative frequency wave 
packet. Note that the new wave packet n^{v,r) is defined also in the interior 
of the black hole and that we have included a new normalizing factor in its 
definition in terms of the analytic continuation of the "old" mode v^J^"*, which 
had its own normalization. 

The normalisation factor A^^ is relevant for the computation of the Hawking 
radiation. Initially modes in (I1.58P were normalised as 

(V^wi^imi, V^a;2^2m2) = +5(l^l — 1^2) ^£1(2 Smim2- (1.69) 

The analytic continuation r ^ r — ie introduces, via the rotation by e"*'^ of 
r-2M in (r-2M)*^*^^, a factor (e-'^Y^^^"^ = 6+^^^^ in the left part (r < 2M) 
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of n^, that is 

na;(r) = Ar^[^(r-2M)y^°'^*(r-2M)+e^"^-'"^(2M-r)^°"*(2M-r)], (1.70) 

where 9{x) is the Heaviside step function. The interpretation of eq. fll.70p is 
that the initial negative-frequency mode rii^, locahsed near the horizon, over 
time, sphts into an outgoing mode <y9°"*(r — 2M), that is observed as a positive- 
frequency mode at infinity, and a second mode y9°^*(2M — r) that is absorbed 
by the black hole. 

Now the calculation of the scalar product (rit^i^imn '"'W2fem2) brings a factor 
INujI"^ [1 — {e^'^^'^'^Y], where the minus sign is due to the essentially negative- 
frequency aspect of <^Ji"*(2M — r) 6{2M — r). Therefore the appropriate nor- 
malisation factor to get 

('^tJi^imi, '"'(.J2^2m2) = — 5(l^l — 1^2) ^£i£2 ^mim2 {^■'^^) 

for a negative-frequency mode is 



l^^l' = ,sJ. _ 1 • (1-72) 



Leaving aside, for the moment the effect of potential Ve{r), the transition 

f = n^.hm, and oc y^^f^ - 



amplitude Tjj - the scalar product between = n^i^^mi and oc ip°^^ 



would be 

{'^°L!drn-,'^^ihrn-^) = N^5{uJ - UJi) Smm^ ■ (1-73) 

Therefore, the number of created particles fll.64p will contain \N^^\^ times the 
square of 5{(jJ — Wi), which, by Fermi's Golden Rule, is simply 5{uj — ui) x 
/ rft/27r. 

To take into account for the gravitational potential Vi{r)., it is necessary 
to add the grey body factor r^(a;), giving the fraction of the flux of which 
ends up at infinity because of the effect of Vn{r). The final result for Hawking 
radiation in a Schwarzschild black hole background is 

A black hole radiates as if it were a black body of temperature Tbh = 1/ (SttM), 
screened by a gray body factor Ti^co). 

From the Planck factor in (ll.74p . it is possible to find the Hawking tem- 
perature, 

T = h^. (1.75) 
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i. Therefore the Bekenstein 



(1.76) 



The extension of the analysis proposed to more general classes of black 
holes backgrounds leads essentially to the substitution of the factor 4M by the 
inverse of the surface gravity of the black hole. The Planck spectrum contains 
a factor (^e27r(a;~p^n-e<i')/re _ guiximarising the combined effect of the general 
temperature and of the couplings of the conserved angular momentum and 
electric charge e to, respectively, the angular velocity VL and electric potential 
$ of the black hole. 

Hawking's radiation is not astrophysically relevant for stellar-mass or larger 
black holes, nevertheless it has been observed [21] that the combined effect of 

2'k{uj — ujq) 

the grey body factor and of the zero-temperature limit of (e « — 1)^^ 
could yield potentially relevant particle creation phenomena in Kerr-Newman 
BHs, associated to the "superradiance" of modes with frequencies fi < u < Uq, 
where /i is the mass of the created particle. 

1.5 Surface electrodynamic properties 

The modern study of the local dynamics of black hole horizons originates from 
a "holographic" description of a black hole's properties: all of the physics of a 
black hole is condensed in the description of a set of surface quantities on the 
horizon - the surface of the black hole - and a set of bulk properties outside of 
the horizon. This approach is called the "membrane paradigm" [2T]). 

A first interesting example of this holographic description pertaining to 
the electromagnetic properties of charged black holes. It is possible to modify 
the field equations of the electromagnetic field F^i, = d^Ay — d^A^ - that 
in principle is expected to permeate the whole spacetime - so as to replace 
the internal electrodynamics of the black hole by surface effects. Maxwell's 
equations 



can be modified by introducing a fictitious field F^jy(x)0H, where 9^ is a 
Heaviside-like step function, equal to 1 outside the BH and inside. The new 



V^J" = 0, 



(1.77) 
(1.78) 
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equation for F^,^ is of the form 

(F^-'e) = {v^F^"") e + F^-'v^e (1.79) 

where it is possible to recognise two source terms, that we can re-interpret 
introducing a surface current defined as 

fn = ^^'^'V.e , (1.80) 

where the term V,^© contains a Dirac 5-function which restricts it to the 
horizon. If we consider a generic scalar function (/? (x) such that (/? (x) = on 
the horizon, with (x) < inside the horizon, and (p (x) > outside it, by the 
properties of the Heaviside step function we get 

O^^Qh = d^e {ip {x)) = 5 {ip {x)) d^ip , (1.81) 

where 5 is the a Dirac delta function with support on the horizon. The modified 
field equation, with the new source term, now reads 

{F'"'Q) = 47r ( J^e + j^) . (1.82) 



To complete the description of black hole electrodynamics it is necessary to 
modify the current conservation equations. In particular, it is useful to define 
a surface current density on the black hole. For this purpose it is important 
to stress that there is an important subtlety in the definition of the normal 
vector to the horizon, due to the fact that, as discussed before, the horizon 
is a null hypersurface. Indeed, the covariant vector such that i^dx^ for 
any infinitesimal displacement dx'^ on the horizon has norm zero ^^^'^ = 
and therefore cannot be normalised in the usual way adopted in the Euclidian 
space. In stationary axisymmetric spacetimes, £^ can be uniquely normalised 
by demanding that the corresponding directional gradient l^d^ be of the form 

with a coefficient one in front of the time-derivative term. For this reason it is 
possible assume that £^ is normalized in such a way to that its normalization is 
compatible with the usual normalization when considering the limiting case of 
stationary-axisymmetric spacetimes. Therefore given any normalization, there 
exists a scalar ou such that 

= cod^V, (1.84) 
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hence it is possible to define a 6 function on the horizon 

SH = -6{ip), (1.85) 

such that 

O^Oh = i,5H. (1.86) 



At this point it is possible to re-write the surface current in the form 

j^^ = K^5^, (1.87) 
where we defined a black hole surface current density as 

Kf" = -^F^X , (1.88) 

making evident that the surface current is due to the presence of the electro- 
magnetic field on the horizon. The conservation of electric current now spells 
out 

V,{QhJ'' + K^^5h) = Q, (1.89) 

which is manifestly the conservation of the sum of the bulk current QrJ^ and 
of the boundary current K^6h- In this way the description of electromagnetic 
phenomena for black holes has be rephrased in terms of surface quantities de- 
fined locally on the horizon. 

It is useful to introduce some more technical tools to describe the surface 
physics of black holes. Performing a change of coordinates using (advanced) 
Eddington-Finkelstein-like time coordinates for which t = x^, is equal to 
zero on the horizon (like r — r_|_ in the Kerr-Newman case), and x^ for A = 2,3 
denotes some angular-like coordinates. In the new coordinates the normalisa- 
tion of the normal vector to he horizon takes the form 

Expression ( ll.90p suggests the interpretation of as the velocity of some 
"fiuid particles" on the horizon, which can be seen as the "constituents" of a 
null membrane. Following this analogy, we can - as is usually done in the study 
of fiuids - introduce the gradient of the velocity field, splitting it into its sym- 
metric and antisymmetric parts, where the antisymmetric part is interpreted 



1.5. Surface electrodynamic properties 



25 



as a local rotation and has no consequence on the physics. The symmetric 
part can be further divided into its trace and trace-free parts, namely 

^ {diVj + djVi) = Gij + ^9 • v5ij (1.91) 

where d is the spatial dimension of the considered fluid - which will he d = 2 
for black holes in 4-dimensions. The usual description interprets the first term 
as the shear, and the second as the rate of expansion of the fluid. Analogous 
quantities will be defined for black holes. 

The local spacetime geometry on the horizon is described by the degenerate 
metric 

ds\i=Q = 'JAB (t, x^) (dx^ - vMt) {dx^ - v^dt) , (1.92) 

where = dx'^/dt. ■jAB{t,x) is a symmetric rank 2 tensor, i.e. a time- 
dependent 2-metric such that the horizon may by viewed as a 2-dimensional 
brane. The element of area of a spatial sections St therefore can be expressed 
as 

dA = y/det-fABdx^ A dx^. (1.93) 

One can decompose the current density into a time component 

and two spatial components tangent to the spatial slices St {t — const.) of 

the horizon, 

K'^d^^andt + K^dA (1.94) 
in which dt — I'^d^ — v'^Oa so that 

K^d^ = ant + {K^ - aHV^)dA. (1.95) 
The total electric charge of the spacetime can be defined as the integral 

Qtot = ^/ ^F^'^d^^. , (1.96) 

which can be rewritten, using Gauss' theorem, as the sum of a surface integral 
on the horizon and a volume integral in between the horizon and oo - that can 
be seen as the charge Qh contained in space - of the form 

The surface element can be expressed in terms of the null vector orthogonal 
to the horizon £^ and of a new null vector n^^ transverse to the horizon and 
orthogonal to the spatial sections St, that is normalised as n^dn = +1, 

dS^^ = ^^e^^pAx" A dx" = (n^4 - nj^) dA . (1.98) 
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Using the definition of the surface current, the black hole's charge is written 
as 

Qh = ^ (TRdA, (1.99) 
where an was introduced before as the time component of the surface current. 

The above definition of the charge for a black hole forces one to consider 
the density an as a charge distribution on the horizon. Therefore one can 
interpret 

an = K^'n^ = ^F^^'n^l, (1.100) 
as the analog of the expression 

a = —E'rii (1.101) 

giving the electric charge distribution on a metallic object. Extending the 
reasoning by analogy to the spatial currents fiowing along the surface one can 
rewrite the current conservation equation fll.82p as 

^ ^ (v^^h) + (V^i^^) = (1.102) 



^7 dt ^7 dx^ 

This form renders manifest the role of the surface current in "closing" the ex- 
ternal current injected "normally" into the horizon in combination with an an 
increase in the local horizon charge density. 

The electromagnetic 2-form restricted to the horizon takes the form 

^F^^dx^ A dx^ln = F^da;^ Adt + B^dA . (1.103) 

Taking the exterior derivative of the left-hand side of the above expression one 

gets 

V X E = -^dti,/^B^_). (1.104) 

which relates the electric and magnetic fields on the horizon. 

An analog of the Ohm's law relating the electric field to the current can be 
written in the form 

Ea + eABB^v"" = An-fAB {K"" - aHv"") , (1.105) 
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that can be recast as 

E + V X = 4:71 (^K - (ThvY (1.106) 

From this expression it is natural to define the surface electrical resistivity for 
black holes to be equal to p = 47r = 377 Ohm ^Si 



1.6 Black hole surface viscosity 

The holographic approach has a second interesting application to Einstein's 
equation on the surfaces of black holes. Indeed defining appropriate surface 
quantities related to the spacetime connection one is able to obtain from the 
equations of gravity a "surface hydrodynamics" described by a Navier-Stokes- 
like equation. The path followed in defining surface electrodynamical quanti- 
ties must be amended due to the non- linear nature of Einstein's equations^. 

To address the issue it is necessary to make a few more technical remarks 
about the near horizon geometry. Given any hypersurface, the parallel trans- 
port along some tangent direction, which we denote by t, of the (normalized) 
vector i normal to the hypersurface yields another tangent vector. In general 
i ■ i = e, where e will be —1 for a time-like hypersurface, +1 for a space-like 
one and zero for a null hypersurface. The directional gradient of the norm of 
the vector i along the arbitrary tangent vector t will give 

(V^-f) -1=0 . (1.107) 

Therefore the vector (V^-f) must be tangent to the hypersurface. More gener- 
ally there exists a linear map K, acting in the tangent plane to the hypersur- 
face, such that Vf£ = K{t). In the case of time-like or space-like hypersurfaces, 
this map is called Weingarten map and is given by the mixed- component Kj 
version of the extrinsic curvature of the hypersurface Kij. The case of a null 
hypersurface is more involved since the extrinsic curvature Kij is not uniquely 
defined. In any case, there exists a mixed-component tensor K^, intrinsically 

defined as the Weingarten map K in V/= ir(t). 

In the coordinate system we defined above ^ {A = 2,3), where 

the horizon is located at xi = 0, a basis of vectors tangent to the horizon 
can be defined containing the null vector i, and the two space-like vectors 



■^For a detailed derivation of the results reported in this section one can refer to [H [TOl 
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= <9a. With this basis, the Weingarten map K is fully described by the set 
of equations 

V^i=gl (1.108) 
V^f = VLAi+ D^cb- (1.109) 

The first equation spells the fact that £ is tangent to a null geodesic lying 
within the horizon, and can be seen as a very general definition of the surface 
gravity g. In the second equation, a two-vector Qa, and the mixed component 
of a symmetric two-tensor Dab, which measures the "deformation" in time 
of the geometry of the horizon, are present. Tensor Dab is defined as 

D,s = -,BcD-, = l^ (1.110) 

where where D/dt denotes the Lie derivative along i = dt + v^Oa of the 
degenerate metric 'Jbc defined in (ll.92p . In explicit form the deformation 
tensor can be written as 

Dab = ^ {dtlAB + v'^dc-fAB + dAv'^-fcB + dBv'^-fAc) (I-IH) 

= ^ {dtlAB + va\b + vb\a) ■ (1.112) 

Here the symbol '|' indicates a covariant derivative with respect to the Christof- 
fel symbols of the 2-geometry 'Jab- The tensor is made of two parts, the ordi- 
nary time derivative of 7ab, and that from the variation of the generators of 
velocity along the horizon. To interpret this object, it is useful to split the 
deformation tensor Dab into a trace-less part and a trace. 

Dab = ctab + ,1ab (1.113) 
where the traceless part (Tab is interpreted as a "shear tensor" , while the trace 

e = Di=^-^^''daAB + vfA, (1.114) 

is an "expansion" term. The remaining component of the Weingarten map, 
namely the 2- vector Qa, is defined as 

nA = n-VA£ (1.115) 

with i ■ n = 1. To give this component a physical meaning it is useful to 
introduce an angular momentum for the black hole. 
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Let us consider, as before, an axisymmetric spacetime. There will exist a 
Killing vector related to this symmetry, 

to which the Noether's theorem will associate a conserved total angular mo- 
mentum, which can be written as a surface integral over the 2-sphere, Soo 



Joo = -^ [ Wk'^dS,,, (1.117) 



where the surface element dS^u was defined above in the case of the charge 
distribution. As before, one can transform the surface integral into the sum 
of a volume integral containing the angular momentum of the matter present 
outside the horizon and a surface integral over the horizon 

J = -/matter + Jh ■ (1.118) 

The second term can be interpreted as the black hole angular momentum. 
Using the expression (11.981) for the surface element in terms of the vectors 
and and observing that we can perform the exchange fWyk^^ = k'^Vi^i'^ 
since the two vectors have zero commutator 

[£, k]=0 , (1.119) 

the black hole angular momentum Jh takes the integral form 



Jh = -3^ / n^FV^rdA. (1.120) 

'Sh 



8n 



We can re- write the above expression as the projection of on to the direction 
of the rotational Killing vector k = d^, so that we have 

Jh = — - (p k^QAdA (1-121) 
Stt Js 

where k^flA is the yj-component of Qa- It is natural at this point to define a 
"surface density of linear momentum" as 

TTA = -^^A = -^n ■ Va/. (1.122) 
Finally we get for the total angular momentum of the black hole the expression 

Jh = [ n^dA . (1.123) 
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It is possible to find a dynamical equation for the local quantities defined 
above by contracting Einstein's equations with the normal to the horizon. In 
this way one relates surface quantities to the flux of the energy-momentum 
tensor T^j, into the horizon. By projecting Einstein's equations along i'^e'^, 
one finds 

where 



Dtta 
dt 



{dt + 9)7lA + v''7lAiB + vfATrB, (1.125) 



(TAB = ^ {dtlAB + Va\B + Vb\a) " ^^lAB, (1.126) 



+ < (1-127) 



correspond to a convective derivative, a shear and an expansion rate respec- 
tively. 

The equation found can be compared to the Navier-Stokes equation for a 
viscous fluid 

(dt + d)7ri + v'^TTi^k = + 2r7<, + CO,i + fi, (1.128) 

where tTj is the momentum density, p the pressure, t] the shear viscosity, 
Cij — I {vij + Vj^i) — Trace, the shear tensor, ( the bulk viscosity, 9 = f *j 
the expansion rate, and fi the external force density. The striking analogy 
between these two equations suggests the description of a black hole's sur- 
face as a brane with positive surface pressure p = +-^, external force-density 
Ja = —C'^T^A which corresponds to the flow of external linear momentum, 
surface shear viscosity r] = and surface bulk viscosity ( — 

It is important to remark the non-relativistic character of the black hole 
hydrodynamical-like equations, a quite surprising feature, in spite of the "ultra- 
relativistic" nature of black holes. 



1.7 Local thermodynamics of black holes 

In previous sections, starting from Maxwell's equations and Einstein's equa- 
tions, a surface resistivity and a surface shear viscosity of black holes were 
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defined. On the same grounds, the existence of a local version of the second 
law of thermodynamics for black holes is of particular interest. Indeed, physi- 
cal systems verifying Ohm's law and the Navier-Stokes equation are expected 
to be also endowed with thermodynamic dissipative equations (the equivalent 
to Joule's law). Naively, the expected "heat production rate" in each surface 
element dA would be of the form 



q = dA 



;i.l29) 



where p is the surface resistivity, and t] and ( the shear and bulk viscosities 
defined above. The corresponding local increase of local entropy s = adA can 
be found to be 

? = I- (1.130) 
dt T ^ ' 

where T = is the local temperature on the surface of the black hole^. 
A precise result can be obtained by contracting Einstein's equations with 
The projection gives the equation 



ds d^s _ dA 
dt~^d^ ~ T 



;i.l31) 



The similarity of the evolution law for the entropy found with eq. fll.l29p is 
quite striking, but there is a relevant difference due to the unexpected second 
derivative of local entropy on the left hand side appearing with a minus sign 
and a coefficient r = g~^, that can be regarded as a time scale. This term is a 
correction to usual near-equilibrium thermodynamics, which involves only the 
first order time derivative of the entropy^. 

In the limit of an adiabatically slow evolution of the black-hole state, eq. 
f ll.lSip reduces to the usual thermodynamical law giving the local increase of 
the entropy of a fiuid element heated by the dissipation associated to viscosity 
and the Joule's law. 

In the approximation of constant r, the rate of increase of entropy is 



di J, T \T 

In other words, the rate of increase of entropy is defined as integral of the 
heat dissipation over the future. This fact points to a very peculiar feature 



'*Even though, as consequence of the zeroth law of black holes thermodynamics, the 
surface gravity is uniform on the horizon of stationary black holes, it is, in general, non- 
uniform for evolving black holes. 

^An interesting observation is that for a — 1/4 one gets r — that corresponds to 
the inverse of the lowest "Matsubara frequency" associated to the temperature T. 
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of black holes: their acausal nature. Indeed, black holes are defined as null 
hypersurfaces which are evolving towards a stationary state in the far future. 

An important value for the black holes physics is the ratio of the shear 
viscosity r] = l/(167r) to the entropy density s = s/dA = a, that for the 
Bekenstein- Hawking value of d = 1/4, is 



V 



1 



1 



(1.133) 



s 



alQn 



This value of the ratio rj/s has been found Kovtun, Son and Starinets in the 
gravity duals of strongly coupled gauge theories, using the AdS/CFT corre- 
spondence f2E\f2^. 



Chapter 2 

Elements of fluid dynamics 



Fluid dynamics is the low energy effective description of any interacting quan- 
tum field theory, in the regime in which fluctuations have sufficiently long 
wavelengths. Hydrodynamics provides a statistical description on macroscopic 
scales of the collective physics of a large number of microscopic constituents. 
A classical reference on the subject is [31], while for a more specific discussion 
of relativistic fluids one can also refer to 

Quantum systems in equilibrium can be described by the grand canonical 
ensemble, given the temperature and chemical potentials for various conserved 
charges. Observables of the system are given by correlation functions in the 
grand canonical density matrix. When the system is perturbed away from 
equilibrium, for fluctuations whose wavelengths are large compared to the scale 
set by the local energy density or temperature, the system can be described 
at a macroscopic level in terms of fluid dynamics. The variables of such a 
description are the local densities of all conserved charges together with the 
local fluid velocities 

The intuition is that sufficiently long- wavelength fluctuations correspond to 
variations that are slow on the scale of the local energy density/temperature. 
Therefore, the system can be considered at equilibrium patchwise and locally 
the temperature can be seen as constant. Therefore while the grand canon- 
ical ensemble remains a valid approach to describe the physics locally, fluid 
dynamics describes macroscopically how different local domains roughly at 
equilibrium interact and exchange thermodynamic quantities. 

Fluid dynamics is a quite surprising field already as a description of classical 
fluids. Although classical fluid dynamics equations have been intensively stud- 
ied for almost two centuries, their extremely rich phenomenology remains not 
completely understood. Well known open problems are the so called "Navier- 
Stokes existence and smoothness problems" |30] , that is the demonstration, for 
non-relativistic incompressible viscous fluids described by the Navier-Stokes 
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equations, of the existence of globally regular solutions. Moreover interest- 
ing phenomena such as turbulence remain to be completely clarified. The 
holographic mapping of the fiuid dynamical system into classical gravitational 
dynamics, that will be the topic of the next chapter, could be a useful tool to 
deal with open problems in both fields. 

In this chapter we will start reviewing relativistic fiuid dynamics [32| [33], 
then we will summarise some aspects of conformally invariant fiuids [31]. Fi- 
nally, we will discuss the non- relativistic limit of fiuid dynamics |3H [37] . 

2.1 Relativistic ideal fluids 

Let us start by defining in a more formal way the hydrodynamic regime. In 
any interacting system there is an intrinsic length scale, the "mean free path 
length" fjnfp, which constitutes the characteristic length scale of the system. In 
the kinetic theory context, £mfp is the mean free motion of the constituents be- 
tween successive interactions. The hydrodynamic regime is therefore archived 
when considering description of the system on at length scales which are large 
compared to imfp- At this scale, microscopic inhomogeneities are sufiiciently 
smeared out. 

The equations of relativistic fiuid dynamics can be summarised as 

V^T'^" = , V^J^ = (2.1) 

where the first states the conservation of the stress tensor T^'^ and the second 
one refers to the conservation of charge currents Jj, where I = {1,2,---} 
indexes the set of conserved charges characterizing the system. To characterize 
a system it is necessary to specify the stress tensor and charge currents in terms 
of the basic thermodynamic variables. 

Let us consider a QFT living on a d-spacetime dimensional background 
Bd with coordinates and non-dynamical metric g^j,^. The thermodynamic 
variables of the system are the fiuid velocity u^, normalized as 

u^g''''u, = -l; (2.2) 

the local energy density p and charge densities qi, that are seen as extrinsic 
quantities; pressure p, temperature T, and chemical potentials fij that are 
intrinsic quantities determined by the equation of state. 

Let us focus first on the volume element of an ideal fiuid which has no 
dissipation, as seen in the reference frame where it is at rest. In this frame 
Pascal's law is expected to hold. The pressure exerted by a given element 
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of fluid is the same in aU directions and is perpendicular to the surface on 
which it acts. The i-th component of the force acting on a surface element 
di can be written as T^^dfj. Therefore in the local rest frame we can write 
T'^^dfj = pdfj, that is Tij = pSij. The component is the energy density p, 
while components T''^'*, that refer to the momentum density, are zero in the 
local rest frame. 

Summarising, in the frame in which the ideal fluid is at rest, the energy- 
momentum tensor has the diagonal form 



p 














p 














p 














p 



(2.3) 



In a generic reference frame the energy-momentum tensor that we wrote in 
the rest frame will be 



(2.4) 



where is the fluid velocity. Similarly, in the local rest frame the components 
of the charge current are the charge density itself .The particle flux will have a 
time component given by the density of particles, while the spatial components 
will form the spatial flux vector. Therefore we can write the charge current as 



ideal 



(2.5) 



Another relevant current to be defined is the entropy current, that similarly 
to the charge current can be written as 



( -^5 ) ideal 



SU<' 



(2.6) 



where s is the local scalar entropy density. The entropy current describes the 
variations of entropy in the fluid. It is well known that from the conservation 
equations and standard thermodynamic relations one flnds that the entropy 
current is conserved for an ideal fluid 



Vm (^5")ideal = 0. 



(2.7) 



A convenient notation can be formulated by observing that the d- velocity 
is oriented along the temporal direction. Therefore, it is possible to use 
to decompose the spacetime into spatial slices with induced metric 



(2.8) 
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where P'^'^ can be seen as projector onto spatial directions, which has the 
following properties 

P^" u, = 0, P^P Pp, = Pt = P^' 9p., P/ = - 1. (2.9) 

The ideal fluid stress tensor may be re- written in terms of the new projector 

as 

iTn^a,^, = Pu'u'^+pP"'- (2.10) 

2.2 Relativist ic dissipative fluids 

Real fluids manifest dissipative effects, resulting in the creation of entropy 
consistently with the second law of thermodynamics. Indeed, in general the 
conservation of entropy current does not hold. These dissipative effects op- 
erate in the fluid to equilibrate it when perturbed away from a given initial 
equilibrium conflguration. From a microscopic point of view, dissipative effects 
are due to interaction terms in the underlying QFT of which fluid dynamics 
is a macroscopic description. Therefore, we expect the terms incorporating 
dissipative effects in the the stress tensor and charge currents to depend on 
the coupling constants of the underlying quantum fleld theory. 

Following [34J we will introduce dissipative terms with a procedure inspired 
by the usual way in which effective fleld theories are modeled, taking into ac- 
count all possible terms that can appear in the effective Lagrangian consistent 
with the underlying symmetry at the order required. In the same way, we will 
introduce in the hydrodynamic approximation all possible operators - con- 
structed as derivatives of the velocity fleld and thermodynamic variables -, 
consistent with the symmetries. This procedure is quite natural in that the 
hydrodynamics can be seen as an effective theory for an underlying microscopic 
quantum fleld theory. 

To clarify the approach we will follow, in preparation for the construction 
to flrst order in the gradient expansion of the dissipative terms of the stress 
tensor, we will look for all possible symmetric two tensors built out of the 
gradients of the velocity fleld and thermodynamic variables. Clearly the con- 
servation equations (12. ip - which are flrst order in derivatives - can be used to 
simplify the expression for the flrst order stress tensor. In this way correction 
terms can be expressed as derivatives of the velocity fleld alone. This process 
can clearly be iterated to higher orders. 
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Let us start introducing generic dissipative terms to the stress tensor, 11^'", 
and charge currents, T'^, that we will have to determine 

(^^Hdissipative = pU^W^+pig'^+U^Un+Il'^ 
(^n^ssipative = 1lU^ + r^. (2.11) 

To proceed in a rather systematic way to the enumeration of dissipative 
terms, it will be necessary - as we have done in the case of the ideal fluid - to 
define a reference frame in which new terms will be constructed. The choice 
of the frame is clearly related to the choice of the fluid velocity. 

In the discussion of ideal fluids we defined a velocity field such that, in the 
local rest frame of a fluid element, the stress tensor components longitudinal to 
the velocity gave the local energy density in the fluid. One can make the same 
requirement in the study of dissipative fluids. The corresponding gauge choice 
is known as Landau frame and is defined by demanding that the dissipative 
contributions be orthogonal to the velocity field, i.e. 

U^''u^ = 0, T^u^ = 0. (2.12) 

The idea in adopting the Landau frame is to define the velocity field u'^ to 
be given by the unique (normalised) timelike eigenvector of T^^,, so that the 
definition of the velocity field is tied to the energy-momentum transport in the 
system. In what follows, we will work in the Landau frame and we will use 
eq. f l2.12p to constrain the dissipative terms of the stress tensor and charge 
currents. 

To warm up, we may start by considering the decomposition of the velocity 
gradient V^u^ into a part along the velocity - given by the acceleration -, 
and a transverse part. The transverse part can be decomposed into a trace, the 
divergence 9, and a traceless part. Symmetric components are to be identified 
as the shear a^^ of the fluid, while the antisymmetric components can be seen 
as the vorticity uj^" . Therefore we can write the velocity gradient as 



d — 1 



(2.13) 
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The divergence, acceleration, shear, and vorticity, have been defined^ as 

e = ^^u^ = P^^-^^u, (2.14) 
a^" = u" y = ^u^" (2.15) 

V^^u""^ + u^^ a"^ - e P^^ = P^" P"^ S/i^aUp) - 9 P^"" (2.16) 

^^u^, ^ y[^,^u] ^ ^[^, ^u] ^ p^a pup y^^^^^ ^ ^2.17) 

and satisfy by construction the following properties: 

a^u^ = f) P^,a^ = a,, (2.18) 

a^''u, = Q, a^pPp, = a^^, < = 0, (2.19) 
uj^''u^ = Q, uj^'Pp, = u^,, < = 0. (2.20) 

Note that above and in what follow we adopt standard symmetrization and 
anti-symmetrization conventions. For any tensor Fab we define the symmetric 
part P(afe) = I [Fab + Fba) and the anti-symmetric part F[ab] = | {Fab - Fba) 
respectively. Moreover we indicate with ^ the velocity projected covariant 
derivative: ^ = u^V ^. 

As mentioned above, the conservation equation (12.1 p can be used to simplify 
the expression of the first order stress tensor. Useful relations can be derived 
projecting the conservation equation, along the velocity field and transversally 

«^V^(TnMeai = =^ (p + p)VX + «'^V^P = (2.21) 
Pua (T^^),,,,i = =^ P/V,P + (p + p) P.„ = • 

Given the identities f l2.2ip that we derived, it is possible in the definition of 
dissipative terms for the stress-tensor to take into account for only symmet- 
ric two tensors built from the velocity gradients. Using the Landau frame 
condition, one is able to single out two candidate dissipative terms: 

nf^) = -2 7]a^"' -Ce P^"" , (2.22) 

that appear associated with two new parameters the shear viscosity, t], and 
the bulk viscosity, (. 



^In a more compact form we can write the projectors as 

pM" pup y^^ ^ PPf^VpU'') , 
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In the same way, it is possible spoiling the conservation of charge current 
to simplify the search for new terms. In particular it is possible to use the 
conservation equations to eliminate the acceleration term, limiting the new 
contribution to ones which are first order in the derivatives of the thermo- 
dynamic variables p and qj and also the velocity field. Another potential 
contribution arises from the pseudo- vector 

^^ = e„;3/«"VV. (2.23) 

This possible term will be responsible for mixing contributions with different 
parity structure at first order^. 

The new terms, at first order, for the charge current are 

'^(1)1 = -^iJ ^^<1J - 7/ P^" ^uP -Ud^ , (2.24) 

where x/j is the matrix of charge diffusion coefficients, 7/ indicates the con- 
tribution of the energy density to the charge current, and 0/ which are the 
pseudo-vector transport coefficients, in terms of the temperature T and of the 
chemical potential the current reads 

^a)/ = -^'iJ P""" V. (^) - r - 7/ P^" V.T . (2.25) 



Summing up, the stress tensor and the current of charge for a dissipative 
fluid at leading order in gradient expansion are 

= g,M'^-x,jP^^V.(^)-0/r-7/P^"V.T. (2.26) 

where the set of transport coefficients {rj, (, xjj, 77, Uj} are to be determined 
to completely specify the relativistic viscous fluid. In principle transport co- 
efficients can be worked out from the fundamental QFT underlying the fluid 
dynamic description. 

The dissipative fluid's entropy is not conserved, therefore the second law 
implies that the entropy current should have non-negative divergence, i.e. 

V^4' > . (2.27) 

^This kind of contribution for charge current appears in fluid dynamics derived from the 
AdS/CFT correspondence for the A/" = 4 Super Yang-Mills fluid and is linked to Chern- 
Simons couplings of the bulk gauge field in the gravitational description [351 ES] • 
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In principle, the discussion can be extended to determine higher order 
terms in derivative expansion. However, usually these terms are irrelevant in 
the hydrodynamic regime. Indeed, these terms are less and less important 
moving towards longer wavelengths and are therefore suppressed in low energy 
description. On the contrary terms at first order are important since their 
presence defines a channel for the fiuid to relax back to equilibrium after a 
perturbation. 



2.3 Conformal hydrodynamics 

In this section we will introduce conformal fiuids, that can be thought of as 
the low energy macroscopic description of field theories which are conformally 
invariant. We will start by discussing the restriction on first order terms in the 
derivative expansion of the stress-tensor for a generic dissipative fiuid arising 
from the requirement of conformal invariance. In the next section, we will 
review the useful Weyl covariant formalism [44J, before singling out operators 
at second order that can appear in conformal hydrodynamics. This discussion 
is motived by the construction of gravitational duals to hydrodynamics that 
will be the topic of the next chapter. 

Having in mind a relativistic fiuid on a background manifold Bd with metric 
Ofiu, we can consider a Weyl transformation of the metric 

9,u = e'^g,,, g^'^ = e-'^r'. (2.28) 

Remembering that the velocity field is normalised as m'^ = — 1, it will scale 
under a Weyl transformation fl2.28l) as 

u>' = e-^u^. (2.29) 

Combining these facts one finds thate the spatial projector also transforms 
homogeneously 

pMi^ = gi^'^ + ut'u" = e-^^ P'"' . (2.30) 

A generic tensor QJ;^^^'''^'^ is said to be conformally invariant if it transforms 
homogeneously under Weyl rescalings of the metric (12.281) , i.e., 

Q = e-''"^Q, (2.31) 

where w is called conformal weight of the tensor and depends on the index 
positions. To have a proper conformally invariant operator in the theory, it 
is also necessary that the dynamical equations for Q remain invariant under 
conformal transformations. 
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Now we can address the question about additional constraints on the stress 
tensor, at first order, due to conformal invariance. From the conservation 
equations for the energy- momentum tensor (12.11) . one finds that it transforms 
homogeneously under Weyl rescalings of the metric with weight d + 2 

rptlU ^ ^-{d+2)^f,,U ^ (2.32) 

provided that its trace is zero T^^ = 0. 

The null trace condition implies for an ideal fluid that 

T; = =^ p=-L_p_ (2.33) 

This relation between pressure and energy density fixes the speed of sound in 
conformal fluids as a function of the spacetime dimension: 

1 



The charge current transforms homogeneously with weight d under confor- 
mal transformations, i.e. 

J^' = e-^'^>. (2.34) 

The last pieces to recollect are the scaling dimensions of the thermody- 
namic variables. It is possible to see that the temperature scales under a Weyl 
transformation with weight 1 

T = e-^f, (2.35) 

and this implies from the thermodynamic Gibbs-Duhem relation, p + p = 
sT + hj qi, that the chemical potentials /i/ also have weight 1. Moreover 
from eq. f l2.32p it follows that the energy density transforms under a Weyl 
transformation as 

p = e~'^'^p. (2.36) 

Using the properties of transformation listed above, one is lead to the 
following expression for the stress tensor of an ideal fluid 

{Tn,^^^, = c^T'i9"' + du^un , (2.37) 

where a is a dimensionless normalization constant fixed by the underlying mi- 
croscopic conformal field theory. 



The next step is to discuss dissipative corrections at first order to the ideal 
fluid stress tensor. The most convenient approach, is to single out all of the 
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operators at order one in the derivative expansion with the right symmetries 
transforming homogeneously under conformal transformations. With a bit of 
work, it is possible to show that the covariant derivative of transforms 
inhomogeneously, i.e. 



(2.38) 



where it is useful as an intermediate result that the Christoffel symbols trans- 
form as 

= + d,<p + 5;: - gx, 5.0 ■ (2.39) 

From eq. fl2.38p it is possible to obtain the transformation laws of fluid 
dynamical quantities 



9 



a 



A 



e-* (0 + (rf-l) 



u^^X = e 



1 



Urye 



d- 1 



(2.40) 



In the last equation one uses the fact that all epsilon symbols should be 
treated as tensor densities in curved space^. 



From eq. f l2.4Up we can deduce the restrictions that one must impose 
on a conformal dissipative fluid at flrst order. Since 9 (and a'^) transform 
inhomogeneously under Weyl transformations, the bulk viscosity should vanish 
for a conformal fluid C = 0- With the same procedure, for the charge current, 
it turns out that the contribution from the chemical potential and temperature 
should appear in the combination yU//T. Moreover since the gradient of the 
temperature P^" V JT transforms inhomogeneously, the coeflicient of the term 
in which it appears will be zero 77 = 0. 

Summing up, one flnds that the stress tensor and the charge current, at 
flrst order, for a viscous fluid are 

= qiU^-^ijP'"'V,[^)-Uji'' , (2.41) 

where we have used the generalized Stefan-Boltzmann expression for the pres- 
sure. 

■^It is possible to show [3^ that the epsilon symbols scale as metric determinants i.e., 
^ap-yS oc y/—g, and e"'^T'5 q<- --^=. Therefore the correct scaling behavior under conformal 
transformations will be ea^js — e'*'^ €afijS and e^^f'J — ~ai3fS^ 
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2.4 Weyl covariant formulation of conformal 
hydrodynamics 

The Weyl covariant formulation of conformal hydrodynamics is a quite pow- 
erful tool to explore second order corrections to the stress tensor [33]. 

On the conformal class of metrics C, defined on the background manifold 
Bd, it is possible to define a derivative operator that transforms covariantly 
under Weyl transformations. Let us start by defining a torsionless connection 
yWeyi^ called the Weyl connection, such that for every metric in the conformal 
class C there exists a connection one-form for which 

VT''9,. = 2A^g,, . (2.42) 

It is possible to define a Weyl covariant derivative as 

^^^^Weyl^^^^ (2.43) 

or in a more explicit form 

+ (gxa - s'x ^« - ^a) q:::. + ■■■ 

- (gxu - 5^ A - S: Ax) Q^::. . (2.44) 

The Weyl covariant derivative transforms in a covariant way under con- 
formal transformations, when acting on a conformally invariant tensor. In 
particular 

VxQt::.=e-'"^VxQ>;:::. , (2.45) 

in other words, the Weyl covariant derivative of a conformally invariant tensor 
transforms homogeneously with the same weight as the tensor itself. 

From the definition, it follows that the Weyl connection is metric compat- 
ible 

Vo^g^, = , (2.46) 

since w = —2 for g^^. 

Requiring that the new derivative of the fluid velocity is transverse and 
traceless 

m"P„m^ = 0, P„m" = 0. (2.47) 
the connection one-form ^ is uniquely determined as 

A = Vau^ - V\ix = a^- 9 . (2.48) 
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Once the Weyl covariant derivative has been defined, one is lead to define 
an associated curvature as the commutator between two covariant derivatives. 
In this case, the standard procedure, can be followed keeping in mind some 
subtleties that we are going to describe. For a given covariant vector field 

= e~'^^V^ , it is possible to define a Weyl- invariant Riemann curvature 
tensor TZ/j^i^xa and a Weyl-invariant field strength J^^,^ as 

[D^,V,]Vx^wJ'^,Vx-n^.x''Va , (2.49) 



where 



and 



(2.50) 



7^ 



+ {gx.Af' - 5iA. - 5Ma] - - 5^^/3] 

gx^A^ - 5iA, - 5^ a] [9p.A^ - ^^A. - S'^Ap] . (2.51) 

The Weyl-invariant Riemann curvature tensor can be rewritten in the more 
compact form 

T^UvXa = Ru.v\a + 4 (^[^S'l/] [A(^f] aAp + AaAp - ^S'a/S^ " ^nv g\a , (2.52) 

where 7-"^^ is the field strength for the Weyl connection defined above. 

Both the Riemann tensor TZf^^xa and the field strength J^^^ are invariant 
under Weyl tranformations 



n 



T = T 



111/ • 



(2.53) 



There are relevant differences between the Weyl covariant Riemann tensor and 
the ordinary one in the symmetries they have: 



T^fMi/Xcr H" T^fMi/aX 2 J~ [jj/ gXcT 1 

T^UvXa — T^Xaixu — ^]^9v\[X^^^ — ^^iv gXa + ^Xa Q^v , 



(2.54) 



In analogy with the ordinary case, given the Riemann tensor one can define 
a Ricci tensor 



TZ = 7^ " 



(2.55) 



= R^^ -{d-2) (V^ A + A^,Au - A^g^iu) - g^,uS/xA^ - 
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and a Ricci scalar 

7^ = 7^<,° = R-2{d- l)VxA^ + {d-2){d- 1)A^ , (2.56) 
that transform as 

n^u = 'ji^.u, n = e-^'t"JZ . (2.57) 

A general class of tensor that appears in the study of conformal manifolds, 
are the the Weyl-covariant tensors that are independent of the background fluid 
velocity. A relevant example of this class of operators is the Weyl curvature 
C^uXa, that can be defined as the trace free part of the Riemann tensor. For 
d > 3 it is 

CfiuXcT = Tl^uXfj + 4:5"^gy][\5^-^Sal3 = C^i^Xa — avQXu = ^^'^ CfiuXa (2.58) 

where we introduced S^^, the Schouten tensor, defined as 

d-2{ 2id-l] 
= S^, - (v,A, + A,A, - ^g^)j - = S,, . (2.59) 

The Weyl Tensor C^uXu = C^uXa + J^^iudXa defined in eq. fl2.58p is a conformal 
tensor. Moreover, it turns out that C^^Xa has the same symmetry properties 
as the Riemann Tensor Rf^^Xa, i-G. 

C^uXa = ~Cu^j.Xa = ~C^waX = Cxa^v (2.60) 
C^aA" = . (2.61) 

From equation f l2.6ip . it follows that C^aupu"^ is a symmetric traceless and 
transverse tensor. 

Given the Schouten and the Weyl curvature tensors, one can recast other 
curvature tensors in terms of these new tensors: 

T^fiuXa = C^uXa — ^[^gu\[X^^„\^aP-, 

7^ = 2(rf-l)5A^ (2.62) 
'Tip.v = {d — 2)S^„ + Sx^g^iiy, 



2.5 Conformal dissipative fluids 

Now we can use the Weyl covariant formalism to provide a classification, up 
to the second order in the derivative expansion, of the operators that could 



46 



Elements of fluid dynamics 



appear in the stress tensor of a conformal dissipative fluid. Second order term 
win be relevant in the discnssion of holographic duals of conformal fluids. 

Let us start by writing the conservation equations for the stress tensor in 
terms of the new covariant derivative, i.e. 

+ {9^aA''-5;Aa-6'^A,) T^^ 

= V^^^^ (2.63) 

where inhomogeneous terms cancel out, remembering that T^^ = and w = 
d + 2. In the presence of a conformal anomaly, that can be present for CFTs 
on curved manifolds in even spacetime dimensions, the dynamical equation for 
Tfxi, spells out 

P^yM^ = V^T''^ + X {T^"^ -W) ^0 , (2.64) 

where W is the trace anomaly. 

Given the deflnition of the Weyl covariant derivative, it is possible also to 
rewrite the building blocks of the stress tensor at flrst order in the covariant 
formulation 

a^"' = VK"^ , uj^''' = -VK"^ , (2.65) 

which have weight w = 3. 

To introduce second derivative terms, we need to pin down all of the pos- 
sible two derivative operators that transform homogeneously under conformal 
transformations. The flst kind of term we can have are those that involve the 
squares of the flrst derivative operators, i.e. 

a'; a"" = e"^ , cu'^ cu"" = e"" cD'^, (ZS"" , < cu"" = e'^^ a^"^ ^""^ 

(2.66) 

Moreover there exist also two derivative terms built from the fluid velocity. 
The possible terms of this kind are 

V^V.u^ = V^a,^ + V^u,^ = e-'^V^V^u^ , (2.67) 
Vxa^, = e'f'Vxa^u, (2.68) 
Vxco^^^e^Vx^f,. . (2.69) 



An other kind of second order terms is one which involves the temperature 
T and various chemical potentials /ij. After some effort, it turns out that all 
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of the possible operators in this group are 



(2.70) 



In addition to the terms hsted above, there is a one last term to consider, 
obtained by contracting the curvature tensors in the Weyl tensor fl2.58p with 
two velocity vectors 

C^^.^u''u^ = C,^,pu''uP (2.71) 
that is a symmetric, traceless and transverse tensor. 

Having classified all possible operators that can enter in the construction 
of dissipative terms for a viscous conformal fiuid, it is possible to proceed in 
isolating the terms that will be useful later on. In particular, for sake of sim- 
plicity, we will construct a dissipative generalisation of a conformal fluid with 
no conserved charges. For this choice, there are five possible symmetric trace- 
less tensors which transform homogeneously under Weyl rescalings at second 
order: 

If = 2 V^a^,, = u^, (2.72) 

If = 4 a'l Tf = 2 a"^'^ co'l (2.73) 
= a;"<^w'^>. (2.74) 

We summarise all the symmetric traceless tensors which transform homoge- 
neously under Weyl rescalings up to the second order in table 2.1. 



Having completed the exploration of possible dissipative terms up to the 
second order, we are finally able to write the correction to the stress tensor 
for a conformal viscous fiuid that, for simplicity, we will assume to have no 
conserved charge. The general contributions to the stress tensor are: 

TJ" = aT'^ (^^" + du>'u'') -2ri a"" 
U"^^^ = -2r/a^^ 

Hg; = r^r/Tr + «:Tf + AiTr + A2 If + A3 ir . (2.75) 

Where {77, r^r, k, Aj} for i = {1, 2, 3} is a set of six transport coefficients char- 
acterizing the viscous fiuid. 
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First order: a^'^ 



Second order: 


If 


= 2 VaO-f^u 












= 4 




IT 


= 2 




■^5 





Table 2.1: Symmetric and homogeneous traceless operator for a conformal 
viscous non-charged fluid up to the second order. 

2.6 Nonrelativistic limit of fluid dynamics 

In this section we will discuss, following the non-relativistic scaling limit 
of the fluid dynamical equations of a non-charged fluid (we will not assume, 
for the moment, conformal invariance), i.e 

V^T^'^ = , (2.76) 

with 

T^^ = pu^'u'' + pP^"" - 27]cx'"' - COP^'' + ... (2.77) 

where p is the pressure^, p the energy density, t] the shear viscosity, ( the 
bulk viscosity of the fluid, to the incompressible non-relativistic Navier-Stokes 
equations 

dv — * — * 

— + Vt7= -Vp + z/V^{T+/ (2.78) 

ot 

V-v = (2.79) 

where v is the fluid velocity, p the fluid pressure, z/ the shear viscosity and / 
an external forcing function. 

In particular, we will discuss the non-relativistic limit 





1 


6x 




Te 




1 


5t 






v' 


~ e 


5v 




e 


^ 



(2.80) 



'^In this section p will indicate the relativistic pressure and p the nonrelativistic one. 
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where 6x is a spatial length scale, St the temporal scale while and 5p have 
to be seen as estimates of the magnitude of velocity and pressure fluctuations 
on a configuration of fluid at rest and in equilibrium. 

The meaning of this limit can be understood by remembering that it is 
necessary to scale to long distances to be in the fluid dynamical regime. Time 
intervals should scale like spatial intervals squared, as it is easy to realise 
looking at the dispersion relation for shear waves, u = iuk"^ . These two scalings 
determine the scaling law for velocities. Finally the pressure fluctuations are 
to be scaled such that they cannot accelerate the fluid velocities outside this 
scaling limit. 

Let us start by considering the flow of a fluid on a spacetime with a metric 
that can be cast in the form 

G^u = gfiu + H^i, , (2.81) 

where H^i, is an arbitrary small fluctuation of a background metric g^^u that 
we want to write as^ 

g^^dx'^dx" = -dt^ + gijdx'dx^ . (2.82) 

The fluid flow on the background space equipped with the metric can 
be reinterpreted as a forced flow on the space with metric g^j_v [4J, mapping 
the velocity field on the space G^i^ to a velocity field on g^^. In con- 
structing the map the normalisation of the velocity field = = —1 must 
be maintained. 

The fluid velocity on the space G^y can be written as 



u^ = ^^[l,V] , (2.83) 



where V is the d—1 spatial vector with components V"^. is a d component 
object with components (1,^) and V'^ is G^uV^V. 

A perturbative expansion of the velocity at the first order in the small 
fluctuation H^^ gives 



1 



^l-g,,V^V^ 



1,V) (2.84) 



^This form of the metric g^i, can be seen as a simply a choice of coordinate system, for 
a large class of metrics. 
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where we have expressed the series in terms of u^^ the new velocity of the fluid 
referring to the metric g^^p. The idea that we want to follow is to consider 
terms in which H^p appears as contributions to an eflective forcing function 
in the non-relativistic Navier-Stokes equation (12 .78 p . 

We want to consider H^p as a small amplitude and long distance fluctuation 
on a uniform fluid at rest - it is always a solution to the equations of fluid 
dynamics on the space with metric g^y. The fluid at rest will be described as 
a fluid with pressure po, energy density po on a manifold that is "close" to (^^j,. 
We set the components of the fluctuation to 

^00 = e'^hoo{ex\eh) 
Hoi = eAi{ex\eH) 
Hij = e^hij{ex\eH) 

and we assume consistently that the velocity and pressure satisfy 

V' = ev'{ex\eH) 

Po + Vo 

where e is a parameter arbitrarily small and one can deflne a covariant vector 
v'^ = The normalisation of the pressure fluctuations, po + -Po, has been 

introduced for future convenience. Energy density p and viscosity u also scale 
in a similar fashion. 

The conservation equations are expected to reduce to the classical Navier- 
Stokes equations in this limit. Let us start by considering the temporal com- 
ponent of V^T'^'^, we get 

V,r-° = [p. + 0{e^) 

Pe = Po + Po 

therefore, for e — )■ 0, this equation reduces to 

Viv' = 0, (2.87) 

where Vj is the covariant derivative with respect to the purely spatial metric 
gij. Spatial components of the conservation equations, after some manipula- 
tion, can be cast in the form 

-2^V, r y ^^.^ _cV.V-t;-r +0(e^) (2.88) 
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where the forcing function / is defined as 

f = pe \^-^ doA, + v^d Ajj . (2.89) 

The coefficient of can be rewritten in a more transparent form. We begin 
using the expression 

to get 

V> + dov' + v.Vv' - V (V^t;* + Ry) = - doA' + Fy (2.90) 

where we defined a field strength Fij = diAj — djAi associated to the vector 
field Ai, and a kinematical viscosity u = rj/pe of the fiuid. It this then possible 
to redefine variables, splitting the gauge field A^ into its pure curl and pure 
divergence parts, i.e. 

Ai = ai + ViX , (2.91) 

where V^a* = 0, and one finds 

fij = dittj - djai = Fij. (2.92) 
It is also useful to define an effective pressure 

Pe=p - ^hoo + X. (2.93) 

In terms of the new variables, equation fl2.88p becomes 

V,pe + doVi + v.Vvi - V {y^Vi + RijV^) = -doa' - fji , (2.94) 

which is indeed the Navier- Stokes equation with a forcing function generated 
by an effective background electromagnetic field on the effectively charged 
fluid. 



In the case of a charged fluid, the discussion follows the same track, adding 
the set of conservation equations for conserved charges (12.11) . It turns out that, 
in the non-relativistic limit, these equations lead to supplementary conditions 
under which the fluid velocity is divergence- free. 

A point that it is useful to assess regards how Cauchy conditions for the 
relativistic fluid are transmitted to the non-relativistic fluid. The Cauchy data 
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of the relativistic dynamical equations consist of d real functions of space: the 
value of the pressure field and the value of the d—1 independent velocity fields 
defined on an initial timeslice ^. In the case of the non-relativistic incompress- 
ible Navier-Stokes equations, one can observe that the derivative of eq. fl2.94p 
results in the condition 

= -V.v^Vjv' - vVR,, + V, [{-uR] + /j) , (2.95) 

that determines the pressure of the fiuid as a function of the fiuid velocity. 
Therefore, the set of independent data necessary to determine the fiow of the 
fiuid is given by o? — 2 real functions that parameterize an arbitrary divergence 
free velocity field. Indeed, two of the degrees of freedom of the equations of 
relativistic fiuid dynamics are projected away in the scaling limit. One can see 
that these two degrees of freedom correspond the fiuctuations of the pressure 
and the divergence of the velocity^. 

2.7 Nonrelativistic symmetries for conformal 
fluids 

As we have discussed in the previous section, the Navier-Stokes equations may 
be obtained as the scaling limit of any relativistic equations of fiuid dynamics. 
For the purpose of our work we are interested in particular in conformal fiuids 
due to their connection with gravity. Therefore it is interesting to assess in 
which form the conformal symmetry descends to a symmetry of the Navier- 
Stokes equations, in the non-relativistic scaling limit we proposed above. 

A dilatation consists of a diffeomorphism 
associated with a Weyl transformation 

= X{uY , f = XT. 



= u^^ f = XT. (2.96) 

^Since the fluid conservation equations are of first order in time, tlie Cauchy data of the 
problem does not include the time derivatives of all these fields. 

As observed in [J] , at the linearized level these two degrees of freedom combine together 
in sound mode fluctuations. Hence the physical fact behind the reduction of Cauchy data 
in the non-relativistic limit is that sound waves are projected out in the scaling limit 



The compound action is therefore: 
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This action is a symmetry of the equations of conformal relativistic fluid dy- 
namics, but it is not a symmetry of the Navier-Stokes equations. Indeed, 
kinematical viscosity u in the Navier-Stokes equations, is proportional to ^, 
hence it is not invariant under dilatation transformations. However as was 
shown in it is possible to modify dilatation transformations into a true 
symmetry of the Navier-Stokes equations by "absorbing" the transformation 
of u into the "anomalous" transformations of time and velocity. 

Let us start by considering special conformal transformations. The scaling 
law for the velocity and temperature fields, under a special conformal transfor- 
mation, may be obtained by combining a diffeomorphism with the appropriate 
Weyl transformation. For an infinitesimal conformal transformation one gets 

Sx'' = -2c.xx^ + x^d" 

= -2 [x'^c^ - x^c'^] - Sx^d^u^" (2.97) 
6T = 2c.xT - dx^d^T. 

To verify the covariance of local equations under these symmetry transforma- 
tions, one has to omit the terms proportional to Sx'^df^, using the following 
expression for the derivative's transformation 

6 (d/s) = 2 [cjjx.d — Xjjc.d + x.cdjj] . 

Under transformation (12.971) . the conformal stress tensor becomes 

5T^^ = 2d{c.x)T^"' + 2{x^c^ -x^c^)T^ + 2{x^c'' -x''c^)T^ -Sx^dxT^" . (2.98) 

Therefore, given the derivative's transformation, any identically traceless stress 
tensor is invariant under special conformal transformations. 

Special conformal transformations (I2.97P give rise to an additive shift in 
the temperature fluctuation, 6T, proportional to x ■ cTq. In order that this 
shift respects the scaling of 6T one has to scale Cq oc and q oc e^. In this 
way the transformations fl2.97p read 

6t = 
Sx' = -I'd 

Sv' = -2cH + t^Cjdjv' (2.99) 
5T = 2{-cH + dx')T + t^CjdjT. 

(2.100) 
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To verify the covariance of the conservation equations under these transfor- 
mations, one has to omit the terms proportional to t^Cjdj and to supplement 
their action with the following derivative transformation 

5{dt)^2tddi 5(di)^0. 

The symmetry generated by c° acts trivially. Indeed it does not act on 
coordinates or velocities, but merely generates a shift, linear in time, of the 
pressure. On the contrary, the symmetries generated by d act nontrivially. 
Under this transformation one gets 

(^Pe = 375" = 2c • x , 
dPo 

and therefore 

5diPe = 2c' . 

Finally S{v^ + v- Vv^) — —2d. Hence, the viscous term - which was responsible 
for the change in Navier-Stokcs equations - is unchanged under the redefined 
transformation. The Navier-Stokes equations are invariant under a conformal 
symmetry group. 

We can conclude this section, for sake of completeness, by giving few more 
details on the conformal symmetry group. 

The group generators of this group are the dilatation D, special conformal 
symmetries Ki, Galilean boosts Bi, the generator of time translations (energy) 
H, momenta Pi and spatial rotations Mjj. The action of these generators on 
velocity fields is given by 



Dv^ = {-2tdt - x™(9„ - 1) V 
KiV^ = -2t5ij + t'^diV^ 
BiV^ = 6ij — tdiV-' 



(2.101) 



Hv^ = —dtv^ 

MikV^ = 5ijv'' - 5kjv' - {x'^di - x'dk)v^ 

While the following commutation relations among generators hold: 

[D, K,] = ~3K, [D, B,] = -B, [D, H] = 2H 

[D, P,] = P, [D, M^j] = [M,j , Pk] = -SikPj + 5,kPi 

[Mij, Kk] = -SikKj + 6jkKi [Mij, Bk] = -SikBj + SjkBi (2.102) 

[M,^,H]=0 [K^,P,] = [Ki,B^]^0 

[K,,H]^-2Bi [H,Bj] = -Pi 
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In addition to the relevant symmetries we enumerated above, the Navier- 
Stokes equations have, as noted before, an infinite dimensional group of trivial 
symmetries that act as a shift of the pressure by an arbitrary function of time. 
These are indeed symmetries since only gradients of the pressure enter the 
Navier-Stokes equations and trivial because the pressure is not really an inde- 
pendent variable, as we discussed above. The action of symmetry generators 
on the pressure is 

Dp, = {-2tdt-x^d.^-2)p, 
Kip, =2x' + ed.pe 

BiPe = - tdiPe 
Hpe = - dtPe 
PiPe = - diPe 
MikPe = - {X^di - x'dk)Pe- 

These generators acting on the pressure do not yield the commutation 
relations, indeed they have additional terms. However spatial derivatives of 
the pressure field correctly transform according to the algebra fl2.102p . For 
this reason, the symmetry algebra (12.1021) is not represented on the pressure 
field itself, but only on its spatial derivatives. 



(2.103) 
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Chapter 3 

Fluid dynamics from gravity 



In this chapter we give a short review of the construction of gravity solutions 
dual to arbitrary fluid flows proposed in [1]. A detailed and comprehensive 
exposition of the subject can be found in [31]. As pointed out in the first 
chapter, the application of hydrodynamic concepts to black holes is part of 
the core intuition of the membrane paradigm [211 ES], wherein one modeled 
the black hole horizon by a membrane equipped with fluid like properties. 
However, in this approach the connection to fluid dynamics is an analogy, 
and therefore provides a qualitative understanding of the physical behaviour 
of black holes via a simpler fluid model. On the contrary, the fluid-gravity 
correspondence - an example of the AdS/CFT correspondence - provides a 
duality between the hydrodynamic description and the gravitational dynamics. 
Indeed, the precise quantitative connection allows one to map fluid solutions 
of the boundary field theory into a black hole solutions in the bulk geometric 
description with regular event horizons. 

The AdS/CFT correspondence conjectures a deep connection between quan- 
tum field theories and theories of gravity [5H1 EH SD]- The correspondence has 
proved to be a powerful tool for the large N limit in which the gravitational 
theory turns classical, and in a simultaneous strong 't Hooft coupling limit that 
suppresses a' corrections to gravitational dynamics. In this limit the conjec- 
ture asserts the equivalence between the effectively classical large dynamics 
of the local single trace operators p„ = A^~^TrOn of gauge theory and the clas- 
sical two derivative equations of Einstein gravity interacting with other fields. 
The AdS / CFT correspondence provides a one to one map from the local and 
relatively simple bulk equations into unfamiliar and extremely nonlocal equa- 
tions for the boundary trace operators p„(x). The authors of [1] proposed a 
construction which can simplify the extremely complicated equations on the 
boundary. 

Two derivative bulk theories of gravity, containing AdS^+i x Mj as a solu- 
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tion, admit a consistent truncation to tlie Einstein equations with a negative 
cosmological constant. The only fluctuating field in this truncation is the Ein- 
stein frame graviton; all other bulk fields are simply set to their background 
AdSd+i X Mi values. This fact implies the existence of a sector of decoupled 
and universal dynamics of the stress tensor in the corresponding dual field 
theories. Indeed, all single trace operators other than the stress tensor may 
consistently be set to zero as the stress tensor undergoes its dynamics. The 
dynamics of the stress tensor happens to be universal since it is governed by 
the same equations of motion in each member of this infinite class of strongly 
coupled CFTs. 

A further simplification in the dynamics of the dual theory can be obtained 
by considering the regime in which the local stress tensor varies on a length 
scale that is large, at any point, compared to a local equilibration length scale 
- intuitively, a "mean free path" - which is set by the "rest frame" energy 
density at the same point. Under this condition, it is expected that boundary 
configurations should be locally thermalized, and therefore well described by 
the equations of boundary fiuid dynamics. In particular, it has been proven 
that the complicated nonlocal T^y dynamics reduces to the familiar boundary 
Navier Stokes equations of fiuid dynamics in this long wavelength limit [1]. 

In |1] a perturbative procedure was proposed, in a boundary derivative 
expansion, to construct a large class of asymptotically AdSs long wavelength 
solutions to Einstein's equations with a negative cosmological constant. These 
solutions are parameterized by a four-velocity field u^{x^) and a temperature 
field T{x^) that have to obey the four dimensional generalized Navier Stokes 
equations V ^jT^^ = where the stress tensor T^^{x^) is a local functional of 
the velocity and the temperature fields. 

The explicit map from the space of solutions of a distinguished set of Navier 
Stokes equations to the space of long wavelength solutions of asymptotically 
AdSrf+i gravity was established up to second order for o? = 4, 5 in ^ and than 
generalized and studied in more detail in [371 El 1121 El SHI [H]. In particu- 
lar, in [12] it was proven that, subject to mild assumptions, these spacetimes 
have regular event horizons. Successively, the construction was extended to 
spacetimes of arbitrary dimensionality [21 115] ■ 

3.1 Black branes in AdS^^+i 

Let us start by analysing the universal sector of a string background of the 
form AdSrf+i x X where X is a compact internal manifold, suitable for obtain- 
ing a consistent string vacuum and whose properties will be irrelevant in the 
following. Einstein's equation of gravity with a negative cosmological constant 
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A follow from the Einstein-Hilbert action 

5buik = ^-j^ I d^+'x (i? - 2 A) . (3.1) 

Choosing RAds = 1; Einstein's equations are given by 

Emn = Rmn — - GmnR — ^ Gmn = (3.2) 

RMN + dGMN = ^. R = -did+l), (3.3) 

where Gm,n is the bulk metric, Emn is the Einstein tensor and M,N = 
l,...,d + l\ 

An obvious solution of these equations is AdS spacetime with unit radius 
ds^ = — + r^r]^,dx^dx''), (3.4) 

where now fi,!/ = 1 . . . d. 

A different solution to Einstein's equations is the planar Schwarzschild- 
AdSd+i black brane which in standard coordinates is written as 

dT"^ 

ds^ = -r^ f(b r) dt^ + , , + 5,, dx' dx^ , (3.5) 

j[br) 

/W = 1 - ^ • (3-6) 

Eq. (13. 5 p can be seen as one-parameter family of solutions labeled by the 
horizon size r+, which sets the temperature of the black hole 

^ = A • (3-7) 

From the black brane solution (13. 5p one can generate a d parameter family 
of boosted black brane solutions by boosting the solution along the translation- 
ally invariant spatial directions x*. In Schwarzschild like coordinates, the new 
solution can be written as 

dv^ 

= -YJT^ + i-fi^K^u + V,,) dx^dx-' (3.8) 



A'kT 



^We adopt the convention of using upper case Latin indices {Af, A^, • • • } to denote bulk 
directions, while lower case Greek indices {/i, i/, • • • } refer to field theory or boundary direc- 
tions. Lower case Latin indices are used, instead, {«, j, • • • } to denote the spatial directions 
in the boundary. Therefore we have M,N — and ^x.v — 1 . . .d. 
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where 

u"" = , (3.9) 
u' = , (3.10) 

and the velocities Pi are all constants with = (]\ and P^'^ = u^m^ + rjf^" 
is the projector onto spatial directions. The temperature T associated to this 
solution is constant. The parameters which characterize the bulk solution 
are now the basic hydrodynamical degrees of freedom, temperature and the 
velocity of the black hole. 

More generally, one can consider solution of the form (13.81) where g^i, is an 
arbitrary constant boundary metric of signature {d—1,1), and u'^ is a generic 
constant unit normalized d velocity, i.e. 

g^uU^u"" = -1, = + u^,Uu, b = • (3.11) 

4:711 

Clearly the new metric is quite redundant, since a d{d + 3)/2 parameter set of 
metrics (13.81) are all coordinate equivalent. Indeed g^i, can be sent in rj^^, by an 
appropriate linear coordinate transformation x'^ — )■ A^^^x^, and can subse- 
quently be set to (1, ... 0) by a boundary Lorentz transformation. Moreover 
b can be set to one by uniform rescaling of boundary coordinates coupled with 
a rescaling of r. 

For our purposes it is useful, to introduce an even more redundant formu- 
lation of solution (13. 8p . obtained with via a transformation f — )■ e~'^r, i.e. 



^ (dr +^^^'^^^^) + f.2 (^-f(lf)u^u^dx^dx'' + P^^dx^dx"^ (3.12) 

where g^j,^, u^, and b defined in the previous equation, while 

^.. = e^^(^^V, «. = e'^(^'%„ 6 = e*(-^)6, (3.13) 

and (f){x^) is an arbitrary function. 

Metrics of the form (I3.12p describe the same bulk geometry but can be 
interpreted as distinct, though Weyl equivalent, boundary description by reg- 
ulating them inequivalently near the boundary. In particular, following the 
references above, spacetimes described by eq. (I3.12p will be regulated on slices 
of constant f and consequently seen as states in a conformal field theory on a 
space with metric g^u{x). 



\2 
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The non-anomalous part of the boundary stress tensor dual to metric f l3.12p 
is given by 

^/^^ = TTTT^ f:^ ^Sf^'' + du^Uu) , (3.14) 

which corresponds to the stress tensor of an ideal conformal fluid with a pres- 
sure p = 1/ {IGnGAdsb'^) and without any vorticity or shear strain rate. 

It is useful to observe that the velocity field is the unique time-like eigen- 
vector of the stress tensor 

T>'^^,A A-^-ii^ (3,15) 

and the inverse temperature field b is simply related to its eigenvalue. 



3.2 Zeroth order ansatz for the metric 

In the previous section, we described locally asymptotically AdS^+i exact so- 
lutions to Einstein's equations f l5.65p which have, as a holographic boundary 
stress tensor, the ideal conformal fluid stress tensor. This is an expected result 
as the solution is stationary and therefore corresponds to the global thermal 
equilibrium. 

Now, to describe hydrodynamics we should perturb the system away from 
global equilibrium. The natural way to do this is to consider thermodynamic 
variables varying along the boundary directions, that is promoting the param- 
eters b {viz. the temperature T), (3i {viz. the velocity m^) to functions of the 
boundary coordinates. In particular, we want discuss solutions to Einstein 
equations with slowly varying boundary stress tensors on a boundary mani- 
fold that has a weakly curved boundary metric. These solutions should be 
locally patch-wise equilibrated, therefore corresponding bulk solution should 
approximately be given by patching together tubes of the uniform black brane 
solutions, extending from local patches on the boundary into the bulk. 

Let us consider a locally asymptotically AdS^+i solution to Einstein's equa- 
tions, whose corresponding boundary stress tensor everywhere has a unique 
time-like eigenvector. Let us promote this eigenvector to a function of 
the boundary coordinates, and let us, in the same way, define the inverse of 
the temperature field as b{x). Now it is useful to define in a precise way in 
which sense we can consider our solution as "slowly varying": denoting by 
Sx{y) the smallest length scale of variation of the stress tensor of the corre- 
sponding solution at the point y, we shell consider 6x{y) ^ b{y). Actually 
b{y) may be interpreted as the effective length scale of equilibration of the 
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field theory at y. Similarly, we say that the boundary metric is weakly curved 
if b{y)'^R{y) ^ 1 (where R{y) is the curvature scalar, or more generally an 
estimate of the largest curvature scale in the problem.) 

An important issue to address is the shape of the radial curves that tubes 
- extending from local patches of uniform black brane on the boundary into 
the bulk - will have. A seeming natural choice is to let tubes follow lines 

= const, in the Schwarzschild coordinates, writing the approximate bulk 
metric as 

= ^"^"^ ^J'^^^^'^^ + (^-f(bf)u^u,dx^dx'' + P^.dx^dx'^) (3.16) 

The problem with this choice is that it has been shown that it leads, in gen- 
eral, to a geometry with non-regular future horizon^. Furthermore, the solution 
f l3.16p has a relevant issue with causality. Let us consider a variation in the 
boundary metric at an arbitrary point y^, inducing an effective force on the 
fluid. The future evolution of the fluid is affected only in the the future bound- 
ary light cone of y'^, that we will denote C{y^Y. The bulk will be affected by 
the variation only in the region of the space-time - B{y^)) - that is the union 
of all the tubes, that originate in the boundary region C{y^). The natural 
requirement for bulk causality is that B{y^) must lie entirely within the fu- 
ture bulk light cone of y^. This is not the case in the geometry described by 
metric f l3.16p . where the tubes run along lines of constant in Schwarzschild 
coordinates. 

The appropriate choice to deal with issues of causality turns out to be to 
promote the parameters of a uniform brane solution, rewritten in Eddington 
Finkelstein coordinates, to slowly varying boundary functions u^{x) and h{x). 
In this way one gets the ansatz for the metric 

ds^ = -2Uf,dx^ (dr + r Audx") + g f,^dx^ dx" + j^u^u^dx^dx" (3.17) 

where as before 

(jf^u = e^'^QiJivix), = e^u^{x), b = e'^6(x), (3.18) 

and g^y{x) is a weakly curved boundary metric. The reason for switching 
to Eddington-Finkelstein coordinates, apart from making issues of regularity 

^Examples of metrics of the form (j3.16p have been discussed in [46l l47l 148) . 

^To be precise this would be strictly true only summing all orders in the fluid expansion. 
Truncation at any finite order could lead to apparent violations of causality over length 
scales of order 1/T. 
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Figure 3.1: (a) Penrose diagram of the uniform black brane where the dashed 
hues denote the future event horizons, (b) The shaded tube extending from 
a local patches of uniform black brane on the boundary into the bulk, and 
running along ingoing null geodesies, indicates the region of spacetime over 
which the solution (13.1 7p is well approximated by a tube of the uniform black 
brane. The two figures are taken from Ref. [31]. 

more transparent, is that for this choice tubes are chosen to run along ingoing 
null geodesies and therefore B{y^) lies entirely within the future bulk light 
cone of the generic point (fig. 13. ip . 

An important observation is that for g^j^^, and b all constants the metrics 
(I3.16P and (I3.17P reduce to equivalent - under a coordinate transformation - 
descriptions of a uniform brane solution. On the contrary, when g^^^, and 
h are functions of the boundaries coordinates x^, the geometry described by 
(13.160 and (I3.17P are inequivalent and can differ qualitatively. Indeed, un- 
der mild assumptions, the metric (I3.17P presents a regular future horizon that 
shields all of the boundary from all future singularities in this space. 

To summarise the general picture, a given boundary patch corresponds to 
an entire tube of width set by the scale of variation in the boundary. In the 
Eddington-Finkelstein coordinates one has to patch together these tubes to 
obtain a regular solution to Einstein's equations and moreover this patching 
can be done order by order in boundary derivatives, just as in fluid dynamics. 
Therefore the metric (13.170 will be the first term in a systematic perturbative 
expansion of a regular solution to Einstein's equations. The perturbative ex- 
pansion parameter is l/h5x. The curvature scale in the metric will be assumed 
to be of the same order as l/5x. 
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3.3 The perturbative expansion in gravity 

The logic of the perturbative expansion of the metric has been described in 
detail in [3] . Let us review the reasoning starting with the ansatz for the bulk 
metric 

Gmn = Cmat + G^M^e + G^^jj^e"^ + O [e^) , (3.19) 

where e is the parameter of the derivative expansion and G^^j^ is the correction 
to the bulk metric at order k that is to be determined with the help of the 
bulk Einstein equation. It will be shown that perturbative solutions to the 
gravitation equation exist only when the velocity and temperature fields obey 
certain equations of motion that will be determined in a perturbative expansion 
order by order in e 

/3, = /3f ) + e /3f ) + O {s') , b = + e b^'^ + O {s') , (3.20) 

where /Sj-™'' and 6*^"^ are all functions of e x^. 

In order to give a precise meaning to the coordinates it is necessary to 
adopt a choice of gauge. In ^ the "background field" gauge was adopted, i. e. 

a, = , Gr^ ocu^ , Tr ((G^o))-^^^")) =0 V n > 0. (3.21) 

A more convenient gauge choice was chosen in [S] 

Grr = 0, = ■ (3.22) 

The meaning of this choice was discussed in |9], where it was observed that, 
under this gauge fixing, lines of constant x'^ are ingoing null geodesies, with 
the radial coordinate r being the affine parameter. 

Having chosen a gauge, e. g. the second one, it is possible to substitute the 
perturbative ansatzs for the metric f l3.19p and for parameters f l3.20p into the 
bulk Einstein equations f l5.65p and expand it order by order in e. It is possible 
to solve the perturbative equation recursively starting form the ansatz for the 
order zero metric defined in the previous section. Let us suppose that we 
have solved the perturbation theory to the {n — 1)*^ order, that is to have 
determined G'-'") for m < n — 1, and have determined the functions /Sj-™'' and 
m < n — 2. The Einstein equation at order n that must be solved can 
be written as 

G(")(x'^) = s„, (3.23) 

where we indicated by EI a linear differential operator of second order in the 
in the variable r alone. An interesting feature of this operator is its ultralocal 
nature in the field theory directions. This depends on the fact that G*^"^ is 
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already of order e", and every boundary derivative appears with an additional 
power of e. Indeed, it turns out that EI is a differential operator only in the 
variable r and is independent of x'^, that it is independent at order n of the 
expansion and depends only on the values of /Sf'^ and 6*^°^ at x^. Hence it is 
possible to solve eq. f l3.23p point by point on the boundary. What changes 
order by order is the source term s„ that can be expressed as a local expression 
of n^^ order in boundary derivatives of and b^^\ as well as of {n — ky^ order 
in /3f \ fe^'^) for all < n - 1. 

As observed in [9j, it is convenient to organize the equations in 

f l3.23p in two classes of equations: equations determining the metric, 

that we call dynamical equations, and a a second set of d equations which are 
essentially constraint equations. 

The constraint equations are those of first order in r derivatives and are 
obtained by contracting the equations with the one-form dr 

4? = EM,ve'^ (3.24) 

where for our considerations C,n = dr. Four of the five constraint equations, 
those with a free index /i (that is with with legs along the boundary direc- 
tion), have a simple interpretation: they are the equations of boundary energy 
momentum conservation: 

^,T^Zi) = , (3.25) 
where T^',^-!) boundary stress tensor dual to the solution expanded up to 

O (e"~^). Recalling that all G^^j^ are local functions of /3j and 6, it is clear that 
'^(n-i) is also a local function (with at most n — 1 derivatives) that respects all 
boundary symmetries. 

The constraint equations can be used to determine and this 

amounts to solving the equations of fluid dynamics at {n — 1)**^ order (I3.25p . 
There is a non-uniqueness in these solutions given by the zero modes ob- 
tained by linearizing the equations of stress energy conservation at zeroth 
order. These may be absorbed into a redefinition of the /3l^\b^^\ and do not 
correspond to a physical non-uniqueness. 

The remaining constraints and E^^ are dynamical equations useful to 
determine the correction that should be added to our initial metric to make it 
a solution of Einstein equations. By exploiting the underlying symmetries of 
the zeroth order solution, specifically the rotational symmetry in the spatial 
sections on the boundary, SO{d — 1), it is possible to decouple the system 
of equations into a set of first order differential operators. Performing this 
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diagonalization of the system of equations one gets a formal solution of the 
form: 

G^^^ = particular(s„) + homogeneous (H). (3.26) 

At this point it is necessary to impose boundary conditions to determine a 
solution uniquely. A possible choice is to require that the solution is normal- 
izable such that the spacetime is asymptotically AdS^+i and also to demand 
regularity at all r 7^ 0. In |3j it was proven that it is possible to require that 
G'-"^ is appropriately normalizable at r = 00 and non-singular at all nonzero 
r, for a generic non-singular and normalizable source Furthermore, if the 
solution at order — 1 is non-singular at all nonzero r, it is guaranteed to 
produce a non-singular source at all nonzero r. This is particularly important 
since implies that the non-singular s„ can be inductively constructed. 



3.4 Outline of the first order computation 

Having described the recursive procedure for constructing a perturbative solu- 
tion of Einstein equations f l5.65p . in this section we can give an outline of the 
derivation of the metric at first order. 

The zeroth order ansatzfor the metric, G^^j^, has been formulated in fl3.17p . 
The equations that determine G^/^v at are ultralocal, consequently it is sen- 
sible to solve the problem point by point. We can pick a point on the boundary 
which by exploiting the Killing symmetries of the background can be 
chosen to be the origin. At Xq we can use the local scaling symmetry to set 
b^^^ = 1 and pass to a local inertial frame so that l^f*^ = 0. 

Making these choices, the metric (13.1 7p can be expanded to the first order 
in the neighbourhood of the origin 

dsfo) = 2dvdr-r^ /(r) dv^ + dxi dx' -2x'' (9^/3f ^ dx' dr 

- 2 d^pf^ (1 - /(r)) dx' dv - d dv^ . (3.27) 

The metric (I3.27P together with Gfj^^ has a background part (the first three 
terms) which is simply the metric of a uniform black brane, plus a correction 
consisting of first derivatives some of which are known (the last four terms in 

^ These requirements do not yet completely specify the solution for G^]^^, since H posses 
a set of zero modes that satisfy both these requirements. A basis for these zero modes is ob- 
tained by differentiating the d parameters class of solutions p.Sp with respect to parameters 
Pi and b. These zero modes correspond to infinitesimal shifts of /jj"'' and and therefore 
may be absorbed by redefining these quantities [J. 
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(I3.27P ) and the reminder of which are to be determined. 

We now need to make an ansatz for the metric correction at O (e), G^^^, 
which we wish to determine. As was pointed out in the previous section one can 
exploit the SO{d — 1) spatial rotation symmetry at Xq to decompose modes 
into various representations of this symmetry. Modes of G^^jj^ transforming 
under different representations decouple from each other by symmetry. We 
have the following decomposition into SO{d— 1) irreducible representations: 



scalars: G^^ , G[^J, ^ G, 



(1) 



vectors: G^'^^ 



tensors: G^f . (3.28) 

Einstein's equations can be solved sector by sector. In the scalar sector, one 
finds that the constraint equations imply that 

^ d4'^ = , (3.29) 



d-1 

while in the vector sector we have 

These two equations are equations of energy momentum conservation. 

The differential operators entering in the dynamical equations for the func- 
tion J-" defining the terms G^^^ are of the form 

vector : Hd-i-F =—[-— — O 

dr yr'^-^ dr 

tensor : Hd(d+i) ( fir) ^ o] (3.31) 

2 dr \ dr J 

where the index of the operator H is the representation label of the SO{d— 1) 
rotational symmetry. The scalar sector involves some mixing between different 
fields and is slightly more involved. As expected the H operators are simple 
differential operators in the radial variable alone (the form of the differential 
operator remains invariant in the course of the perturbation expansion) and 
can be inverted to find the function J-' once the source s„ is specified. 

Extending this procedure, the calculation can in principle be carried out to 
any desired order in the e expansion. Clearly it is necessary to compute at any 
given order the source terms s^. In addition one always has to ensure that the 
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lower order stress tensor conservation equations are satisfied. For instance, in 

(2) 

order for the source terms which appear in the determination of G)Jj^ to be 
ultra- local at our chosen boundary point Xg, we have to ensure that the first 
order fluid equations of motion are satisfied. 



3.5 Outline of the second order computation 

Once the procedure has been implemented at first order, it is then possible to 
move on to find a solution to Einstein's equation at the second order. Before 
proceeding, it is necessary to pause for a moment to notice that the first order 
computation was performed as a Taylor expansion in the neighbourhood of a 
special point x'^ that we chose as the origin. This was enough to obtain G^lj^, 
while in order to solve the second order problem it is necessary to ensure oneself 
that the first order constraint is solved to the second order in Taylor expansion 
of the fields b^^^ and about the that we chose. Therefore the proper 
requirement is 

dxd^T^^^ix'^) = , (3.32) 

this is equivalent to requiring that T^-q^ satisfies the conservation equation to 
the order e'^ before attempting to find the second order stress tensor. In gen- 
eral the right procedure would have been to satisfy the conservation equation 
globally, however the ultralocality of the set-up allows one to have less strin- 
gent requirements, that amounting to checking that conservation holds only 
at the order at which we are working.^ 

Once we have checked the requirement (13.321) . it is possible to reexpress 
Einstein equations schematically as 



G(o)(/3f),6(o)) G(2)(x^) = + (3.33) 



where for convenience the source term was written as a sum of the pieces. 
Sa is a local functional of /3f'^ and b'^^^ up to the second order in field theory 
derivatives. Terms entering in Sa originated from two field theory derivatives 
acting on the metric G^^^^^ or one field theory derivative acting on G'^]^. 
arises from first order derivatives of the velocity and temperature corrections 
pf''^ and 6*^^^ The second part of the source term is new with respect to the 
first order computation. 



good way to rephrase this requirement is to say that before proceed to the second 
order expansion it is necessary to control that the background defined with the first order 
calculation solve the first order fluid dynamics. 
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Having defined the new source terms, the second order step of the pertur- 
bative calculation proceeds similarly to the first order computation^. 



3.6 Weyl covariance of the metric 

An important aspect of the bulk metric dual to fluid dynamics is that it has 
to transform covariantly under the Weyl transformation [9j. Under the gauge 
choice f l3.22p the bulk metric can be cast in the general form 

ds^ = -2n^(x)da;''(dr + V^(r, x)dx'') + 0^^(r, x)dx^'dx'' (3.34) 

where (3^^ is transverse to the velocity field u'^, i.e., 

u'^'S^, = 0. (3.35) 

All the Greek indices are raised and lowered using the boundary metric gfj_u 
defined by 

9f.u = Ji^^"^ [^M^^ - (3.36) 

and Ufj, is the unit time-like velocity field in the boundary, i.e., g^^u^u^ = — 1. 

Let us consider a bulk diffeomorphism of the form r = e~'^f supplemented 
by a scaling in the temperature b = e'^b, where (p = 4>{x) is defined as a function 
of the boundary co-ordinates. Components of the metric scale as 
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(3.37) 
(3.38) 



and {(3-y = {&-y (3.39) 
dr + Vudx" = e-^{df + V^""). (3.40) 

Therefore we conclude that consistency demands that shall transforms 
like a connection under a Weyl transformation while (5^^ remains invariant. 
These properties of the metric components yield a prescription for constructing 
these two objects as appropriate sums of objects with the right transformation 
properties under Weyl scaling. Therfore (6^^ will be a linear sum of Weyl 
invariant forms, while — rA^ will be a linear sum of Weyl-covariant vectors 
with weight unity. The form of the coefficients is not determined by symmetry 
requirements and have to be fixed by direct calculation. 



detailed derivation of the results for d = 4, 5 at first and second order can be found 
in the original work [?] 



70 



Fluid dynamics from gravity 



3.7 The metric dual to hydrodynamics 

Following the procedure described in the previous sections up to the second 
order, one find the metric dual to the hydrodynamics on the boundary. In this 
section we report the results of the calculation that yield, up to second order, 
a map identifying fluid solutions in the boundary field theory with black brane 
solutions in the bulk geometry. 

The metric can be cast in different forms, using different conventions. For 
the sake of completeness we give in the following a brief overview of some of 
them. 

A very neat and explicit form is 



ds^ = - 2u^dx^ {dr + r A^dx") + [r^Qfxv + U(^^j,S„)\v} — u^^u\^\ dx'^dx" 



+ J^('^^ ~ -ujapuj"'^)ui^u^dx^dx'' 



+ 2{brfF{br) 



-CT^, + F{br) 



- 2{hrf 



dx^'dx" 
{hrY 2{d 



L{hr) px^ 
2{brfHi{br) 



dx^dx"" 



+ 2{brfH2{br) [u^V^a^^ + cj^Va^ + u^^a^x] dxi'dx" , 



dx'^dx'' 
(3.41) 



Where the functions appearing in the metric are defined by 



F{br) 
H,{br) 
H2{br) 



~7~d — TT^^' 
br l/(r - 1) 

00 yd-2 _ I 

br 2/(r - 1) 

d^ /•« 



y'^-My [l + {d- l)yF{y) + 2y''F'{y)\ 



-F{brf 



r dj r 

Jbr m''-'^)Ji 



g /-^ - 1 

y{yd - 1) 



dy. 



Ki{br) = 



1/ dyy'F'iyf, 

br ? 



(3.42) 
(3.43) 
(3.44) 
(3.45) 
(3.46) 
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L{hr) 



ri[l-^(^-l)^'(0-2(rf-ir-^ 

Jbr S 

poo 

+ {2{d-l)i''-{d-2)) I dyy'F'iy) 
dy 



y3^yd _ 1) • 



(3.47) 
(3.48) 



The integral functions have the following asymptotics [5]: 



F{br) 
Hi{hr) 

Kiibr) 



L{hr) 



+ 



# 



hr d{hrY {d + l){hrY+^ (br) 



2d 



+ ... 



1 



+ 



2(6r)2 d{brY {d + 2){brY+^ {brY'^ 
1 1 



oo yd-2 _ ^ 

2(6r)2 d{brY y{y'^ - 1) 
1 2 



2(6r)2 d(d+l)(6r 



+ 



\d+l 



+ 



(c/+ l)(c/ + 2)(6r)'^+2 (5^)2d 



+ ... 



(rf-3)(rf-l) 2(rf-2) 
(i(6r) 



2(rf+ l)(6r)2 
1 



+ 



+ 



# 



rf(2c/- l)(6r)'^ {brY+^ 



+ ... 



+ 



ci(ci + 2)(6r)2 (ci + l)(6r) 
1 



{d + l){2d + l){brY+^ 
1 

2(d+ l)(d + 2)(6r)'='+2 ' ibr) 



2d 



(3.49) 
(3.50) 
(3.51) 

(3.52) 



(3.53) 



(3.54) 



The metric f l3.4ip can be formulated in an equivalent form using the equal- 



ity 



1 



1 



n 



d- 2^^"^^^ + d - 2^^""^^ 2{d -l){d- 2) 



u^ + ... . (3.55) 
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The resulting expression for the metric is 



ds^ = —2u^dx^ (dr + r A^dx'^) + r^g^ydx^dx" 



7^ 



{d-l){d-2) ^ 



dx^dx" + (r^ - -ujaRu"^)unU^dx^dx'' 
[or)"- 2 



+ 2{hrfF{hr) 



dxf'dx" - 2(6r) 



Ki[hr)^ — —P, 



d-1 



a/3 



-2{hrYHi{hr) 



A-n , , „ A, ^ ^ 



+ 2{hrfH2{hr) [u^Vxa^^ + w^^^a^ - a^^uoxu] dx^'dx" . 



(3.56) 



We will instead use the much more compact form (13.341) that renders ex- 
plicit the Weyl symmetries, adopting the Weyl covariant formalism - with a 
slight redefinition of the integral functions - that is 



ds' = Gmn dX^^ dX^ 

= -2 u^{x) dx^ (dr + V^{r, x) dx") + <3^^{r, x) dx^ dx'' 



(3.57) 



where the fields Vfj, and 0^,^ are functions of r and x'^ which admit an expansion 
in the boundary derivatives 



V^ = rAf,- S^x u^-Mbr) Vxo\ + u, 







^r"^ f{br) 



+ \ (1 - f{br)) uj^pu-^ + t)2(6r) 



(3.58) 



r^P^,-UJ^^UJxu + 2{hrfQi{hr) 
-Q2{br) Pf^y-gsibr) 



1 



(y^iu + Qi{br)a^ axu 



(3.59) 



The tensor iS^j, is transverse, since it is built out of operators that are orthog- 
onal to the velocity, and it can be inverted via the relation (0~^)^" = P^- 
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The new integral functions are defined in terms of the old ones as 
Qi{br) = F{br), 02 (&r) = 2 (rbfKiibr), 

Qsibr) = (brfHiibr), Q^{br) = {brfH^ibr), 



t)i(6r) 



-L(6r), X)2{br) 



;K2{br) 



Finally, we recall the following definitions (that are reported in table 2.1) 



(3.60) 



pi/ 



= 2 (t"<^ cu'l 



(3.61) 



3.8 The boundary stress tensor 

The dual stress tensor corresponding to the metric at the second order pre- 
sented in the previous section can be determined quite straightforwardly using 
the holographic prescription in [19| 150 J. To perform the computation one has 
to define a regularisation of the asymptotically AdS^+i spacetime at some cut- 
off hypersurface r = and consider the induced metric on this surface, which 
up to a scale factor involving Ac is our boundary metric g^u- The holographic 
stress tensor is given in terms of the extrinsic curvature K^j^i, and metric data 
of this cut-off hypersurface 

Kpu = gpp V'n, (3.62) 

where is the outward normal to the surface. Now, for asymptotically AdSs 
spacetimes the prescription in [50j is 



lim 



A 



d-2 



167rG' 



(d+i) 



N 



K"" ~ K gf"" - {d ~ 1) g 



d-2 



W --R g^"" 
2 ^ 

(3.63) 

Implementing this procedure, one is able to determine the dual stress tensor 
corresponding to the metric in the previous subsection 



- flU 



P {gpv + du^Uy) - 2r]a^^ 



27]b 



a/3 



-P, 



(3.64) 
(3.65) 
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where 



b 



d 



1 



(3.66) 



47rT' 



P = 





dy. 



(3.67) 



Substituting the standard entropy density, one gets 



V 



s 



1 

4^ ' 



(3.68) 




the celebrated ratio of shear viscosity to entropy density [28] . 



3.9 Horizon location for the dual metric 

A last point to assess is the regularity of the horizon of the black brane metrics 
that we have derived^. The event horizon T-L^ of a given spacetime is defined 
as the boundary of the past lightcone of future null infinity. This is a for- 
mal statement of the physical fact that the spacetime events inside the event 
horizon of the black hole cannot communicate to the asymptotic region. The 
future null infinity X+ is the set of points which are approached asymptotically 
by null rays which can escape to infinity and it is time-like for asymptotically 
AdS spacetimes. Since 'H'^ is the boundary of a causal set, it is a null surface 
which is, in particular, generated by null geodesies in the spacetime. 

The event horizon of the spacetimes dual to the hydrodynamics on the 
boundary is the unique null hypersurface that tends, at late times, to the known 
event horizons of the late time limit of our solutions. As seen above, the metric 
we found can written, in the gauge grr = 0, = —u^, in the form f l3.57p that 
makes explicit the invariance under boundary Weyl transformations. 

Let us suppose that the event horizon is given by the equation 



where is the radial position of the horizon. The vector normal to the 
horizon is defined by the one-form 




(3.69) 



dS = ^Ady^ = dr — d^rndx^ , 



(3.70) 



that can be written in a manifestly Weyl covariant form 




(3.71) 



'^In this section we will follow the derivation proposed in [S] . 
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Also its dual normal vector can be written in Weyl covariant form as 

^^dA = G^^'dASdB = n^{d^ - V^dr) - u^n^dr 

= {d^ + d^r^dr\ = [d,l^^^ . (3.72) 

In the above expression we introduced two new Weyl-covariant vectors k'^ = 
e'^k^ and ra^ = e~^n^ defined via the expressions 

The subscript H indicates that the functions are to be evaluated at the event- 
horizon. 

The induced metric on the horizon l-i^y^x) is defined as 

dsl = [GAB{y)dy^dy\^^^^^^ = 'H,,{x)dx^dx'' , (3.74) 

where the boundary coordinates x'^ have been used as the coordinates on the 
event horizon. Therefore we have 

^-fiu = ^fjiu - U(f,i^u) (3.75) 

and the null-condition on the horizon, [Gab]h(,^(,^ = 'H^u^^^'^ = corre- 
sponds to 

{(&-yK^K, = 2u^K^ . (3.76) 
Now it is possible to construct a Weyl-covariant derivative expansion^ for 



= r(f)+r(f) + ... (3.77) 



®As noted in [5] since there is no first order Weyl-covariant scalar, there are no corrections 
to at the first order in the derivative expansion. A detailed classification of the possible 
Weyl-covariant tensors is contained in [5]. 
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Computing Kn one gets 



2(rf- 1^ 



+ ... 



= u 



1 + (rf - l)6r(^2) ^ _ 62^,^a,^a"^ 



(3.78) 



Now substituting the expression for in fl3.76p . and after some manipu- 
lations, one finds the position of the event horizon as 

1 



(3.79) 



where 



hi 
h2 



2{cP + d-A) K< 



2H 



(P{d-l){d-2) d{d-l) 
d + 2 



(3.80) 



'2d{d-2) 



and 



d{d-l){d-2) 



and K2H is 



= m - i)^'(e) - 2{d - 1)^'-' 

POD 

+2{{d-l)^'-{d-2)) J dyy'F'iyf 



(3.81) 



3.10 Boundary entropy current 

Dissipative fluids are characterized by entropy production. On the other side of 
the duality, for black holes at equilibrium there is a natural deflnition of entropy 
associated with the area of the event horizon which is the usual Bekenstein- 
Hawking entropy of a back hole. Unfortunately, in general, when we consider 
deviations from equilibrium, there is no unambiguous notion of entropy. In- 
deed, the substantial requirement for the entropy current is to satisfy the 
second law 

V^J^ > , (3.82) 
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and it is possible, in principle, to adopt any local function having positive 
divergence to characterize the irreversibility of the fluid dynamical flow. The 
only constraint is that the candidate Boltzmann H-function must agree with 
the thermodynamic notion of entropy in global equilibrium. However since in 
the case of stationary black holes, as discussed in the first chapter, the area of 
the event horizon is associated the entropy of the dual field theory, it seems 
quite natural to associate the entropy of the field theory with the area of the 
event horizon^. 



In the previous section the geometry of a regular future event horizon 
"H"*" for the spacetimes dual to boundary fluid dynamics was derived. It is 
possible, starting from the entropy associated to the area of this horizon, to 
define an entropy current directly for the field theory in the boundary. To this 
end, consider spatial sections of the event horizon, which are co-dimension 
two surfaces in the spacetime, which we label as T-t^. We are working in 
a coordinate chart where the coordinates a* for i = {1, ■ ■ ■ d — 1} define a 
chart on the spatial section and we use as an affine parameter, the boundary 
coordinate v, to propagate these surfaces forward along the horizon generator 
On the surface "H^ it is natural to define a (d — l)-form whose integral gives 
the area of the spatial section. The entropy of the black hole will be - assuming 
that the null energy condition satisfied - proportional to this area. Using the 
Bekenstein-Hawking formula, in terms of the area of the event horizon, the 
entropy density is found to be 

I I , d , , 

(3.83) 



To get the entropy current of the field theory, one has to pull the entropy 
of the black hole back to the boundary. In [12] the authors proposed that one 
pull-back the area form on the horizon using radially ingoing null geodesies, 
which provide an isomorphism between the spatial sections on the boundary 
and the corresponding Ti^ on H"*". The general procedure detailed in ^2] yields 

^In [33] it has been noticed that this point is quite delicate. Indeed, a key feature of 
the event horizon is its "teleological nature", i.e., the fact that the entire future evolution 
of the spacetime is needed in order to determine its location. This may clash with the 
attempts to determine the horizon perturbatively. Associating an entropy current with the 
event horizon may lead to a non-local and acausal definition of entropy [5V. Moreover in 
the case of the conformal soliton flow [5 2) , a simple hydrodynamic flow on R'^-i'i, it was 
shown that the event horizon area does not capture the entropy of the dual field theory 
[53j . Therefore under certain one should use the area of apparent horizons (more precisely 
dynamical horizons), to define the entropy current. 
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the boundary entropy current 



dett\^ 



Js = (3-84) 

[w^ - {(SJ^Tk,] . (3.85) 



dS 



To define det^l^T^ one has to spht the boundary co-ordinates into (f,x*). 
After this split, the components of the n^^ are further split into {n'",n^). We 
denote the d — 1 dimensional induced metric on the constant v submanifolds 
of the event horizon by Finally, one can define 



where Qfj^i, is the boundary metric and the expression on the right hand side 
has been assumed to be pulled back from the horizon to the boundary via the 
ingoing null-geodesics. 

Following this procedure one finds that the entropy current takes the form: 

+ s {Bi Vxa^^ + B2 Vxuj^^) +■■■ (3.87) 

where s is the entropy density and ^1,2,3, B12 are arbitrary numerical co- 
efficients. Requiring positivity of the divergence one finds that two of the 
coefficients must solve the following linear constraint 

Bi + 2^3 = 0. (3.88) 

An interesting result is that the entropy current in d dimensions is a Weyl 
covariant vector of weight d. This is consequence of the fact that the entropy 
density scales like the inverse spatial volume, since the total entropy is dimen- 
sionless, and the velocity field scales according to fl3.40p . Using the geometric 
data of the space-time described by the metric fl3.57p the coefficients are fixed 
to be 
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We conclude this section by observing that, writing the divergence of the grav- 
itational entropy current as 

2 



d 



4:71 T 



;i + Aid)- 



+ 



(3.90) 



which is accurate up to third order in the derivative expansion, it is possible to 
check that it satisfies the requirement of non-negative divergence at this order. 



3.11 The gravity solution dual to hydrodynam- 
ics in the non-relativistic regime 

In chapter 2 we discussed the non-relativistic scaling limit ofthe hydrodynamic 
equations. In particular it was shown that, under this scaling, the incompress- 
ible non-relativistic Navier-Stokes equations are equivalent to the relativistic 
equations of fiuid dynamics dual to gravity up to 0{e^). In this section we 
will report some results about the gravitational dual to a solution of the non- 
relativistic Navier-Stokes equations constructed as a small fiuctuation about 
a black brane background that solves all of Einstein's equations up to order 

In computing the bulk metric up to C(e^), it turns out that only contribu- 
tions from the zeroth order and the first order in the derivative expansion of 
the gravitational solutions of are relevant. Keeping in mind the discussion in 
section 2.7, the metric up to first order in derivatives is 



where 



ds"^ = dsl + ds\ (3.91) 

dsl = —2u^dx^dr + ^_^ u ^Uydx^ dx" + r'^ g ^^dx^ dx" 

dsl = -'^ruu (w^Va) u^dx^dx" 
2 

+ - — -r (V„m") u^Uydx^dx" + 26r^F (6r) a^^^dx^'dx" 
and anu, b, and T can be expanded as 



^ {V^u, + V,u^) (3.92) 

1 - - 1 - 

b = -^ = bo + Sb (3.93) 

T = T0 + 6T. (3.94) 
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V denotes the covariant derivative with respect to the full boundary metric 
Gf^u which is equal to a background 5^^^ plus perturbation if^i., while Tq is the 
temperature of the background black brane. 

Contributions of the metric in which the derivative V appears can be re- 
written as covariant derivatives of the d—1 velocity Vi and the metric pertur- 
bation Ai = with respect to the spatial part of the background metric gij, 
i.e. 

ViUj = ViVj + 0{e'^) 

V^uo + VoUi = do{vi + Ai) - -dihoo- -di{vjV^) -v^Fij + 0{e^) 

VX = V,V + 0(e^) 
m'^V^Mo = 0{e^) (3.95) 

M^V^Mi = do{vi + Ai)-^dihoo + {v^Vj)vi-v^Fi, + 0{e^) 

Fij = diAj djAi. 

Now V is the covariant derivative with respect to the background metric Qij. 
The raising and lowering of the i,j indices is to be performed using the metric 
Qij. To simplify the expression of 0"^,^ in f l3.96p the constraint Vjf* = has 
been used. 

Substituting in the metric, the first order part takes the form 
dsj = bor^F{bor) (ViVj + VjVi) dx'dx^ - 2bor^F{bory (ViVj + VjVi) dt dx' 

+ 2r (^do{v, + Ai) - ^dihoo - v^F,, + iv^Vj)v,^ dt dx' + O(e^) 

(3.96) 

where the first term is of order and the last two terms are of order e^. The 
constraint equation Vjf* = 0, was used to cancel out contributions from the 
scalar sector. 

Also the zeroth order metric can be expanded in powers of e. One finds 
that to solve Einstein equation up to order e^, it is sufficient to expand the 
zeroth order metric up to order in the fluctuations. The result is 

dsl = yf-p^dt^ + [-dt^ + Qijdx'dx^) + 2dt dr 
2 

- Trf-^ {A + Vi) dt dx' -2{Ai + Vi) dx'dr 

o"" ^ (3.97) 

H — {Ai + Vi) {Aj + Vj) dx^dx^ — {—VjV^ + /iqo) dt dr 
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where the first hne is of order e°, the second line is of order and the third 
and the fourth are of order e^. 

It is possible to show by a direct calculation that the proposed metric 
solves Einsteins equations up to the third order in e, provided that the velocity 
and temperature fields above obey the incompressible Navier-Stokes equations 



An important observation arises from the equation that determines 5T (and 
hence Sb) 



ST is a spatially nonlocal but temporally ultralocal functional of the velocity 
fields f Therefore, even though the bulk metric at x'^ is determined locally as 
a function of temperatures and velocities at x^, it is not determined locally as 
a function of velocities at x^. This nonlocality is a consequence of the infinite 
speed of sound in the scaling limit. 




To 
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Chapter 4 

Matter and Gravity 



In this chapter we will discuss some techniques and issues that are relevant to 
the study of black holes and singular hypersurfaces in gravity. 

In particular we will start by reviewing various possible energy conditions 
that enter in black holes physics [Ml 155] . Then we discuss the issue of the 
violation of the loosest among these condition, the null energy condition, and 
its physical consequences [561 EZl EH] • 

Later on we will introduce some geometric notions useful in the description 
of hypersurfaces and in particular the extrinsic curvature. These concepts will 
be used in applying Israel's seminal work on singular hypersurfaces [TT] to the 
configuration that we will present in the next chapter. 

4.1 Energy conditions 

In general relativity, an energy condition is one of various alternative conditions 
which can be applied to the matter content to define some of its properties. 
Since almost every spacetime is a possible solution to Einstein's equations 
for some particular choice of energy- momentum tensor T^j,, it is important to 
understand what general restrictions can hold on the energy-momentum of a 
physical system and which are their implications. Indeed, energy conditions 
play a critical role in the study of black holes and cosmological singularities, 
entering in important theorems such as the no hair theorem and the laws of 
black hole thermodynamics [SlJ [55]. Moreover, even though a violation of an 
energy condition does not necessarily imply a pathology in a physical system, 
the violation of the null energy condition has been related to the insurgence 
of instabilities and superluminal propagation (one can see e.g. [56[ [57 [ 158]). 
Let us start with the weak energy condition. 
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Weak energy condition: The energy-momentum tensor at each point p 
on a manifold J^, p E Ai, obeys the inequality 



T^uwV > , 



(4.1) 



for any timelike vector w in the tangent space at p, w E Tp. By continuity this 
will also be true for any null vector w & Tp. 

The weak energy condition is equivalent to the statement that the energy 
density of any matter distribution, as measured by any observer in spacetime, 
must be nonnegative. Indeed, for an observer whose world line at p has unit 
tangent V, the local energy density is Tp^yV^V" . 

Following [M] , we will explore more in depth the meaning of this assump- 
tion using the fact that it is possible to cast the 4-dimensional stress tensor Tp^ 
in four canonical forms, with respect to an orthonormal basis {eo, ei, e2, 63}. 



Type I 



rpfj.1/ 



p 














Pi 














P2 














P3 



(4.2) 



In this case, the stress tensor has a unique timelike eigenvector (unless p = —pi, 
(z = 1,2,3)). p is the energy-density measured by an observer whose world 
line has unit tangent vector eo at p, while the three eigenvectors Pa are the 
pressures in the three spacelike direction spanned by Si {i = 1,2,3)^. The 
energy-momentum tensor can be cast in this form for fields with non-zero rest 
mass and for some zero rest mass fields, except those of the Type II. 



Type II 



rpfMU 



U + K 


V 








V 


V — K 














P\ 














P2 



±1 



(4.3) 



In this case, the stress tensor has a double null eigenvector (eo + ei). The 
stress tensor has this form in the case of radiation traveling in the direction 
Go + Gi. In this case pi, pi and k are zero. 



^ As discussed in chapter 2, if the stress tensor is that of a perfect fluid, then the pressures 
in the spacehke directions are equal V\ — Vi — Vz — V- 
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lype III 



rpfJ,L> 



V 














—u 


1 


1 





1 










1 





p 



(4.4) 



The stress tensor of this type has a triple null eigenvector (eo + ei). Actually 
there is no observed field of this kind. 



Type IV 














V 


—K 














Pi 














P2 



(4.5) 



This stress tensor corresponds to the case in which there is no timelike or null 
eigenvector. No field of this type has been observed. 

The week energy condition holds for type I under conditions 



p>0 p + Pi>0 for i = l,2,3; 

and for type II if 

Pi>0, P2>0, K>0 , = 
while it does not hold for type III and IV. 



(4.6) 



(4.7) 



Dominant energy condition (A): The energy-momentum tensor at 
each p & A4, obeys the inequality 

Tf,,wV > , (4.8) 

for any future directed timelike vector w e Tp, and —Tj^w" is a future directed 
non-spacelike vector. 

The dominant energy condition embodies the notion that matter should 
fiow along timelike or null worldlines. It may also be interpreted stating that 
to any observer the local energy density should appear nonnegative and the 
local energy fiow should be non-spacelike. It is possible to restate the domi- 
nant energy condition as: 



Dominant energy condition (B): In any orthonormal basis the energy 
dominates the other components of the stress tensor 
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This condition holds for type I if 

p>0 -p<Pi<p for i = l,2,3; (4.10) 

while it holds for type II under the conditions 

u^+l, K>0, 0<Pi<K for i = l,2; (4.11) 

and again it does not hold for type 111 and IV. Therefore the dominant energy 
condition, in addition to requirement of the weak energy condition, states that 
the pressure should not exceed the energy density. 

Strong energy condition: The energy-momentum tensor at eachp e A4, 
obeys the inequality 

(t^, - ^Tg^,^ wi'w^ > , (4.12) 
for any future directed normailsed timelike vector w E Tp. 
Since the Einstein equation implies 

T^. - ^Tg^u = (4.13) 

the strong energy condition is a statement about the Ricci tensor. 
The strong energy condition holds in type 1 for 

3 

P+'^Pi>0, P + Pi>0 for i = l,2,3, (4.14) 

i=l 

and for Type II if 

i/ = +1, K>0, Pi > 0, P2 > . (4.15) 
It is useful to notice that this condition does not imply the weak one. 

In presence of a cosmological constant term, the strong energy condition is 
known as the null convergence condition. In this case the statement is simply 
modified to 

(t^, - ^Tg^, + ^^g^^.^ w^w" > ■ (4.16) 
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The conditions for type one now are 

P + J2p'-^^-^ ^ P + Pi^^ for i = 1,2,3, (4.17) 

i=l 

and for Type II read 

i/ = +l, Av>0, Pi>0, p2>0, pi+p2-i^A>0 . (4.18) 

47r 

The weakest requirements on the stress tensor are formulated in the null 
energy condition. 

Null energy condition (NEC): The energy-momentum tensor at each 
p G Ai, obeys the inequality 

T^uk^k" > , (4.19) 
for any future directed null vector k & Tp. 

The null energy condition holds for type I under the condition 

p + Pi>0 for i = 1,2,3, (4.20) 

and for type II if 

= + ^ + >o i = l,2. (4.21) 

It is easy to convince oneself that the null energy condition is implied 
by both strong and weak conditions, while the weak one is implied by the 
dominant energy condition. Relations among different energy conditions are 
summarised in the following diagram: 





Strong E.G. ' 


> =^ Null E.G. 


Dominant E.G. = 


Weak E.G. 





In table 4.1 we summarise energy conditions and their associated con- 
straints on type I stress tensors. The discussion above - that we carried on in 
4 dimensions - can be easily generalised to higher dimensional spaces. 
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Name Statement Conditions for Type I 

Weak Tf.^w'^w'' >0 p>0,p + pi>0 

Dominant —T}^w^ future directed P > 0, |p| > pi 

Strong {T^u-\Tg^,)w^'w^ >Q P + ELi > 0, p + Pi > 

Null T^^A;^r > P + P-. >0 

Table 4.1: Energy conditions. 



4.2 Violation of the null energy condition 

Energy conditions are typically fulfilled by classical fields. However it is well 
known that they can be violated by quantized matter fields. The most com- 
mon example is the Casimir vacuum energy between two conduction plates 
separated by a distance 5 

"=-720?- f^-^^' 



In this kind of situation ione usually formulates an averaged version of 
the energy conditions. Indeed, for each of the energy conditions listed in the 
previous section it is possible to formulate a corresponding averaged version. 
Average energy condition requirements of non- negativity must hold only on 
average along the flowlines of the appropriate vector fields. For example the 
averaged null energy condition states that the integral along a null geodesic 7 
of T^yk^k^ must be non-negative 

T^^k^k'^dX > . (4.23) 

This condition holds for noninteracting scalar and electromagnetic fields in 
arbitrary quantum states; this is true even though T^j^^k'^k" may be locally 
negative. 

Even though, in some cases energy conditions are violated {e.g. the strong 
energy condition may be violated in some inflation scenarios), it is generally 
believed that any well behaved physical systems should respect at least the 
null energy condition. In various settings it has been been shown that the 
violation of the null energy condition results in the insurgence of instabilities 
in the system sourcing T^i, [561 EI] • 
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For example, in [SZ] the authors showed that classical solutions of both 
minimally and non-minimally coupled scalar-gauge models which violate the 
null energy condition are unstable. Moreover it was also proven that perfect 
fluids which violate the NEC are unstable. Finally quantum states in which 
the expectation of the energy-momentum tensor violates the NEC cannot be 
the ground state, including models with fermions. On the basis of these re- 
sults, it seems that physically interesting cases of violation of the null energy 
condition are likely to be ephemeral. 

In [58] the authors studied systems of derivatively coupled scalar fields with 
coordinate dependent condensates, without taking into account gravitational 
back-reactions. In this setting it was suggested that a violation of the null 
energy condition does not strictly imply instabilities - such as ghosts or imag- 
inary frequencies - in the system. Indeed, it was proven that there may exist 
an acceptable effective field theory with a positive definite Hamiltonian for 
quadratic perturbations around a background whose energy momentum vio- 
lates the null energy condition. It turns out that a necessary feature of such 
a model is the anisotropy of the background and the presence of superlumi- 
nal modes. More generally it was suggested that for systems that are either 
isotropic or do not feature superluminality, a violation of the null energy con- 
dition always implies an unescapable instability. In the next chapter we will be 
interested in two nonisotropic backgrounds which both violate the (averaged) 
null energy condition on a hypersurface. 



4.3 Some notions on hypersurfaces 



Let us introduce some notation and some concepts related to the study of 
hypersurfaces in General Relativity |59j. We will introduce concepts in four 
dimensions. These can be readily generalised to higher dimensions. 

Consider a curve parametrised by a parameter r on a four- dimensional 
Riemannian manifold Ai of class C^. The covariant derivative of a smooth 
vector function A'^ defined on this curve will be denoted (we adopt the notation 
of [33j) as: 

DAf'(T) dA^ dx" 

(4. 24) 



Dt 



dr 



+ ^ dr 



If A^^iu^v) is defined on a two-space 
of the Riemann tensor spells out 



x^(M,f), the standard definition 



D D 



A^ 



. „ dx^ dx'^ 
^^'^^ du dv ' 



(4.25) 
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Now, let us consider a smooth hypersurface T, in M. with unit normal 

. s r +1 spacelike n . . _„s 

n-n-e(n)^^_^ ^.^^^^^^ . (4.26) 

An infinitesimal displacement on the hypersurface can be expressed in 
terms of the intrinsic coordinates ^* as 

ds = e(,)de , (4.27) 

where e(j) are a natural set of independent tangent vectors associated with the 
intrinsic coordinates In terms of the system of coordinates {x'^} defined on 
V, the vector e(i) can be expressed as 



= (4-28) 



In what follows we will write explicitly 4-dimensional indices only when nec- 
essary, otherwise leaving them implicit to avoid confusion. 

Given a tangent vector function lying in the hypersurface, it is possible to 
define intrinsic components as 

A^A'- e(,) , (4.29) 

The scalar product of this vector with the base vector e(j) is 

(A ■ e(,)) = A\e^i) ■ e(,-)) = A^g^j = Aj , (4.30) 

where the metric three tensor Qij was introduced. 

The covariant derivative of the vector A in the direction e(i) is 

(^)Ve(,^ =(^) V.A =(^) Ve,,{e,Ar) = e,^ + {^'%e,)A^ (4.31) 

where we used the notation to make evident the dimension of the space- 
time in which objects are defined. Indeed, the above definition has components 

"out of the hypersurface" . 

To get a covariant derivative defined intrinsically to the 3-dimensional hy- 
persurface it is necessary to project the above derivative orthogonally onto S. 
In this way one gets 

An\i = en -^'^ Ve(,/ = -Q^- ^""^mM . (4.32) 
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The connection Tjji,hi in three dimensions can be expressed in terms of the 
basis vector and of the metric 3-tensor Qij = e{i) ■ e(j) 

^'^r^,hi = e(„) -(3) V,ei^H) ■ (4.33) 

Given this definition, it is evident that intrinsic covariant derivatives and the 
associated intrinsic Riemann 3-tensor R^j^n do not depend on the nature of 
the embedding and are therefore invariant under changes of the embedding 
space that preserve the intrinsic metric on S. 

We are now in a position to define the notion of extrinsic curvature, that 
can be seen as a description of the curvature of a shce of the 3-geometry relative 
to that of the 4-geometry. This is measured by the covariant variation ^^^Vj^ 
of the unit normaL In components the extrinsic curvature is defined^ by 

(^V,n = ir/e(,). (4.34) 

Taking the scalar product of this definition with the basis vector e(m,) and 
remembering that 



Hm) ■n = 



by definition, we get 



K^m = K'gjm = em ■ ^ = -n ■ ^ = -n ■ ^ = ir,, , (4.35) 

that spells out the symmetry properties of Kij. 

The Gauss- Weingarten equations can be derived directly form fl4.33p and 

^ = -e{n)K,,n + r5,.e(.) . (4.36) 

It is possible to work out a different expression relating the extrinsic cur- 
vature to the standard Riemann tensor. Particularly useful are the equations 
of Gauss and Codazzi: 

^'^^R^p.se'l^^elflfl,) =^'^ Rabcd + K^cKm - K,,Ka, , (4.37) 
^"^^R^fs^sn'^el.flfl^) = Km\c - i^6c|d . (4.38) 

Equations fl4.38p can be contracted using respectively g^'^g"-'^ and g^'^ and 
observing that 

9%fl) = 9^' - e{n)n^n^ . (4.39) 



^In the definition of Kij we adopt tlie sign convention of [TT|, which is opposite to that 



in [Sg 
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The contracted versions of the equations of Gauss and Codazzi are 

-2e(n)(^)G'„^n"n^ =(3) R + KabK^^ - , (4.40) 
^'^G^pe-^^^n^ = Kl^, - , (4.41) 

where K = g'^^Kab and all terms are independent of the coordinates. 

Consider now a 3-tensor field Sab living on the hypersurface S, it is possible 
to define an associated 4-dimensional discontinuous vector 4-tensor as 

10 off S ' ^ ' ' 



then we have 



Vm5'"'' = eL^^V - ^(^)'5"'^abn" (4.43) 



where is the standard covariant derivative defined on the coordinates x^. 



4.4 Singular hypersurfaces 
in General Relativity 

An important issue in gravitational theory is the formulation of correct junc- 
tion conditions at surfaces of discontinuity. These surfaces can be characterised 
by a jump discontinuity in the density, such as boundary layers or shock waves; 
or by the fact that the density become infinity, as happens for surface layers. 
In Newtonian gravitation, one has to set a system of coordinates which has to 
be a priori well defined, supplemented by the appropriate continuity and jump 
conditions connecting the potential and its first derivatives across the surface. 
In Einstein's theory of gravitation, the issue is more complicated. Indeed, the 
smoothness of the gravitational potential Qajs is determined by the smooth- 
ness of the physical conditions and by the smoothness of the coordinates one 
chooses to describe the space-time manifold [TT]. The problem is therefore 
made more subtle by the necessity of distinguishing between the physical dis- 
continuity and spurious "bumps" that may arise from an unhappy choice of 
coordinates. 



Let us define in a more precise form the classification of surfaces of dis- 
continuity. Let us focus on a hypersurface S which separates the spacetime 
in two four-dimensional manifolds V~ and V~^. Moreover we assume that V~ 
and are both of class and contain E as part of their boundaries. The 
normal n to E, that we can choose to be directed from to V^, is assumed 
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to be spacelike everywhere. Let K~j and be the two extrinsic curvatures 
of S associated respectively with embeddings V~ and V~^. 

Given extrinsic curvatures, we can give a definition of singular hypersur- 
faces independent of the choice of four-dimensional coordinates. We define, 
singular hypersurfaces of higher order surfaces E for which everywhere 

Kr, = K± , (4.44) 

and singular hypersurfaces of order one or histories of a surface layer as those 
surfaces S for which 

Ki^^K^- (4.45) 

In particular, (the histories of) boundary layers are part of the class of the 
singular hypersurfaces of higher order that contain also all of the regular hy- 
persurfaces. 

The case of boundary surfaces can be rapidly assessed. From eq. fl4.4ip we 
find 

G^/jn^n^P = G^pn^nf'\' , G^^^ef^^n^f = G,;3ef,)n^|+ , (4.46) 

where we considered the limits of the function Gapn'^nl^ and Gape'^^-^nl^ at 
E from and V~ , and we used the fact that the expressions fl4.4ip have 
a continuous right hand side across E. Equations (I4.46P establish identities 
between terms that may be evaluated distinctly in two independently chosen 
4-dimensional coordinate systems on V~ and . It is interesting to notice 
that there is no requirement about the continuous matching of charts x"^ and 
on E. The only requirement about the two charts is that, given the set of 
intrinsic coordinates the equations of E 

= f{i\ e, e) , < = g'^ie, e), (4.47) 

must be formulated in terms of two functions g". 

The case of a surface layer is more interesting. We start by defining the 
non-vanishing 3-tensor 

7., =K^-K^ . (4.48) 

In terms of jij it is possible to define the surface energy tensor of the layer, 
as^ 

SnSij = --fij + g^j-i^ . (4.49) 

■^As before we set the Newton constant and the speed of hght in vacuum equal to one, 
G = c = 1. 
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Eq. (14.49 1) can be reformulated as 

7ij = -Svr (^Sij - ^gtjS^ , (4.50) 

where 5" = S/f'. 

For a surface layer in vacuo Gq/? = both in V~ and V~^, therefore equa- 
tions fl4.4ip are reduced to constraints on the extrinsic curvature: 



^^^R + K^h^f -Kl = , (4.51) 

Ma 



After some manipulations eq fl4.52p can be recast as the following group of 
constraints: 

5""^ = , (4.53) 

-K\a = 0, (4.54) 

KabS""' = , (4.55) 

^'^R + KabK'^' -K' = -IGn^SabS"' - ^S^) , (4.56) 

where 

KabliKbKb). K = g'^'K^,. (4.57) 

Given two independent charts on V~ and , respectively and x", we 

can define, using eq. (14.421) . 4-dimensional extensions S*^^ of Sab in the two 
domains V^. Hence, we can reformulate condition (14.531) using relation fl4.43p 
as 

e(a)V,5/|^ = , n-W,Si\^ = -K^^S^' . (4.58) 
Finally, using equations fl4.50p and fl4.55p we find 



(4.59) 



Stt ( SabS'^' - ) , (4.60) 



that can be seen as the relativistic analogues of the Newtonian formula for the 
normal components of mechanical force due to self-attraction of the two faces 
of a plane layer of surface density a: 

aF' ■ n = -aF+ ■ n = 2a^ . (4.61) 
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The tangential force, on the contrary cancels out, in agreement with the first 
equation of fl4.58p . 

We can conclude this section by giving a heuristic argument for the def- 
inition of the surface energy tensor (14.49 p . Let us consider a layer of finite 
thickness e separating the two vacuum regions V~ and and having respec- 
tive boundaries S~ and S+. It is useful to introduce Gaussian coordinates x", 
based on E~ by setting x* = ^\ and letting x° equal the geodesic distance 
normal to S~, taken with sign plus or minus respectively for points in and 
V~ . With these coordinate, the equations for are respectively = and 

= e. Moreover, for a 3-space x^ = cost., one finds the following expression 
for the extrinsic curvature 

K, - . (4.62) 



while 



where Zij is defined as 



'''R^. = ^ + % , (4.63) 



% =(3) i?,^. - KK,, + 2K^Khj . (4.64) 
At this point it is possible integrate Einstein's field equation over the layer 

Rap = -Stt (t^i^ - ^gafsT^ ■ (4.65) 

Remembering expression fl4.63p one gets 



- Stt ^ (^T^p - ^ga^T^ dx' = K± - R-- + Z^jdx' 



(4.66) 



In the limit £ — )■ 0, for given values of K^, if Kij remains bounded in the layer 
then the integral of Zij over the layer will converge to zero. Hence Sij defined 
in flCTjl is 

Sij = lim / Tijdx^ , (4.67) 
Jo 

or in other words, Sij is the integral of the stress tensor Tij over the surface 
layer. 
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Chapter 5 



The Surface Layers Dual to 
Hydrodynamic Boundaries 

The AdS/hydrodynamics correspondence provides a 1-1 map between large 
wavelength features of AdS black branes and conformal fluid flows, so long as 
the fluid variables are constant over distances large compared with the inverse 
temperature. For example progress towards gravity duals of shock waves and 
vortices has appeared in Refs. [60] and [61]. A case of particular interest is 
represented by gravity duals to turbulent flows. Turbulence is generic in fluid 
flows under a wide range of conditions. The dual of these fluid conditions then 
provides some condition on a gravity solution under which it to generically 
decays into a turbulent conflguration. An example of such a situation was 
presented in Ref. [9]. It would of course be interesting to characterize the 
gravity duals of turbulent flows, and of the conditions under which turbulence 
may be expected. In hydrodynamics, even the most basic scaling laws are 
altered by turbulence. If gravitational solutions near, for example, spacelike 
singularities (where indeed chaotic evolution is expected [7j) or certain event 
horizons do generically decay to turbulent solutions, it would be difficult to 
overstate the potential consequences for, for example, the horizon problem. 

Perhaps the best understood turbulence is steady state turbulence, in 
which energy is injected into a system at the same rate at which it dissipates. 
Richardson's cascade model [8] of steady (3-|-l)-dimensional turbulence is as 
follows. Energy is injected into a system at large characteristic distance scales, 
for example, a lake warms the air. This creates large vortices, which decay 
into smaller vortices. Thus the energy flows to smaller distance scales. At suf- 
flciently small distance scales, higher order derivative terms in the equations 
of motion become relevant, such as viscosity terms. These lead to dissipation 
of the energy in sufficiently small vortices. Thus energy cascades from the long 
length scale in which it is introduced, down to the dissipation scale. 
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To realize steady state turbulence, one needs to inject energy into a sys- 
tem. There are two principal ways to do this. First, one may deform the 
fluid via external perturbations. Second, one may apply boundary conditions, 
for example one may consider fluid flow in a pipe or wind tunnel. The flrst 
approach was applied to the AdS/hydrodynamics correspondence in Ref. [D], 
where it was argued that a laminar fluid flow and the dual gravity solution 
decay to turbulent conflgurations. This approach has the disadvantage that 
solutions are quite complicated, due to the necessarily inhomogeneous forcing 
and to the geometric implementation of the forcing on the gravity side. 

In what follows we will take a preliminary step towards a realization of 
the second approach to creating steady state turbulence, we will investigate 
boundary conditions in the AdS/hydrodynamic correspondence. For simplic- 
ity, we will consider nonrelativistic, incompressible flows. Consider the surface 
which separates a solid object from such a fluid. The normal velocity of the 
fluid into the solid must vanish. If furthermore the fluid is viscous, as fluids in 
the AdS/hydrodynamics correspondence are [2], then the tangential relative 
velocity of the fluid must also vanish. 

What does this correspond to on the gravity side? The answer to this 
question is not necessarily unique, one may deflne a dual and then attempt to 
understand its dynamics. One interesting case, which is already suflicient to 
generate turbulence, is a solid which is a thin, inflnite sheet with a stationary 
fluid on the left side and a moving fluid on the right. In this case a natural 
choice would be to consider the gravity duals of both fluids and then to attempt 
to glue them together. Equivalently one may choose to think of the entirety 
of the left side as a solid wall, fllling the left half of spacetime, and a liquid 
flUing the right half. The wall is stationary and so one chooses the dual to 
be a stationary black brane in half of AdS. Whatever one chooses to think, 
the logic is that one imposes that the left half of the gravity dual be a static 
black brane in AdS, and that the right side be the gravity dual given by the 
prescription of Ref. [1] . 

So how does one glue these two vacuum gravity solutions together? Clearly 
there are many inequivalent choices. One possibility is to simply attach them 
and then use the Israel matching conditions [llj to determine the stress tensor 
on the surface layer that separates the two sides. This is equivalent to letting 
the gravitational solution continuously interpolate between the two solutions 
over a flnite distance d and then taking the limit as this distance tends to zero. 
While there are many ways of performing this interpolation, so long as the 
extrinsic curvature is kept flnite, they all lead to the same stress tensor as the 
interpolation distance c? — )■ 0. 

Another possibility is to let the fluid conflguration continuously interpolate 
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between the two solutions, and then take the dual using the prescription of 
Ref. [1]. As the fluid is not a solution of the Navier-Stokes equation in this 
region, the dual will not be a solution of the vacuum Einstein equations in this 
region. Instead it will solve Einstein's equations with a nonvanishing stress 
tensor supported on a surface layer. The ultralocality of the duality map 
implies that the vacuum Einstein equations will however be solved way from 
the surface layer. In this case, one cannot take the interpolation distance d to 
zero, because the dual is not defined when derivatives are large with respect to 
the inverse of the temperature T. Thus the minimum size of d will be of order 
1/T. Again there are many inequivalent ways of performing the interpolation. 
But we will see that, at least for the quantities at we are able to calculate, 
when d is large with respect to 1/T, the difference between these prescriptions 
is suppressed by powers of dT and so, like Israel's method, there is a single 
answer. 

The perhaps surprising result is that the two methods yield bulk stress ten- 
sors which differ by a finite amount. They did not need to agree, indeed one is 
derived at small d and the other for large d. The reason that they disagree is as 
follows. The construction of the metric from the fluid flow proceeds order by 
order in the derivatives of the fluid's velocity. The boundary conditions imply 
that the velocity of the fluid is the same on both sides of the wall, however the 
first derivatives differ. Therefore, whatever regularization scheme one uses on 
the fluid side, the second derivative of the velocity diverges at small d. This 
means that the metric corrections derived using the map of [1] will diverge at 
small (i, invalidating the flniteness assumption in Israel's derivation. In fact, 
we will see that the disagreement between the two calculations of the stress 
tensor differ only in these higher derivative terms. Of course the fluid map is 
not deflned at small d, as it yields a divergent series, and so no divergences 
appear within the range of validity of either approach. 

We will begin in Sec. 15.11 by describing the flow of interest. The velocity 
will be kept sufficiently arbitrary to allow a general interpolation between the 
fiows on the two sides of the wall, and in Sec. 15.21 the naive gravity dual will 
be calculated using the prescription of [3] . We will see that those higher order 
derivative corrections which we calculate are indeed suppressed by factors of 
dT. Then in Sec. 15.31 we will calculate the bulk stress tensor of the interpo- 
lation between the two gravity solutions. First it will be calculated for the 
interpolation dual to a continuously interpolating fluid flow. It will be seen 
that contributions from the second derivative of the velocity are d-independent, 
while higher order contributions are suppressed by powers of dT. Thus the re- 
sult is independent of the interpolation scheme when d is sufficiently large. The 
stress tensor will then be calculated directly from the Israel matching condi- 
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tions on the two solutions of the vacuum Einstein equations. It will be seen 
that the two stress tensors agree up to terms corresponding to a divergence 
in the extrinsic curvature at small rf, and that only the second stress tensor 
contains a nonvanishing stress. 



5.1 The Flow 
5.1.1 The ansatz 

We will consider a hydrodynamic flow in 4-dimensional Minkowski space, using 
a (— , +, +, +) metric. To highlight the essential features of the boundary 
condition, we will consider the simplest possible flow. The liquid will only 
move in the y direction, with a velocity v = v{x) that only depends on the 
coordinate x. The velocity will be taken to be small, and we will drop all terms 
which are quadratic in v. In fact, as described in Refs. [3Z1 EI] we will work 
in the nonrelativistic, incompressible limit. More precisely, we will show that 
our flow satisfles both the full relativistic equations of motion at order 0{v) 
and also the incompressible Navier-Stokes equation. 

We will set c = 1. The conformal fluid which is dual to Einstein gravity 
with a negative cosmological constant is very particular. Being conformal, all 
of its transport coeflicients may be expressed in terms of a single dimensionful 
quantity, such as the temperature T, and certain constants which may be 
calculated from the gravity dual. In the case at hand for example the shear 
viscosity t], pressure p and density p have been found in Ref. [1] 

(^-^ 

where Gat is the dual Newton's constant. 

The relativistic velocity 4-vector u is, to linear order in v, simply 

u, = (^=^,0,-=^,0) ~ {1,0, v,0). (5.2) 
V 1 — f V 1 — 

We will be interested in the fluid velocity in three regions, as illustrated in 
Fig. 15.11 First, on the left, where v = 0. Second, we will be interested in the 
velocity on the right, where v will be linear in x. We will show momentarily 
that this is a solution to the hydrodynamic equations of motion and so will be 
dual to a vacuum solution of Einstein's equations. Finally, we will be interested 
in an interpolating region where v will be arbitrary and we will not impose 
the equations of motion, therefore the dual metric will not solve the vacuum 
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v(x) 





x=0 



x=d 



X 



Figure 5.1: The fluid velocity v is in the y direction, and it depends on the x 
coordinate. On the left the fluid is stationary, on the right the fluid velocity 
is linear. These two regions solve the fluid equations of motion at linear order 
in V. There is an interpolating region of width d, which must be larger than 
the inverse temperature, in which v does not satisfy the equations of motion. 
V and its first derivative v' are continuous at x = and x = d. 



Einstein equations but, like any metric, will solve Einstein's equations with 
some stress tensor. 

Clearly the left, v = 0, satisfies the fluid equations of motion. We will now 
verify that the region on the right satisfies the relativistic equations of motion, 
which are simply the conservation of the stress tensor 



In accordance with the usual fluid approximation [3T], we will work at large 
enough distance scales that only the velocity v and its first derivative v' need 
be considered in the stress tensor. This approximation in general is problem- 
atic, leading for example to superluminal propagation [52]. However, as we 
will be interested in velocities well below the speed of light, no problems will 
arise. Later, when we will consider the interpolating region, where the second 
derivative may be large, we will make no such approximation. We will con- 
sider the bulk stress tensor to higher order, calculating all terms up to two 
derivatives and several terms up to three or four derivatives to check that they 
are subdominant. However we do not impose that the interpolating region 
satisfies the equations of motion, indeed that would lead to a vanishing bulk 
stress tensor. 



= d^T^^r 



(5.3) 
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5.1.2 Relativist ic and nonrelativistic equations of mo- 
tion 

As discussed in Chapter 2 (see Sec. 12.31) . dropping all higher derivatives of the 
velocity and using the fact that the fluid is conformal to eliminate the bulk 
viscosity and replace p with 3p, the hydrodynamic stress tensor is 



where the shear strain rate a'^'^ is defined as 



2?7(T' 



(5.4) 



(5.5) 



Here parenthesis denote, as usually, symmetrization with a factor of one half 
and P^'^ is the projector onto the spacelike directions in the reference frame of 
the fiuid 
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Substituting the velocity ansatz (15. 2p into the definition (15. 5p one easily 
finds the shear strain at linear order in v 
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(5.7) 



The constants of proportionality (15.11) in this particular fiuid can then be 
inserted into the general formula (15.41) for T'^" to express the stress tensor in 
terms of the temperature T and the velocity v 
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The velocity only depends on the coordinate x. Let us choose boundary 
conditions so that the temperature T also only depends on x. Then the equa- 
tions of motion (15. 3p are simply 



= d^T""". 



(5.9) 
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However all of the components T^'^ are constants except for T^^ and T^y. 
Therefore the only nontrivial equation of motion at linear order in v is 



When V is linear in x, its second derivative vanishes. Therefore this equation 
of motion, the conservation of momentum in the y direction, implies that the 
temperature is constant. More precisely it implies that {dT)/T is negligible in 
this approximation. In light of the relations (15.1 p this is consistent with the 
incompressibility assumption that we have imposed on our fluid. 

Therefore we recover the fact that at sufficiently small velocities, a constant 
temperature and a velocity which is linear in x solve the equations of motion 
(15. 3p . Intuitively this is clear. Further to the right, the fluid is moving faster. 
Therefore the viscous force on a unit of fluid exerted by the faster fluid on 
its right (in the +x direction) will accelerate it in the +y direction, whereas 
the slower fluid on its left will exert a viscous force that decelerates it. A 
steady flow occurs when these two forces cancel, which implies that the second 
derivative of the flow vanishes. Had there been a temperature gradient, then 
the viscosity to density ratio would also have been stronger on one side by 
(15. ip . and so this balance could only be maintained by introducing a second 
derivative of the velocity. 

Clearly a linear velocity also satisfles the nonrelativistic Navier-Stokes equa- 
tion for an incompressible, Newtonian fluid 



In fact, each term vanishes independently. Note that the vanishing of the dp 
term is not merely a consequence of incompressibility. In incompressible flows 
it may be of the same order as the viscous term. It vanishes in this case 
because this provides a solution to (15. lip and it is consistent with the various 
nonrelativistic, small gradient and incompressible limits taken above. 

In conclusion, we have considered fluid flows with an x-dependent velocity 
V in the y direction. We have verifled that, to linear order in v, these satisfy 
both the relativistic and nonrelativistic equations of motion when v is linear in 
X and the temperature T is constant. Our flows of interest will have f = on 
the left, V linear on the right and an interpolating region inbetween. Thus the 
equations of motion will be satisfled on the left and the right but not in the 
interpolating region, leading to a dual gravity solution which solves Einstein's 
vacuum equations on the left and right, but inbetween requires material de- 
scribed by a nontrivial stress tensor. We have also found formula (15. 6p and 
(15. 7p for the projector P^'" and the shear tensor a'^^ for general functions v, 
and so these results may be applied to the interpolating region. 




(5.10) 
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Clearly if the fluid velocity is linear over a large enough distance, it will 
eventually approach the speed of light and the nonrelativistic approximation 
will break down. Therefore our analysis is only relevant near the boundary. 
The solution may be made global by introducing a second boundary, such that 
the velocity is constant on the other side of the second boundary. We will see 
below that the stress tensor on the second boundary, to linear order in f , will 
be minus the stress tensor of the flrst boundary. Of course at higher orders one 
may expect an attraction between the two boundaries, and so this solution will 
not be stationary. The underlying assumption in this note is that the walls 
on the gravity side are built of a solid which remains flxed. The fact that the 
null energy condition is violated may be a sign that such a material would be 
inconsistent, as some null energy violating conflgurations lead to superluminal 
propagation or instabilities [57], although some do not [63], in which case this 
conflguration should only be considered over a sufliciently short timescale. It 
may be that this timescale is never suflicient for turbulence to develop. 

5.2 The Gravity Dual 

In Chapter 3 we gave a detailed presentation of the AdS/hydrodynamics cor- 
respondence, that yields a black brane metric dual to arbitrary flows in very 
particular conformal fluids, which for example obey the relations (15. ip . If the 
flow satisfles the hydrodynamic equations of motion (15. 3p . the dual satisfles the 
vacuum Einstein equations, in this case with a negative cosmological constant. 

5.2.1 A note on ultralocality 

As mentioned in the construction of the gravity dual to a generic conformal 
fluid, the correspondence, at least in the incarnation in Refs. [4] and 0, is 
ultralocal. It is useful at this point to explain a bit further what this means. 
Consider the set of ingoing null geodesies which run from the boundary to the 
black brane horizon. Clearly each point in the bulk is on precisely one such 
geodesic. Also each point on the boundary is on one such geodesic. Therefore 
these geodesies can be used to associate a flxed single boundary point to each 
bulk point. This association is not one to one, there is an entire geodesic worth 
of bulk points associated to each boundary point. 

To implement the map, one identifles the boundary with the Minkowski 
space on which the fluid lives. The metric and its derivatives at a point in the 
bulk are determined entirely by the fluid velocity and temperature and their 
derivatives at the associated boundary point. One does not need to know the 
behavior of the fluid elsewhere. This is the ultralocality of the correspondence. 
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In particular, the fact that the fluid satisfies the hydrodynamic equations of 
motion on the left and right (at x < and x > d) implies that the dual metric 
will satisfy the vacuum Einstein equations on the left and right, so long as 
the characteristic distance over which these quantities vary is greater than the 
inverse temperature. 

Our fluid does not satisfy the hydrodynamic equations of motion at < 
X < d. This means that the dual gravitational configuration will not satisfy 
the vacuum Einstein equations, instead it will only satisfy Einstein's equations 
with a nonzero stress tensor, which we will calculate in Sec. 15.31 

This ultralocality is somewhat different from the ultralocality that one 
encounters in classical field theories, or in the BKL limit of gravity theories, 
in that the bulk geometry is ultralocal in terms of null and not temporal 
evolution. This ingoing null identification identifies the temporal evolution of 
the fiuid with outward radial evolution for the gravity theory. That is to say, 
the metric at larger radii but the same time is determined by the fiuid in the 
future but at the same location. In particular, a timeslice of the bulk geometry 
is determined by the evolution of the boundary fiuid during a fixed interval of 
time. This interval is of order the inverse temperature, and so no appreciable 
evolution may occur during this interval if the temperature is large enough 
for the correspondence to hold. In this sense any fixed timeslice of the gravity 
dual contains only as much information as a fixed timeslice of the fiuid, despite 
being one dimension greater. 

As each event in the fiuid is identified with an inward null geodesic in the 
bulk, the metric corresponding to this event appears to be falling towards the 
black hole at the speed of light. This is not at all to say that there is a Killing 
vector in the inward null direction, the metric changes in that direction, but 
in a fashion which is fixed by the map. Thus a disturbance on the boundary 
creates gravity waves which fiy inward at the speed of light to the horizon. 
Similarly, pasting together an infinite sequence of bulk timeslices which are 
separated by time intervals 1/T, one obtains a pattern which falls from the 
boundary in to the horizon at the speed of light. Although each individual 
timeslice is too small to see any evolution, the entire pattern is dual to the 
entire history of the fiow. Like a movie reel, the pattern in turn allows one to 
reconstruct the gravity dual, as it contains the timeslices. 

One may use this identification to speculate on the gravitational dual of 
decaying turbulence. For example, the inverse cascade of decaying (2+1)- 
dimensional turbulence consists of a chaotic period during which the fiuid is 
subjected to random external forces followed by a relaxation period, charac- 
terized by the merging of well-separated vortices [6H [65]. This would then 
be dual to a kind of forest of gravity waves falling from the boundary to the 
horizon, beginning when the boundary is subjected to a random perturbation. 



106 



The Surface Layers Dual to Hydrodynamic Boundaries 



First tlie canopy, representing the chaotic period, falls out of the boundary 
into the horizon. When the external perturbation is turned off, it is followed 
by the branches representing the vortices, and then the branches merge into 
trunks as the vortices merge. The usual cascade [66l|67] in (3+l)-dimensional 
turbulence may be similarly described, but when the boundary is randomly 
perturbed, the trunk falls out first. This leads to the rather bizarre observa- 
tion that black brane geometries in AdS4 and AdSs respond very differently 
to random perturbations of their boundaries. Needless to say, it would be in- 
teresting to make this picture precise, or to see whether it is inconsistent with 
the various approximations involved in the duality. 



We will now calculate the metric dual to the flow fl5.2p using the map in Ref. [1] 
with the simplified notation of Refs. [SI El]. This map takes regions in which 
the flow satisfies the fluid equations to regions in which the metric satisfies the 
source-free Einstein equations. Acting on the region in which the fluid does 
not satisfy the hydrodynamic equation, the map is not known to have any 
special properties other than continuity, which will produce an interpolation 
between the vacuum Einstein metrics on the two sides. Therefore the choice of 
this map corresponds to a rather arbitrary choice of interpolation. However we 
will see that this interpolation has two nice properties. First, it is reasonably 
independent of the interpolating velocity function chosen. In particular, the 
third and higher derivatives of the velocity will yield contributions to the inte- 
grated stress tensor which are suppressed by powers of dT, while the leading 
contribution is independent of d. Second, the resulting stress tensor is simpler 
than the Israel stress tensor, it will have zero stress, whereas Israel's stress 
tensor has shear stress. 

Ultralocality in the ingoing null direction implies that the simplest co- 
ordinates in which to express the metric are Gaussian null coordinates, in 
which r parametrizes the ingoing null lines. In these coordinates, the bulk 
5-dimensional metric f l3.57p corresponding to an x-dependent 4-dimensional 
fluid flow is [5] 



rfs2 _ GMNdX^dX^ = -2u^{x) dx" {dr + V^(r, x) dx") + 6^^(r, x)dx^'dx' 



5.2.2 The metric 



(5.12) 



where, up to second derivatives in v, Vu and are defined as 



A 




(5.13) 
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and 



(g = r^P - 



Q2{br]— —P, 



d- 1 



(5.14) 



The functions that appear in this definition were introduced in Sec. 13. 7[ (one 
can refer also to the Appendix). 

Eq. (15.121) is the bulk metric dual to a fluid flow in an arbitrary curved 
space. The fluid is conformally invariant, and the conformal invariance has 
been used to write the metric in a compact form using objects which trans- 
form covariantly under the conformal symmetry. We are interested in a flat 
boundary, and so many of these objects will vanish. In fact, since u only de- 
pends on X, but only has nonvanishing components in the t and y directions, 
even the gauge fleld for a Weyl transformation will vanish 

1 



d 



-u 







(5.15) 



as u^Vxu^ = and V\u'^ = 0. This imphes that the Weyl-covariant derivative 
reduces to the ordinary derivative 



(5.16) 

The Weyl-covariant Schouten tensor S (deflned in eq. fl2.59p ) is propor- 
tional to the Weyl-covariant curvature of the boundary. As the Weyl-covariant 
derivative is just the ordinary derivative, this is just the ordinary curvature. 
As the boundary is Minkowski space, the curvature vanishes, and so the Weyl- 
covariant Schouten tensor also vanishes 







(5.17) 



Similarly the Weyl-covariant Weyl curvature C is the sum of the ordinary Weyl 
curvature and the curvature of the Weyl tensor (see eq. fl2.58p ). which both 
vanish and so 

C^uXa = . (5.18) 

The vorticity u does not vanish, however like the shear strain a it is of flrst 
order in v. Therefore w^, ua and cr^ terms are all of order 0{v'^) and so do 
not contribute at order 0{v). Thus only the third and fourth terms of f l5.13p 
contribute to V^. The third term is easily evaluated 

/ \ 





\ y 



(5.19) 
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Adding the third and fourth terms one finds at order 0{v) 

/ \ / 1 \ 



X)i{br) 





V / 



+ 





V 



- r 
2 



(5.20) 



The functions T, reported in Table 2.1, are easily expressed in terms of the 
shorthand notation <>, which symmetrizes and contracts with the projectors 
P^j,. The first two are identically zero, while the second two are of order ^) 

(5.21) 

(5.22) 

(5.23) 

(5.24) 







Therefore only the first and third terms of (15.141) contribute to (5f^u 
f V 0\ /0000\ 



(5.25) 



10 , 2 n ^ t;' 
^; 1 +^5i{br) 

\0001/ \0000/ 

Finally, inserting Eqs. f l5.20p and f l5.25p into fl5.12p . one finds the final form 
of the metric (TU] 



(l - ^) dt' - 2dtdr + (2^ + Mbry) dtdy 

,2/ T 2 , .2 , ,2\ , r,„2i 



+ r^{dx^ + dy^ + dz^) + 2r% Qi{br)v' dxdy - 2vdydr (5.26) 

where b = l/nT and the functions Di and Qi are defined in Eqs. fl5.69p and 
fl5.70p . Using the basis {t,x,y, z,r), we may write the metric in matrix form 



as 





^ + it)i(6r)t;" r%g,{br)v' 



\ 





-1 








^ + it)i(6r)t;" -1\ 


—V 

-V 0/ 



{,4^2 I 2 " 

r^60i(6r)f' 
,2 



while the inverse metric is 
/ 




9 



V 



1 







b 01 (br)v' 
-2 






b^x{br)v' 



1 









-1 



•Oi{br)v" 





Di(fer)i'" 



2r2 







+ v 
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5.2.3 Christoffel symbols 

In Sec. 15.31 we will see that the leading contribution to the stress tensor comes 
from the second derivative of the velocity. Contributions at that order come 
from the second derivative of the velocity in the curvature, which in turn 
contains contributions from first and second derivatives of the velocity in the 
Christoffel symbols, as well as from Christoffel symbols which are velocity inde- 
pendent as these are multiplied by velocity-dependent terms when calculating 
the curvature. We will now calculate all of the Christoffel symbols up to first 
order in v, v' and f", although the v terms will not contribute to the stress 
tensor. 

We will begin with the terms at order these are just the Christoffel 

symbols of the static black brane 



tt 






-pr 
tt 


— ^3 _ 


1 


pt 

XX 


= r* 

yy 


= r* = r 

^ ZZ ' 1 


tr 


pr 

~ rt 


= r + 


pa; 

xr 


-px 

rx 


1 

1 

r 


ry 

yr 


— ry 

ry 


1 

r 


zr 


= 

rz 


1 
r 


XX 


= r 
yy 


= r 

ZZ 



The new terms, at order 0{v), are 



ty '-yt ^"^U"w^ ^4^3) 



v 

ry =ry = 

^ tx ^ xt 2h^r^ 

y ^ y ^ l X)[{hr)hv" 

^ tr ^ rt A 

nr = nx = lb'g[ibr)v'-^^, 
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yy 



X 



■pr "pt 

yr ~ '- ry 



2r 

1 



ty yt y ^4^2 j \ ^ij. 



- lbt>[ibr)v"^ , 



v' 



^ ty — ^ yt — ~2bh^ 



r 



-^1 = 1 {v' + 2rbQ^ibr)v' + rH'Q[ibr)v') . (5.28) 



5.2.4 The Riemann tensor and the Ricci tensor and 
scalar 

Using the Christoffel symbols one can now easily compute the Riemann tensor. 
The order 0{v^) terms again are just those of the static AdS black brane 



(5.29) 

1 







^trtr 


= 1 3 


^txtx 


— Rtyty — Rtztz — 


Rfxrx 


— Rtyry Rtzrz 


^xyxy 


^xzxz Ryzyz 



(5.30) 



^^ + 1. (5.32) 



The bulk stress tensor is entirely determined by the contributions to the 
Riemann tensor which do not solve the fluid equations of motion, as it is these 
that do not solve the vacuum Einstein equations. If f is a constant, this yields 
the boosted black brane, which satisfies the vacuum Einstein equations. If v 
is linear, then again this is a solution of the linear order fluid equations as we 
have checked above, and therefore as we will check below v' will not contribute 
to the gravitational stress tensor at order 0{v). Therefore the flrst nontrivial 
contributions to the stress tensor arise from the Riemann tensor at linear order 
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in v" (v — v' — 0) 



Rtrty = + ^ ) ^^(br)v" (5.33) 



v" 



Rxrxy = ^ (2 + 2t)i(6r) + 2h\^ 0;(6r) - hrX>[{hr)) (5.34) 



v" 



Rzrzy = -J (2t)i(6r) - 6rD'i(6r)) (5.35) 

Ryxt. = -4^(2 + (b'r' - br)t,[(br)) (5.36) 

Ryrtr = ^U,--bX]- (5-37) 

As a check on our calculation, we will also calculate the contributions to 
the various tensors at linear order in the nondifferentiated velocity v 

Rtxyx = Rtzyz = ^ (^1 " ^i^^ V (5.38) 

Rtytr --(r' + T^]v (5.39) 



Rtryr = ( 1 " ) ^ (5-40) 



R-xyxr R-zyzr ~l~ ^4^,2^ ^ (^"^l) 



and in v' 



Rtxty = {b'r' + b'r' + (6V - l){2g,{br) + brQ[{br)))^ (5.42) 

i?t.,. = (-2 + + (6V^ + 6r)(20i(6r) + ftrflUftr)))^ (5.43) 

Rtyxr = -(2 + b^r' + (6^5 + 6r)(20i(6r) + brQ[{br)))^ (5.44) 
_ 2v' 

= -{b'r' + {b'^r^ - l)(20i(6r) + hr^,{hr)))^ (5.46) 

i?....--^6VW.(^r) + W;(6r)). (5.47) 

The Ricci tensor is now easily calculated. Again the order 0{v^) terms are 
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those of the static black brane solution 

= 4r^ - ^ (5.48) 

Rtr = Rrt = 4: (5.49) 
Rxx = Ryy = Rzz = '^T^ • (5.50) 

Contributions to the stress tensor will arise from the v" terms, {v = v' = 0) 

Rty = Ryt = -4^(2 + 4(1 + 6V)t)i(6r) + (6^5 - hr){>o\{hr) + hr>o'i{hr))) 

(5.51) 



v" 



R^y = Ry^. = —{2 + 4Di + hr{2br ^-^{br) - t>[{br) + hr-K)'[{hr)) . (5.52) 

The terms in the Ricci tensor proportional to v are those of a rigidly boosted 
black brane 

R'S = = 4- (5-53) 

which provides an exact solution both to the hydrodynamic equations and to 
Einstein's equations with a negative cosmological constant. Again the terms 
linear in v' yield a solution to the fluid equations and so Einstein's equations, 
although only to linear order 0{v) 



v' 



R = -{8b^r^Q,{br) + {bh^r^ - l)^^{hr))-— . (5.54) 
Using the large r asymptotic expansions of Ref. 

we flnd that the asymptotic behaviors of the v" terms in the Ricci tensor are 

13 v" 

K ~ ToV^ (5.57) 

1 v" 

The Ricci scalar is 

R=-20. (5.59) 

There is no contribution at order 0{v) to the Ricci scalar. This is guaranteed 
for any solution of the vacuum Einstein equations with cosmological constant 
A = —6, and so there could not have been any corrections from the v and v' 
terms. There are no corrections from the v" terms at linear order because the 
corresponding components of the inverse metric are themselves of order 0{v), 
and so the contributions to the Ricci tensor are of order (9(f^). 



Rry-^^T^ ■ (5.5^ 
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5.2.5 Contributions to the Ricci tensor at 0{v^^^) and 

(9(^(4)) 

Before continuing with the calculation of the bulk stress tensor, we will pause 
to discuss some of the approximations that we have made. We have made two 
truncations. First, we have calculated everything at order 0{v). As we are 
working in units in which c = 1, v is small for nonrelativistic speeds and so 
this is a valid approximation in a region in which the flow is sufficiently slow. 

A more dangerous truncation is that of higher derivatives of the velocity. 
The gravity/hydrodynamics correspondence is a one to one map between grav- 
itational and fluid solutions in a derivative expansion. More precisely, the kth 
order map relates the truncation of the fluid equations to k derivatives and 
that of the gravity equations to {k + 1) derivatives. The iterative procedure 
described in Ref. |1] in principle determines this map for all k, however in 
practice this map has only been determined to order k = 2. In other words, 
it provides a metric as a function of v, v' and f ", however a perfect matching 
with Einstein's equations would require also corrections involving the higher 
derivatives f which are not known. 

General arguments based on dimensional analysis suggest that these correc- 
tions become smaller at higher k. In general one expects that each derivative 
leads to a contribution which is subdominant by a factor of Tl with respect the 
previous derivative, where / is the distance scale of the derivative. Ideally one 
would like to check this claim for all terms with, say, three or four derivatives. 
However this would require a knowledge of the map at orders k = ?> and = 4. 

The map at order k = 2, which we have used, does produce some terms 
in the curvature which depend on the third and fourth derivatives of v. In 
this subsection we will verify that two of these have the expected convergence 
scaling, and determine the corresponding condition on our fluid flow. In other 
words, we determine a necessary condition for the derivative expansion to 
apply to our flow. 

The Ricci tensor components Rxy and Rty have corrections from the third 
and fourth derivatives of the velocity respectively 

^(3) _ t;(3)(a;) {brv[{br) + ^i{br)) 

<' - -HW^ , (5,61) 



We want to determine the condition under which R^^^ is subdominant to R^xy- 
As the higher derivatives of v define an interpolating function between two 
solutions over an interval of length d, each derivative is larger than the previous 
one by about 1/d. In other words, dx ^ 1/d. 
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To test the subdominance of Riy\ it is sufficient to compare it to the similar 
term in Rxy , which contains Di. The ratio of these terms is 

< . ^ . 1 (5,62) 

R'^y rv^^^x) rd 

therefore the fourth order term is subdominant if ^ 1/r in the entire bulk. 
The bulk extends from the horizon at r = 1/6 = ttT to the boundary at r = oo. 
Therefore convergence requires 

ci»^. (5.63) 

This fourth order term is suppressed by vrdT with respect to the third order 
term, in line with the above expectations from dimensional analysis. This 
means that the gravity duality procedure is only convergent when d is suf- 
ficiently large. Of course, the duality never yields a solution of the vacuum 
Einstein equations, and so one may argue that its convergence is immaterial. 
Nonetheless, it is only well-defined as a series when d satisfies (15.631) . 



5.2.6 The static black brane solution 



As a check on our calculation and conventions, we recover that the static 
{v = 0) black brane satisfies the vacuum Einstein equations with cosmological 
constant A = — 6 
















-6\ 







6r^ 



















6r^ 



















6r^ 
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-6 











/ 



(5.64) 



5.3 Two Calculations of the Stress Tensor 

In this section we will calculate the bulk stress tensor of the surface layer in- 
terpolating between the vacuum gravity solutions using two different methods, 
corresponding to two different metrics. First, we will apply the duality map 
of Ref. to a fluid flow which interpolates between the two solutions, the 
stationary solution on the left and the linear velocity solution on the right. 
In this case, as we have seen, the interpolating region is necessarily larger 
than the inverse temperature. Next, we will directly interpolate between the 
gravitational solutions using the Israel matching conditions [TT] . This method 
requires the interpolating region to be very thin, and uses the assumption that 
in this limit the extrinsic curvature remains bounded. 
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5.3.1 Interpolating between the hydrodynamic flows 

The duality map of Ref. |1] takes a fluid flow and yields a dual metric. This 
dual metric solves the vacuum Einstein equations when the fluid flow satisfies 
the hydrodynamic equations of motion (15. 3p . If the flow does not satisfy the 
equations of motion, the dual metric does not satisfy the vacuum Einstein 
equations. Thus apparently there is no benefit in using this map over any 
other map. However we will use the map, and observe the consequences. The 
resulting dual metric will necessarily solve Einstein's equations with some value 
of the stress tensor 

SttGatT^i, = R^^ - ]^Rg^u + ^QtJiu ■ (5.65) 

We will determine this value. 

We saw in Eq. (I5.64p that there is no contribution to the stress tensor at 
order C(t>°). We have argued that, at order C(f ), the dominant contributions 
to the stress tensor are proportional to v" . These are easily found from fl5.65p 
to be 

_ v"{x) (4 (6V - 1) Oi(6r) - hr {b^r^ - 1) {v[{br) + 6rt)'/(6r)) - 2) 



(5.66) 

_ v"{x) (4t)i(6r) + br {2brg[{br) - tj[{br) - hr>o'{{hr)) + 2) 



There appears to also be a contribution proportional to v' 

^ _ v\x) (br- {(b^r^ - 1) Q'jibr) + 3&r) + (5&V - 1) g^&r)) 

At order v' one expects no contributions to the stress tensor, as a solution with 
a linear velocity satisfies the fluid equations at order 0{v). Therefore a non- 
trivial contribution would be in contradiction with the gravity/hydrodynamics 
correspondence. We will see shorty that this contribution is in fact equal to 
zero. 

The functions t)i(r) and 0i(r) are defined as 

Oi(r) = -/ dxx' dy—^—— (5.69) 



x^ 



01 W= / dx (5.70) 

X{X^ — 1) 
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Integrating we [5] obtain analytical expressions for t)i(r) and Qi{r) 

«i = "i + i + ^(-' - 1) (l°g 1^ + 2 tan-^(r) - vr) (5.71) 
1/ /(l + r)2(l + r2) 



Si = - (^log (^^ '-^^ '-j - 2tan-^(r) + nj . (5.72) 

The derivatives of these expressions are 

// ^'^ ^ 111 

01 = "(r2 + l)2 + (H + 1)2 + 2(r2 + l) + H " 2(r + l)2 ^^"^^^ 

and 

= (vrr^ + 2 (r^ + l) log(r + 1) - 2 (r^ + l) tan-i(r) 

_4^3 ^ - (r^ + l) log (r^ + l) + tt) (5.75) 



+ (r(r + 1) (r'l + r 2 - 3) - 3) r (2 log(r + 1) 
- log (r^ + l)) - 2(r + 1) (r^ + l) (r^ - 3) tan i(r) 

-2 (2r^ + + + r + 3) - 6 log(r + 1)) . (5.76) 

The explicit formula Eqs. (I5.73p and ( I5.74p for the derivatives of Qi can be 
combined to show that 

r ((r^ - 1) g'lir) + 3r) + (5/ - l) g[{r) = . (5.77) 

This combination is proportional to formula f l5.68p for T^y, therefore 

T^y = (5.78) 

and there are no contributions proportional to v'. 

Similarly one may evaluate the combination of functions that appears in 



T 

±ry 



r {2rg[{r) - r<(r) - \j[{r)) + 40i(r) + 2 = 0. (5.79) 
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This implies that 

Try = (5.80) 

leaving only Tty, the momentum in the y direction. Thus the bulk stress tensor 
contains no stress, only momentum. 

We may use the exact expressions for the functions Qi and Oi to simplify 
the only nonvanishing component of the stress tensor [10] 



v"(x] 



Using the fundamental theorem of calculus, this may be integrated over the 
interpolating region to obtain 

■ v' 

where v' is the derivative of the velocity in the region x > d. In particular, at 
this leading order the integrated stress tensor of the surface layer is independent 
of the interpolation and independent of d. Of course it still depends on the 
map that we used to generate the dual metric. 

Had the v' term been the dominant contribution, the stress tensor would 
have been constant, and so the integral would be have proportional to d. 
Similarly a v^^^ term would have led to a stress tensor proportional to 1/d, and 
higher powers of v to other scalings. Therefore it is somewhat nontrivial that 
the leading contribution to the integrated stress tensor is in fact d-independent. 
Clearly this li-independence is desirable, as d is not a physical quantity but 
merely an artifact of the scheme that we used to regularize the divergent second 
derivative of the fluid velocity. 

The bulk stress tensor does not satisfy any of energy conditions we discussed 
in Sec. 14. 1^ not even the null energy condition. As the only nonvanishing 
component is Tfy, the only nonvanishing product of a null vector w and the 
stress tensor is 

w^Tw = 2w%ywy. (5.83) 

As Tfy is already of order 0{v), at order 0{v) one need only consider the terms 
in w of order 0{v^). That is to say, w only needs to be null with respect to 
the static black brane metric. Consider for example the null vectors w± 



wi = r, wl = ±r\ll-^^. (5.84) 



The product (I5.83P is 
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which is nonzero. However and w_ yield opposite signs, as incidentally 
do the two choices of signs of v. Therefore at least one of these will yield 
a negative product, and so the bulk stress tensor does not satisfy the null 
energy condition. This may or may not mean that no external matter may be 
consistently added which produces such a surface layer. 



5.3.2 Israel's matching conditions on the gravity duals 

We will now calculate the bulk stress tensor in a different geometry. Following 
Ref. [11] and discussion in Sec. 14. 4[ we will consider the vacuum Einstein 
solution corresponding to a static fluid on the left and that corresponding to 
a linear velocity flow on the right. These solutions will be glued together 
by interpolating continuously between the two metrics over a distance d and 
taking the limit d — )■ such that the extrinsic curvature remains bounded. 
In Ref. [TT], Israel has shown that the resulting conflguration contains two 
solutions separated by a surface layer whose bulk stress tensor is independent 
of the interpolation used. 

Following Ref. [TT], the flrst step in the calculation of the stress tensor is 
the deflnition of the unit normal vector to the hyperplane 

n^ = {0,r,0,0,0} , (5.86) 

which satisfles the normalization condition 

Uf^g^-'n, = ^{n,f = 1 . (5.87) 

The surface layer E extends along all of the directions except for the x direction. 
A basis of tangent vectors to E is 

ds = e(i)dx' (5.88) 

where 



e(i) = {1,0,0,0,0} (5.89) 

e(,) = {0,0,1,0,0} (5.90) 

e(,) = {0,0, 0,1,0} (5.91) 

e(^) = {0,0,0,0,1} . (5.92) 

In terms of these tangent vectors the extrinsic curvature may be calculated as 
K., = eo) ■ V,n = ^ - n^r^,, = ^ - n^F^ . (5.93) 



5.3. Two Calculations of the Stress Tensor 



119 



On the left, where the fluid is static (f = 0), substituting (I5.27P into f l5.93p 
one finds no extrinsic curvature 

k[-^ = -rVl = (5.94) 
K^-^ = -rr^, = . (5.95) 

On the right, where the fluid velocity is linear, the Christoffel symbols of 
Eq. fl5.28p yield a nontrivial extrinsic curvature K^^\ 

< = -'^l = ^ (5-96) 
i^W = -rr^, = 1^ (1 + b'r'Q[{br)) . (5.97) 

The tensor fl4.48p is defined to be the difference between the extrinsic 
curvatures on the two sides of the surface layer 

7., = Ki;^-Ki-K (5.98) 

The bulk stress tensor integrated over x is equal to the tensor Sij, defined by 

- SnG^S,, = 7., - 9^J7Z ■ (5.99) 

The expression fl5.99p for the integrated bulk stress tensor was derived in 
[TT] for a 4- dimensional space with no cosmological constant. While several fac- 
tors in the derivation change in our current 5-dimensional situation, Eq. fl5.99p 
remains unchanged. The cosmological constant term yields a contribution pro- 
portional to the integral of A times the metric integrated over the thickness d 
of the surface layer. As the metric is taken to be finite, this term vanishes in 
the (i — )■ limit. 

The trace of 7 is O(f^), therefore fl5.99p yields the integrated bulk stress 
tensor 

These are equal to the integrals over the x direction^ of the stress tensors T^^^ 
of Subsec. 15.3.11 at order /c = 1, in other words, without the Oi term that 



^Note that, following Ref. [TT], the measure of this integral must be that of x rescaled to 
normal coordinates. Therefore the integral contains an additional factor of r = y/g^ 



ixx • 
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entered into tlie metric fl5.2.2p multiplied by v" 



v 



,11 



^(1) 

ty 



(5.102) 



1) 



(l + 6Vg;(6r 



)) 



(5.103) 



yr 



The Oi terms arose from the dualization of the interpolating region, which did 
not satisfy the equations of motion. It therefore cannot enter into the Israel 
calculation, which uses only the solutions of the vacuum Einstein equations. 
Indeed, the Oi terms in fl5.2.2p are singular in the limit (i — )■ as v" diverges 
as and therefore the boundedness of the extrinsic curvature assumed in 
Israel's derivation fails for the metric interpolation fl5.2.2p . 

Like the stress tensor fl5.8ip calculated by interpolating the hydrodynamic 
flow, the Israel stress tensor does not satisfy the null energy conditions. Again, 
to linear order in f , one may consider vectors which are null with respect to the 
static black brane metric. Therefore, again one may consider the null vectors 
w± of Eq. fl5.84p . As Tty is, at least for any flnite (i, equal to that of Subsec. 15.3.11 
divided by the positive combination 6r, the sign of the inner product (15.851) is 
unchanged. Therefore the null energy condition is also violated by this stress 
tensor. 

The main difference between the two stress tensors is then that Ty,. does 
not vanish for the Israel tensor. Remembering that in our Gaussian null coor- 
dinates the r direction is the sum of a spatial and temporal piece, the spatial 
component implies that there is a nonzero stress. More precisely, while both 
Israel's thin surface layer and the thick fluid surface layer have a nonvanishing 
y momentum, the Israel surface layer also has a flux of this y momentum in 
the radial direction, from the boundary into the horizon of the black brane. 
As the black brane is inflnite in the x direction, this is not problematic for the 
time-independence of the solution. 



Conclusions &; future directions 



Turbulence often arises as a result of the boundary conditions placed on a fluid. 
As a preliminary step towards an understanding of turbulence in gravity, we 
have proposed two gravitational duals of such boundaries. Both of these duals 
involve the addition of a surface layer of matter, with a certain stress tensor. 
These proposals are in a sense trivial, as the dynamics of the duals is defined 
not by any known equations of motion, but by the duality map itself. It remains 
to be shown whether such matter can exist. For example, even if the equations 
of motion which it obeys can be found, the existence of a UV completion of 
the matter theory may be fundamentally obstructed as in Ref. [68]. Or the 
failure of the null energy condition may imply that, whatever the ultraviolet 
theory may be, the wall simply disintegrates before it has any significant effect 
on the fluid. 

Of course, an ultraviolet completion is not necessarily a prerequisite for 
learning something interesting about whatever the gravitational dual to tur- 
bulence may be. After all, no ultraviolet completion of Einstein gravity is used 
in this correspondence. The surface layer implies the existence of equations 
of motion which are distinct from the Einstein vacuum equations and perhaps 
pathological. However the interesting part of the fluid, the turbulent part, is 
not at the wall. For example, if we consider the motion of a fluid in a pipe, the 
flow may be turbulent throughout the interior of the pipe. The ultralocality 
of the duality map implies that, at a distance greater than 1/T from the pipe, 
the vacuum Einstein equations are still satisfied by the gravity dual. Thus in 
a sense the ultralocality decouples the problem of understanding turbulence in 
gravity from the problem of defining a gravity dual of a boundary. 

Besides trying to characterize the gravitational dual of turbulent flow, the 
other interesting question is to find the gravitational dual of the conditions 
under which turbulence can occur. In nonrelativistic, incompressible flows, 
turbulence is expected when the product of a system's characteristic scale L 
times the characteristic velocity w of a fluid is much greater than the kine- 
matic viscosity. In Ref. [9], the authors claim that for the conformal fluids 
dual to AdS black branes, turbulence is expected when LTv ^ 1, where T 
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is tlie temperature of tlie fluid. The AdS/hydrodynamics correspondence is 
expected to be reliable at scales L such that LT ^ 1. Therefore since w < 1, 
it appears that whenever turbulence is expected, LT > LTv ^ 1 and so the 
correspondence can be trusted at least for quantities that vary over a distance 
L. (3 + l)-dimensional turbulence is characterized by vortices of various sizes 
from L down to the dissipation scale [S]. Thus the duality appears to be reli- 
able at least for the largest vortices in a turbulent flow. The dissipation scale 
is a function of L, T and v, and so in principle one may determine whether or 
not the duality is reliable for vortices all of the way down to this scale and so 
for the entire flow. 

Understanding the gravity duals of turbulent flows, as described above, 
may yield new insights into the dynamics of black branes in AdS space, per- 
haps revealing a surprising difference between branes in AdS4 and AdSs, or 
indicating that generically they come with funnels attached as in Refs. [69]. 
The main weakness of this program is the dependence on asymptotically AdS 
geometry in the duality map of Ref . [1] . There was no such restriction in the 
original correspondence of Ref. [2] , nor in other identifications of black holes 
and viscous fluids such as the blackfold program of Refs. [TOl [71] and the Wilso- 
nian identiflcation of Ref. [72]. An extension of turbulence to asymptotically 
Minkowski space could relate (3-1-1 )-dimensional fluid dynamics to wealth of 
studies of asymptotically Minkowski 5d black objects, such as Refs. [73]. More 
importantly, relaxing the asymptotically AdS condition may mean that fluid 
mechanics, perhaps in only 2+1 dimensions, has something to teach us about 
real world gravity. 



Appendix A 

Conventions &; Notation 



We work in the ( — h + • • •) signature, /i, v denote space-time indices, j = 
\ . . .k label the k different conserved charges. The dimension of the spacetime 
in which the conformal fluid lives is denoted by d. In the context of AdS/CFT, 
the dual AdS^+i space has d + 1 spacetime dimensions. 

We adopt standard symmetrization and anti-symmetrization conventions. 
For any tensor Fah we define the symmetric part as 

F(afe) = \ {Fab + Fba) (A.l) 

and the anti-symmetric part as 

F[ab] = - {Fab — Fba) ■ (A. 2) 

We also use ^ to indicate the velocity projected covariant derivative: 

^ = M^V^. (A.3) 

For any two tensor T'"^ we denote the symmetric traceless projections trans- 
verse to the velocity field as: 

^it^u) ^ pi^a pup ^^^^ _ pM. p«/3 _ 4^ 

(Jj -L 

A.l Integral functions 

We list in the following the two different conventions on integral functions 
entering in the second order metric dual to fluid solutions on the boundary. 
The definitions adopted in many works {e.g. /^O ^^^^ 
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F{br) = 
Hiihr) = 
H2{hr) = 



OO yd-1 _ I 



br 



y(^yd _ 


1) 


yd-2_ 


- 1 


y{jjd _ 


1) 








1) 



dy, 



dy, 



(A.5) 
(A.6) 
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br 



y'-^dy [l + {d- l)yF{y) + 2y^F'{y)\ (A.7) 
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+ (2{d - Ijf" -(d- 2)) / dy y'F'iy) 
Definitions adopted in [31] are 
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(A.9) 



d-l 



(A.IO) 
(A.ll) 



Qiibr) 



oo „,d~l 



y"-' - 1 
y{y^-l) 



dy, 



Q2{hr)^2{rhf j §1 dyy'F\yf 

'br S 



(A.12) 
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and tlierefore one can easily move across definitions noticing that 



A.2. List of symbols adopted 
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01 (6r) = F{br), 
03(6r) = {hrfHi{hr), 

A. 2 List of symbols adopted 

We summarise in table A.l some of the symbols we adopted. The relevant 
equations defining symbols are denoted by their respective equation numbers 
appearing inside parenthesis. 



02 (or) 
X)2{hr) 



2{rbyKi{br 
1 



2{hr) 



d-2 



K^ihr) 



(A.18) 



Symbol 


Definition 


Symbol 


Definition 


d 


dimensions of spacetime 


P 


Pressure 


s 


Proper entropy density 


P 


Energy density 


T 


Fluid temperature 




Chemical potentials of the fluid 


c 


Bulk viscosity 


ry 


Shear viscosity measured at 


qi 


charge density 




zero shear and vorticity 


11 


Contr. en. density to charge curr. 




Pseudo- vector transport coeff. 




Matrix of charge diffusion coeff. 








Energy-momentum tensor 




Entropy current 




Charge currents 




Fluid velocity (u^u^ — ^1) 




Spacetime metric 




Projection tensor, g^"^ + u^^w" 




Fluid acceleration (12.151). (12.401) 
Shear strain rate (15.71) and (12.401) 
See (|2.23() and (|2.40p 


d 


Fluid exoansion (|2.14p. (|2.401) 
Fluid vorticitv (|2.171). (|2.40p 




Wevl-cov. deriv. (12.43P andl2.44p 
Lorentz-covariant derivative 


r ^ 


See (|2.48|) 

Christoffel connection 




Riemann Curvature 




See(|2.52|) 




Ricci tensor/scalar 
Einstein tensor 
Schouten tensor (]2.59p 
See (|2.59|) 


n 


See (|2.56|) 
See (I2.5fi|) 

Wevl Curvature (12.581) 
See (|2.58|) 
See Table 2.1 



Table A.l: List of principal symbols & definitions adopted. 
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